Introducción 
a las ecuaciones 
diferenciales ordinarias 


Leónidas Cerda Romero 
Janneth Morocho Yaucán 



INTRODUCCIÓN A LAS ECUACIONES 
DIFERENCIALES ORDINARIAS 



INTRODUCCIÓN A LAS ECUACIONES 
DIFERENCIALES ORDINARIAS 


Leónidas Cerda Romero 
Janneth Morocho Yaucán 



DIRECCIÓN DE 
PUBLICACIONES 











INTRODUCCIÓN ALAS ECUACIONES 
DIFERENCIALES ORDINARIAS 


© 2018 Leónidas Cerda Romero y Janneth Morocho Yaucán 


© 2018 Escuela Superior Politécnica de Chimborazo 

Panamericana Sur, kilómetro 1 V 2 
Dirección de Publicaciones Científicas 
Riobamba, Ecuador 
Teléfono: (593 3) 299 8200 
Código Postal: EC060155 

Aval ESPOCH 

Este libro se sometió a arbitraje bajo el sistema de doble ciego 
{peer review). 

Corrección y diseño: 

La Caracola Editores 

Impreso en Ecuador 

Prohibida la reproducción de este libro, por cualquier medio, sin la previa 
autorización por escrito de los propietarios del Copyright. 


CDU: 517.91 

Riobamba: Escuela Superior Politécnica de Chimborazo 
Dirección de Publicaciones, año 2017 
173 pp. vol: 17 x 24 cm 
ISBN: 978-9942-35-642-0 

1. Análisis matemático 

2. Cálculo diferencial 

3. Cálculo integral 

4. Ecuaciones diferenciales ordinarias 


3 




Contenido General 


1 Ecuaciones diferenciales de primer orden 7 

1.1 Modelos matemáticos. 7 

1.1.1 Ejemplos de modelos matemáticos. 10 

1.2 Ecuaciones diferenciales ordinarias de primer orden separables .... 28 

1.3 Ecuaciones diferenciales ordinarias que se pueden reducir a ecuaciones 

separables. 36 

1.4 Ecuaciones diferenciales exactas y factor de integración.42 

1.5 Ejercicios del Capítulo 1. 56 

2 Ecuaciones Diferenciales lineales de orden superior 60 

2.1 Teoría básica de las ecuaciones diferenciales lineales. 60 

2.2 Ecuaciones lineales homogéneas con coeficientes constantes. 76 

2.3 Ecuaciones no homogéneas con coeficientes constantes. 84 

2.3.1 Método de los coeficientes indeterminados. 85 

2.3.2 Método de variación de parámetros. 92 

2.4 Ecuaciones lineales homogéneas con coeficientes variables.100 

2.4.1 Solución alrededor de puntos ordinarios.102 

2.4.2 Solución alrededor de puntos singulares.109 

2.5 Aplicaciones de las ecuaciones diferenciales lineales de orden superior 116 

2.6 Ejercicios del capítulo 2 123 

3 Sistemas de ecuaciones diferenciales lineales 125 

3.1 Generalidades sobre los sistemas lineales.127 

3.2 El operador diferencial.131 

3.3 Método matricial para sistemas normales homogéneos con coeficientes 

constantes.137 

Bibliografía .155 


4 




















Indice de figuras 


1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 

1.10 

1.11 

2.1 

2.2 

2.3 

2.4 

2.5 

2.6 
2.7 


dP 

Función qn e satisface la ecuación — = kP . 

Curvas ortogonales en su punto de intersección. 

Familia xy — c\ = 0 y sus curvas ortogonales. 

Representación de la entrada-salida de un medicamento en un órgano. 

Representación de un tanque cilindrico de almacenamiento de agua . 

Representación de una curva de persecución. 

Aproximación de la rapidez de un objeto. 

Tanques de mezclado conectados . 

Gráfico de la función L(t) para algunos valores de k . 

dy 

Curvas integrales de la ecuación x + y— = 0. 

dx 

Solución al problema — = f(x,y),y(x o) = yo . 

Resorte suspendido desde un techo, longitud natural L . 

Masa en equilibrio, longitud del resorte L + l . 

Masa a una distancia x por debajo de la posición de equilibrio; longitud 

del resorte L + l + x . 

Circuito en serie. 

Viga horizontal. 

Aplicación de una carga a una viga. 

Momento flexionante en a:. 


10 

10 

12 

16 

19 

21 

22 

24 

26 

31 

32 

117 

117 

118 
120 
122 
122 
123 


5 


















Prólogo 


El presente texto ha sido elaborado pensando en entregar material para un curso 
introductorio de ecuaciones diferenciales ordinarias. En éste se tratan los tipos más 
comunes de ecuaciones diferenciales, así como sus métodos de resolución. 

El texto ha sido diseñado de tal forma que los lectores tengan primero una 
panorámica general de los temas tratados. Se han incluido las demostraciones de 
algunos resultados que son propias de los autores pero, se las puede hallar en la 
mayoría de textos introductorios a las ecuaciones diferenciales. Se han realizado 
ejemplos detallados para un entendimiento adecuado de los procesos que están 
involucrados en la resolución de cada tipo de ecuación. 

La obra ha sido desarrollada en tres capítulos, cada uno de los cuales comienza 
con una visión general de los temas tratados en el transcurso de los mismos. El 
capítulo 1 realiza un tratamiento bastante completo de las ecuaciones diferenciales de 
primer orden. El capítulo 2 trata las ecuaciones diferenciales de orden superior, sus 
métodos resolutivos y algunas aplicaciones de las ecuaciones diferenciales de segundo 
orden. El capítulo 3 está dedicado al estudio de los sistemas lineales de ecuaciones 
diferenciales. 

Se piensa tratar temas referentes a la transformada de Laplace y su aplicación a 
la resolución de problemas con condiciones iniciales, así como la teoría de las series 
de Fourier en un segundo tomo. 
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Capítulo 1 


Ecuaciones diferenciales de primer orden 

La teoría de las ecuaciones diferenciales es una aplicación directa del cálculo diferencial 
e integral. El cálculo diferencial es una herramienta poderosa para la construcción 
de modelos matemáticos de problemas en los cuales están involucrados la variación 
de una o varias variables, las variables dependientes, con respecto a otra u otras 
variables, las variables independientes. 

Determinar el crecimiento poblacional, calcular el decaimiento radioactivo, hallar 
la edad de un fósil, son problemas en los cuales se utiliza el cálculo diferencial para 
el planteamiento de un modelo matemático. 

El modelo matemático que se construye para hallar la solución de un problema 
particular deja planteado un nuevo problema: ¿qué herramientas se deben utilizar para 
hallar la solución del modelo matemático? El cálculo integral es quizá la herramienta 
matemática más usada en la resolución de modelos matemáticos que describen la 
variación de ciertas cantidades con respecto a otras. 

Este capítulo esta dividido en cuatro secciones. La sección 1.1 muestra algunos 
ejemplos de modelos matemáticos que describen la variación de una variable con 
respecto a otra. En la sección 1.2 se defininen conceptos básicos relacionados con 
las ecuaciones diferenciales, se realiza una clasificación de estas y se estudian 
las ecuaciones diferenciales separables. La sección 1.3 está dedicada al estudio 
de ecuaciones diferenciales que se reducen a ecuaciones separables por medio de 
sustituciones. Finalmente, la sección 1.4 estudia las ecuaciones diferenciales exactas 
y aquellas que se hacen exactas al multiplicarlas por un factor adecuado, llamado 
factor de integración. 


1.1 Modelos matemáticos 

La parte más complicada al usar la matemática para estudiar una aplicación particular 
es la conversión de fenómenos reales al formalismo matemático. Esto presupone 
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tácitamente la conversión de hipótesis imprecisas (guiadas por la creencia de que 
cierto fenómeno real tiene un determinado comportamiento) en fórmulas muy precisas 
(normalmente fórmulas matemáticas). 

Un modelo matemático es una representación, por medio de fórmulas matemáticas, 
de algo real. El objetivo no es tener una copia exacta de lo real, sino más bien 
representar las carcaterísticas que a uno le interesa estudiar de la cosa real. La 
construcción de modelos matemáticos siempre es una tarea complicada. Blanchard 
|BP, pag. 16-17] sugiere que los pasos básicos en la elaboración de un modelo 
matemático deben ser los siguientes: 

Paso 1. Establecer claramente las hipóteisis en que se basará el modelo. Estas deben 
describir las relaciones entre las cantidades que se van a estudiar. 

Paso 2. Definir completamente las variabes y parámetros que se usarán en el modelo. 

Paso 3. Usar las hipótesis formuladas en el paso 1 para obtener ecuaciones que 
relacionen las cantidades del paso 2. 

Observación 1.1. Cuanto más son las hipótesis involucradas en el paso 1 mayor es 
la precisión del modelo matemático donde “precisión del modelo” se dehe entender de 
la siguiente manera: supongamos que tenemos dos modelos Mi y M 2 de A. Decimos 
que Mi es más preciso que M 2 cuando M\ describe mejor el comportamiento de A. 
Esto no quiere decir que M 2 no sea útil, su utilidad dependerá del interés que tenga 
quien construye el modelo. 

Un ejemplo bastante sencillo de construcción de un modelo matemático es el 
crecimiento ilimitado de la población. La hipótesis en la que se basa este modelo es 
la siguiente: 

“La velocidad de crecimiento de la problación es proporcional al tamaño 
de la población”. 

Notamos que esta hipótesis no toma en consideración las limitaciones de espacio o de 
recursos; sin embargo, es bastante razonable para poblaciones pequeñas en entornos 
grandes, por ejemplo los primeros brotes de moho en un pan. Hasta aquí liemos 
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completado el paso 1 de las indicaciones dadas anteriormente por Blanchard. El paso 
2 pide determinar las variables y parámetros que serán utilizados en el modelo, en 
este ejemplo utilizamos las siguientes variables y parámetros: 


• t representará el tiempo (variable independiente). 


• P representará el tamaño de la población (variable dependiente). 


• k representará la constante de proporcionalidad (parámetro) entre la tasa de 
crecimiento de la población y el tamaño de esta. 


Finalmente, el paso 3 pide usar las hipótesis del paso 1 para establecer relaciones 
entre las cantidades descritas en el paso 2. Utilizando el cálculo diferencial es muy 
fácil darse cuenta de que: La única relación de ligadura entre las cantidades descritas 
anteriormente es: 


dP 

dt 


kP. 


Además podemos pedir que P(to ) = Pq lo cual nos indica que la población en 
el tiempo t = t 0 es igual a P 0 . Luego el modelo pide resolver el siguiente problema 
matemático: 


d Z = kP 

dt 

P(to) = Po. 


En la sección 1.2 llamaremos a problema de este tipo “problema con condiciones 
iniciales”. 

La solución de este problema es una función P que satisface tanto la ecuación 
dP 

= kP como la condición P(t 0 ) = P 0 . 

La siguiente figura muestra una posible solución a este problema. 
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Figura 1.1: Función qu e satisface la ecuación — = kP 

dt 

1.1.1 Ejemplos de modelos matemáticos 

A continuación se presentan algunos ejemplos de construcción de modelos. Donde es 
posible, se presenta la solución al problema matemático que genera el modelo. 

Trayectorias ortogonales 

Desde la geometría analítica sabemos que dos curvas C\ y C 2 son ortogonales en su 
punto de intersección si y solo si sus tangentes Tí y 71 son perpendiculares en este 
punto (ver figura 1.2). 



Figura 1.2: Curvas ortogonales en su punto de intersección 
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La noción de ortogonalidad se puede extender a dos familias de curvas de la 
siguiente manera: 

Definición 1.1. Sean G(x,y,ci) y H(x,y,C 2 ) dos familias de curvas. Se dice que 
las familias son, cada una, trayectorias ortogonales de la otra cuando todas las 
curvas de G(x,y,Ci) cortan ortogonalmente a todas las curvas de H(x,y,C 2 ). 

Las trayectorias ortogonales aparecen en aplicaciones de electricidad y 
magnetismo. Por ejemplo en el campo eléctrico que rodea a dos cuerpos de carga 
opuesta, las líneas de fuerza son perpendiculares a las curvas equipotenciales. 

Un problema típico sobre trayectorias ortogonales consiste en lo siguiente: Dada 
la familia G(x, y , cf) hallar la familia H(x, y , cf) tal que G(x, y , cf) y H(x, y , C 2 ) son 
trayectorias ortogonales. 

Para hallar el modelo matemático que permite hallar la familia H(x,y,C 2 ) nos 
damos cuenta de que la única hipótesis con que contamos es la ortogonalidad de 
las familias G(x,y,c\) y H(x,y,C 2 ), este sería el paso 1 descrito anteriormente. Las 
variables que intervienen en este problema son x, y que representan las coordenadas 
que describen las curvas. Este problema no tiene parámetros (ci y C 2 son parámetros 
de las curvas, no parámetros del modelo), este es el paso 2. Finalmente, para el paso 
3, procedemos como sigue: 

i. Hallar las tangentes a la curva G(x, y , cf). Para hallar estas tangentes calculamos 
dy 

— = f(x,y), derivando implícitamente G(x,y,Ci) en caso de ser necesario. 
dx 

ii. Exigimos que las tangentes de la familia H(x, y , C 2 ) sean perpendiculares a 

= f(x,y ). Esto se consigue poniendo = — —- r. 

dx dx j (x, y) 

Por tanto, el modelo de este problema es: 

dy = 1 

dx f(x, y)' 

Para ejemplificar lo descrito arriba nos permitimos presentar el siguiente ejemplo: 

Ejemplo 1.1. Hallar las familia de curvas ortogonales a la familia de hipérbolas 
xy — C\ = 0. 
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dy 

Solución. Despejamos la variable y de la expresión xy — C\ = 0 y hallamos —. 

dx 

Estos cálculos nos dan como resultado — = — —jr. Luego, reemplazando el valor de 

ax x ¿ 

Ci, tenemos 


dy_ _ _y 

dx x 


Por tanto, el modelo que describe la familia de curvas ortogonales a la familia de 

hipérbolas es -j- = —. 

dx y 

Para hallar la solución al problema matemático que deja planteado el modelo, 
procedemos de la siguiente manera: escribimos la última ecuación (el modelo 
matemático) en la forma ydy = xdx ; integrando cada lado de esta igualdad se 
tiene y 2 = x 2 + C 2 donde C 2 es una constante. Luego la familia de curvas ortogonales 
a la famila de hipérbolas xy — c\ = 0 es la familia de hipérbolas y 2 — x 2 — = 0. 

El siguiente gráfico muestra algunas curvas de estas familias de curvas ortogonales. 



Figura 1.3: Familia xy — ci = 0 y sus curvas ortogonales 


Determinación de edades 

La desintegración radiactiva se entiende como la disminución (variación) con el paso 
del tiempo de la intensidad de la radiación de cualquier material radiactivo. 
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A partir de datos experimentales, se ha llegado a concluir que si N es la cantidad 
de un material radiactivo (intensidad de radiación), entonces la desintegración de N 
es proporcional a N ; luego, el modelo que expresa esta disminución está dada por 


dN 

dt 


= -XN, 


donde el signo ” indica disminución y A es la constante de desintegración radiactiva. 

El carbono 14 (C 14) es un isótopo radiactivo del carbono que permite estimar 
la edad de fósiles y otras materias orgánicas a partir del modelo de desintegración 
radioactiva y del conocimiento que la semivida 1 de C 14 es de 5600 años. Este método 
de determinación de edades ha contribuido grandemente al conocimiento que tenemos 
del tiempo pasado. Su empleo se ha generalizado habiendo alcanzado un alto grado 
de precisión en muestras de edad conocida, demostrando ser un excelente método de 
determinación cronológica. 

Como ejemplo de la utilidad de este método de determinación de edades nos 
permitimos presentar el siguiente ejemplo: 


Ejemplo 1.2. Se encuentra que un hueso fosilizado contiene la centésima parte de 
la cantidad de C 14 encontrada en la materia viva. Determine la edad del fósil. 


Solución. Sea N(t) la cantidad de C 14 presente en el hueso al tiempo t y Nq la 
cantidad inicial de C 14. Se tiene que el modelo matemático que describe la variación 
de la cantidad de C 14 en el hueso es: 


dN 

— = -XN 
dt 

N( 0) = N 0 . 


dN 

Desde la primera ecuación se tiene que = —A dt. Luego, recordando que la 


semivida del C 14 es 5600 años, tenemos que: 


'N 0 


No/2 dN _ 

w~~ x 


-5600 


dt. 


0 


: La semivida o edad media de un material radioactivo es el tiempo que tarda este material en 
desintegrarse a la mitad de su intensidad original 
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Realizando los cálculos (hacerlo) nos damos cuenta de que A 


ln(l/2) 

5600 


Nuevamente, desde la primera ecuación, tenemos que: 

[No/ioo dN ln(l/2) rt 


IN 0 


N 


5600 


dt, 


donde t es el tiempo en el cual la cantidad de CIA es una centésima parte de la 
cantidad original presente en el hueso. Haciendo cuentas se tiene que t « —37206 
años. 


Observamos que nuestros cálculos nos conducen a un tiempo negativo. Este 
resultado se debe interpretar de la siguiente manera: El valor “í —37206” es el 
tiempo que debemos recorrer hacia atrás para determinar la edad del fósil. Por 
ejemplo, si t es el año 2000, entonces la edad del fósil sería aproximadamente de 
35206 años. 


Problemas de mezclas 


El siguiente problema ha sido tomado de Ross [RS, pag. 109]. Se permite que una 
cierta sustancia S fluya con una cierta rapidez en un recipiente que contiene una 
mezcla, manteniéndose ésta uniforme por medio de un dispositivo que la agite. 
Además, esta mezcla uniforme fluye simultáneamente hacia el exterior del recipiente 
con otra rapidez (generalmente diferente). Determinar un modelo que nos permita 
calcular la cantidad de sustancia S presente en la mezcla en el instante t. 

Las condiciones del problema ya nos indican cuales son las suposiciones que 

vamos a tomar en cuenta en la construcción de nuestro modelo. Si x representa la 

dx 

cantidad de S presente en el instante t, entonces la derivada — mide la rapidez de 

dt 

variación de x con respecto a t. Si denotamos “ent” la razón a la que S entra en 
la mezcla y por “sal” la razón a la que sale, entonces el modelo que expresa esta 
situación es 

dx 

— = ent — sal. 
dt 

Para tener una idea más clara de la aplicación de estos modelos matemáticos 
presentamos el siguiente ejemplo: 
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Ejemplo 1.3. La corriente sanguínea lleva un medicamento hacia el interior de un 
órgano a razón de 3 cm 3 /seg, y sale de él con el mismo caudal. El órgano tiene un 
volumen líquido de 125 era 3 . Si la concentración del medicamento en la sangre que 
entra en el órgano es de 0.2 gr/cm 3 : 


a. ¿Cuál es la concentración del medicamento en el órgano en el instante t si 
inicialmente no había vestigio alguno del medicamento? 

b. ¿Cuándo la concentración del medicamento en el órgano será de 0.1 gr/cm 3 ? 

Solución. Sea M(t) la cantidad de medicamento en el órgano al segundo t. Puesto 
que el órgano tiene un volumen líquido de 125 cm 3 , tenemos que la concentración de 
medicamento que sale del órgano es igual a -^^gr/cm 3 . 

Utilizando el modelo matemático para la cantidad de una sustancia en un 
recipiente se tiene que (ver figura 1.4 para una representación gráfica de este ejemplo): 

dM 3 M 

~df ~ °' 6 _ 125’ 

M( 0) = 0 . 


Desde la primera ecuación se tiene que 


dM dt 

-= -. Integrando a cada lado 

75 - 3 M 125 


y realizando los cálculos necesarios se tiene que la cantidad de medicamento en el 
órgano al segunto t está dada por la expresión: 


M(t) = 25 ^1 — e pkV 


Observamos que las integrales se han tomado desde 0 hasta M{t ) y desde 0 hasta 
t respectivamente. 

Finalmente, si escribimos Con(í) por la concentración del medicamento al segundo 


t tenemos que Con(í) = 


M(t) 
125 ' 


Es decir: 


Con(t) = 


1 — e 125 


Para responder al literal b del ejemplo tenemos que resolver la ecuación 


i - M - 

1 — e i25 


= 0 . 1 . 
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Realizando los cálculos respectivos tenemos que t « 29 seg. 



Figura 1.4: Representación de la entrada-salida de un medicamento en un órgano. 

Ley de enfriamiento/calentamiento de Newton 

Si un cuerpo, cuya temperatura es T, se encuentra rodeado o en los alrededores de 
otro cuerpo, cuya temperatura es Tm, termina alcanzando una temperatura igual a 
T m sea por perdida de temperatura (cuando Tm <T) o por aumento de ésta (cuando 

T m > T ). 

Si denotamos por Tm a la temperatura del medio, entonces, de acuerdo con 
la ley empírica de Equilibrio Térmico de Newton, la rapidez con la que cambia la 
temperatura de un cuerpo es proporcional a la diferencia entre la temperatura del 
cuerpo y la del medio que lo rodea. Si además suponemos que la temperatura del 
cuerpo es To cuando t — to, entonces el modelo matemático que expresa esta situación 
es: 

(f = *< r - r -> 

|r(ío) = T 0 , 

donde k es una constante de proporcionalidad. 

A continuación presentamos un ejemplo que muestra una posible aplicación de 
este modelo matemático. 

Ejemplo 1.4. Justamente antes del mediodía, el cuerpo de una víctima de homicidio 
es hallada en un cuarto que se conserva a temperatura constante e igual a 21 °C. A 
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mediodía, la temperatura del cuerpo es de 26 °C ya la una de la tarde, es de 24 °C. 
Considere que la temperatura del cuerpo en el momento de la muerte era de 37 °C y 
que el cuerpo se ha enfriado de acuerdo con la Ley de Equilibrio Térmico de Newton. 
¿Cuál fue la hora de la muerte? 


Solución. Fijemos el mediodía como la hora cero para nuestros cálculos. Si í 0 denota 
la hora de fallecimiento, entonces el problema pide hallat el valor de t 0 tal que se 
satisfaga el siguiente problema: 

dT , 

-=*(T- 21 ) 

T(t„) = 37. 


Desde la primera ecuación, con los datos proporcionados por el problema, se 
tiene 

f U dT r f 1 r 

J 26 T — 21 Jq 


De la última igualdad se tiene que (realizar los cálculos) k — ln ( - 

5 


Nuevamente, desde la primera ecuación tenemos 


p37 dT (3 

= ln - 


'26 


T - 21 


fío 


Calculando la integral vemos que ¿o = 
aproximadamente las 10 de la mañana. 


Mf) 

111 (!) 


dt. 


—2. Luego, la hora de muerte es 


Proliferación de enfermedades 


Una enfermedad contagiosa, por ejemplo, la gripe, se propaga a través de una 
comunidad por personas que han estado en contacto con otras personas enfermas. 
Denotemos por x(t) el número de personas que han contraído la enfermedad y por 
y(t) el número de personas que aún no han sido expuestas al contagio. Si se supone 
que la razón con la que se propaga la enfermedad es proporcional al número de 
interacciones entre estos dos grupos de personas, entonces la ecuación que describe 
esta situación es 

dx 

— = kxy, 
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donde k es una constante de proporcionalidad que se la puede llamar “constante de 
propagación de la enfermedad”. Si además suponemos que, en una comunidad de 
N personas se introduce un número de N 0 enfermos, entonces podemos argumentar 
que x, y están relacionadas por la ecuación x + y = N + N 0 . Luego y = N + N 0 — x. 
Finalmente, el modelo matemático para describir la propagación de una enfermedad 
es el siguiente: 

\~r~ — kx(N + N 0 — x) 

[a(ío) = N 0 . 

Observación 1.2. El modelo matemático presentado para la propagación de 
enfermedades no toma en cuenta las migraciones. Un modelo que tome en cuenta 
este hecho es mucho más complicado (al menos es lo que se espera) que el presentado 
en este texto. 

El siguiente ejemplo es una modificación del problema 7 de Zill |ZD, pag. 28]. 
Este muestra una aplicación elemental del modelo de propagación de enfermedades 
expuesto más arriba: 

Ejemplo 1.5. Suponga que un alumno es portador del virus de la gripe y regresa al 
apartado campus de su universidad de 1000 estudiantes. Si al tercer día el número de 
enfermos es de cinco estudiantes, halle una expresión para el número de estudiantes 
x(t) que contraerán la gripe si la razón con la que la enfermedad se propaga es 
proporcional al número de interacciones entre el número de estudiantes que tienen 
gripe y el número de estudiantes que aún no se han expuesto a ella. 

Solución. Para este problema, consideramos que N = 999, el número de estudiantes 
en el campus antes de reintegrarse el estudiante enfermo, N 0 — 1. Podemos suponer 
que ¿o = 0- Luego el modelo matemático que expresa esta situación es: 

\ ( ~jf = ^(1000 — x) 

[x(0) = 1. 

Desde la primera ecuación y desde la información proporcionada por el problema 
tenemos 

dx 

a:(1000 — x) 
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realizando los cáclulos (hacerlos) se tiene que 


3000 


— ln 



Nuevamente, desde la primera ecuación (realizar los cálculos), tenemos que la 
expresión para la número de estudiantes que contraerán la gripe está dada por: 



Modelos matemáticos generados por problemas diversos 

A continuación se presentan algunos problemas que ilustran la construcción de un 
modelo matemático a partir de situaciones físicas que se presentan en determinados 
problemas de carácter físico y/o técnico. 

Problema 1.1. Hallar un modelo matemático que describa la variación del nivel de 
agua en el tiempo, en un tanque cilindrico de almacenamiento de agua que tiene un 
diámetro D y una altura inicial de agua /i 0 , a partir del momento en que se abre la 
tubería de desfogue inferior con un diámetro d. 

Solución. La figura 1.5 muestra gráficamente las condiciones del problema 1.1. 

D 


h 


d 


Figura 1.5: Representación de un tanque cilindrico de almacenamiento de agua 
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Para hallar el modelo matemático tomamos en cuenta lo siguiente: 


1. Por un lado, el caudal de salida de agua es igual a la velocidad de salida de 

agua multiplicada por el área de la sección transversal. Sea v la velocidad de 

salida de agua. Por el Teorema de Torricelli se tiene que v = y/2gh. Además, el 

d 2 n 

área de la sección transversal es igual a ——. Denotemos con q s el caudal de 
salida del agua. Así, tenemos que q s = —- a/2 gh. 

dV 

Por otro lado, se define q s = —. donde V es el volumen de agua en el cilindro. 

dt d d 2 

Igualando las dos expresiones para q s tenemos = —j— \J 2 gh. Luego 


dV = ——\f2ghdt. 


( 1 . 1 ) 


2. El diferencial de volumen de agua en el tanque cilindrico es igual al área de la 
sección transversal multiplicada por un diferencial de altura. Sea h la altura de 
agua en el tanque, luego se tiene que 

D 2 7T 

dV =-—fifi, (1.2) 

donde el signo ” indica que el volumen va disminuyendo. 


¿ 2 ^ _ Y )2 tt 

Igualando las ecuaciones 1.1 y 1.2 se tiene \j2ghdt =-— dh. Eliminando 

términos semejantes, la última ecuación se puede escribir como 


dh 

dt 


D' 2 


V^gh- 


Si pedimos que h(to) = fio, entonces el modelo matemático que describe la variación 
del nivel de agua en el tiempo, en un tanque cilindrico de almacenamiento de agua 
que tiene un diámetro D y una altura inicial de agua fio, a partir del momento en 
que se abre la tubería de desfogue inferior con un diámetro d está dado por 

(dh d 2 

J dt = ~D 2 

[fi/o) = fio- 
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Problema 1.2. Se denomina curva de persecución a la curva que describe un objeto 
que se desplaza a velocidad constante w, y que persigue de manera óptima a otro que 
se desplaza en línea recta a una velocidad v también constante. 

Hallar un modelo matemático que describa una curva de persecución. 


Solución. Sea A el objeto que describe la curva de persecución (la curva en rojo de 
la figura 1.6) y sea B el objeto que es perseguido por el objeto A. Además, suponemos 
que la distancia entre A y B es c. 



Figura 1.6: Representación de una curva de persecución 


En el instante en el que el objeto B se encuentra en la posición Q, el objeto A se 
encuentra en la posición P. Luego, por la definición de curva de persecución, tenemos 

dy vt — y 
dx c — x 

Desde la física elemental se conoce que la rapidez se define como la variación del 
arco de la curva descrita por un objeto respecto al tiempo. 


(1.3) 
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Figura 1.7: Aproximación de la rapidez de un objeto 


Desde la figura 1.7 se puede ver que 

As « \J (Air) 2 + (Ay) 2 



Dividiendo cada termino por A t se tiene que 


As 

Ai 


límite cuando At —» 0 se tiene 




2 Ax 

A t 


. Tomando 


ds 

dt 


V^W-% 


(1.4) 


Observación 1.3. Dejamos como ejercicio para el lector la justificación del hecho 

que —y — cuando At —y 0. 

Ax dx 


dy 

En la Ecuación 1.3 ponemos p = — y despejamos la variable t. Realizando 

dx 


cuentas se tiene 


t = 


y + p{c — x) 
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Ahora derivamos t respecto a x: 


dt 

dx 


1 ( d y , ( , x , / x dp 


, / x d P 
p - p + (c - x) — 


1 

V 

c — x dp 
v dx 


En la ecuación 1.4 ponemos b = — y despejamos el valor de —. Realizando cuentas 

LLL \JjJU 

dt 

(hacerlo) se tiene — = -. Igualando las dos últimas expresiones, tenemos 

dx b 


c — x dp yA — p 2 

v dx b 


Finalmente, suponiendo que p(0) = 0, el modelo matemático para describir una curva 
de persecución esta dado por: 

{ dp v \J\— p 1 
dx b c — x ’ 

P( 0) = o. 

Problema 1.3. Considere dos tanques A y B conectados entre sí. Suponga que el 
tanque A tiene un volumen Va de agua en el que se han disuelto La libras de sal. 
Suponga que el tanque B tiene un volumen Vb de agua pura. Un mecanismo bombea 
agua pura hacia dentro del tanque A y salmuera hacia afuera del tanque B, existe un 
intercambio de líquido entre los tanques. Construir un modelo matemático que describa 
la cantidad de libras x\(t) y X 2 (t) de sal en los tanques A y B, respectivamente, en 
el tiempo t. 


Solución. La figura 1.8 muestra las condiciones impuestas al problema 1.3. 


Con un análisis similar al que se realizo en el ejemplo 1.3 (ver p. 15) vemos que 
la razón de cambio de sal en el tanque A esta dado por: 


dxi 

dt 


razón de entrada de la sal 

C A -0 + C BA -^ - 

Vb 

Cab Cba 

-Xi + ——£ 2 - 


razón de salida de la sal 

,- A -s 

Xi 


C ™-T a 


V a 


V b 


(1.5) 
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De manera similar, para el tanque B , la razón de cambio de sal es 


dx 2 

dt 


n X 1 n X 2 n x 2 
C AB ---C BA ---C B -~ 


C_BA 

~Va 


■Xi - 


Cba + C, 


B 


Vi 


B 


X 2 . 


( 1 . 6 ) 


agua pura caudal Cba 



Figura 1.8: Tanques de mezclado conectados 


Tomando en cuenta que aq(0) = La y £ 2 ( 0 ) = 0, el modelo matemático que 
describe la cantidad de libras aq(í) y x 2 (t) de sal en los tanques Ay B está dado 
por el siguiente conjunto de ecuaciones: 


r dx\ 
dt 

dx 2 
dt 


Cab Cba 

X\ + —rz — x 2 


V A 

Cba 


V a 

£i(0) = L a 


Xi - 


V B 
Cba + Cb 
Vb 


%2 


x 2 (t) = 0. 


Problema 1.4. El aprendizaje es un proceso extremadamente complejo. Sin 
embargo, es posible construir modelos matemáticos de ciertos tipos de memorización. 
Por ejemplo, considere una persona a quien se le da una lista para estudiar y 
posteriormente se le hacen pruebas periódicas para determinar exactamente qué tanto 
de la lista ha memorizado. Sea L{t) la fracción aprendida de la lista en el tiempo t, 
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donde L = 0 corresponde a no saber nada del listado y L = 1 corresponde a saber 
todo el listado. Evidencia experimental muestra que la variación de L es proporcional 
a la fracción que queda por aprender. 


1. Hallar un modelo matemático para describir L(t). 


2. ¿Para qué valor de L ocurre más rápidamente el aprendizaje? 


Solución. Puesto que la fracción aprendida de la lista en el tiempo t es proporcional 
a la fracciona! que queda por aprender y L = 1 corresponde a saber todo el listado, 
se tiene que 


dL 

dt 


k(l-L), 


donde k es una constante de proporcionalidad. Si además suponemos que L(0) = 0, 
entonces el modelo matemático qne describe la fracción aprendida de una lista está 
dado por 


Tt= k ' l ~ L) 
L( 0) = 0. 


Para responde al segundo literal partimos de la ecuación — = k{ 1 — L ). 


a partir de esta igualdad se tiene que 


dL i.n , dL 

= kdt, luego 


1 — L 


1 — L 


Realizando los cálenlos (hacerlo) se tiene que 


En efecto, 



L(t) = 1 - e~ kt . 


(1.7) 


Un análisis cuantitativo de la ecuación 1.7 nos indica que el aprendizaje ocurre más 
rápidamente cuando L = 0. La figura 1.9 muestra la curva L{t) para algunos valores 
de la constante k. 
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Figura 1.9: Gráfico de la función L(t) para algunos valores de k 


Observación 1.4. El lector puede trazar algunas curvas L(t) para valores negativos 
de k y darse cuenta que L{t) tiene sentido únicamente para valores positivos de k. 

Problema 1.5. Hallar una curva que pase por el punto (0, —6), de tal forma que 
la pendiente de la tangente en cualquiera de sus puntos sea igual a la ordenada del 
punto más siete unidades. 


Solución. Sea y(x) la curva. Puesto que la pendiente de la tangente en cualquier 
dy 

punto es igual a —, el modelo matemático que hallará la curva solicitada es 
dx 


dy _ , 7 

~r~ — y + 7 

dx 

2 /( 0 )) - 6 . 


Para hallar explícitamente y(x) resolvamos el problema matemático planteado por 

dy 

el modelo matemático. Es decir, hallamos y(x) a partir de la ecuación — = y + 7. 

dx 

dy 

Se tiene que - = dx. Luego, integrando a cada lado, ln (y + 7) = x + k, donde 

y+ 7 

k es una constante de integración. Para hallar el valor de la constante k utilizamos 
la condición y = —6 cuando x — 0. Haciendo cuentas, se tiene que k — 0, luego 
ln(y + 7) = x. Despejando la variable y de la última expresión tenemos que la curva 
está dada por la ecuación y(x) = e x — 7. 
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Para más ejemplos de modelos y los problemas matemáticos que generan estos, 
se recomienda consultar la bibliografía. 

Ej ercicios 

El siguiente bloque de ejercicios ha sido tomado de Becerril y Elizarraz [BE], 

1. Encontrar las trayectorias ortogonales de la familia de parábolas y = ex 2 . 

2. Determine el miembro de la familia de trayectorias ortogonales de 3 xy 2 = 2+3 ex, 
que pasa por el punto (0,4). 

3. En 1950 se hicieron excavaciones en Nipur (Babilonia), en las cuales se 
encontraron muestras de carbón que reportaron 4.09 desintegraciones por 
minuto y por gramo. Una muestra actual reportó 6.68 desintegraciones por 
minuto y por gramo. Se sabe que la primer muestra se formó en la época 
del reinado de Hammurabi. Con estos datos, determine hace cuanto tiempo 
Hammurabi reinó en Babilonia. 

4. Un cuerpo a una temperatura desconocida se pone en un refrigerador a una 
temperatura constante de 1 °F. Si después de 20 minutos la temperatura del 
cuerpo es de 40 °F y después de 40 minutos la temperatura del cuerpo es de 
20 °F, hallar la temperatura inicial de este. 

5. En un cultivo de bacterias, se tenían x número de familias. Después de una 
hora, se observaron en el cultivo 1000 familias de la bacteria y después de 
cuatro horas, cuatro familias. Encontrar la expresión para el número de familias 
de la bacteria presentes en el cultivo al tiempo t y el número de familias de la 
bacteria que había originalmente en el cultivo. 

6. Un tanque contiene inicialmente 60 galones de agua pura. Entra al tanque, a 
una tasa de 2 gal/min, salmuera que contiene 1 libra de sal por galón, y la 
solución (perfectamente mezclada) sale de él a razón de 3 gal/min. Obtenga el 
número de libras A(t) de sal que hay en el tanque en un instante cualquiera. 
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¿Cuánto demorará el tanque en vaciarse? ¿Cuál es la máxima cantidad de sal 
que llega a tener el tanque? 


1.2 Ecuaciones diferenciales ordinarias de primer 

orden separables 


En esta sección estudiamos el primer y quizá más sencillo tipo de ecuación diferencial. 
Notamos además que el método usado para hallar la solución de estas ecuaciones 
también fundamenta la forma en que resolvimos algunos de los problemas matemáticos 
planteados por los modelos de la sección 1.1. 

Comenzamos definiendo qué se entiende por ecuación diferencial y dando algunos 
conceptos relacionados con las ecuaciones diferenciales. 

Definición 1.2. Se llama ecuación diferencial a cualquier ecuación en la cual 
aparezcan las derivadas de una o más variables independientes, las incógnitas de la 
ecuación, respecto a una o más variables dependientes 2 . 


Son ejemplos de ecuaciones diferenciales las siguientes: 
i d y a c 

1. — — 1- oxy — sena;. 
dx 

d 2 y dy 

2. —r + 5— + 7x = xy. 

dx dx 

d 2 u d 2 u 
dx 2 dy 2 

. du 

4. — = ux. 
dx 


Sobre el conjunto de todas las ecuaciones diferenciales se puede hacer una primera 
clasificación de acuerdo a la naturaleza de las derivadas que aparecen en la ecuación. 


Definición 1.3. Una ecuación diferencial se llama ecuación diferencial ordinaria 
(EDO) si en ella aparecen únicamente derivadas ordinarias. Si en la ecuación aparecen 

2 En esta definición de ecuación diferencial están excluidas ecuaciones de la forma 

d(uv) du dv , _.. _ , ...... 

—-— = v— —I- u—, (e ) = e , etc. que son identidades. 
dx dx dx 
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derivadas parciales la ecuación se dice que es una ecuación diferencial en derivadas 
parciales (EDP). 

En los ejemplos anteriores, las ecuaciones 1 y 2 son EDO mientras las ecuaciones 
3 y 4 son EDP. 

Una segunda clasificación sobre el conjunto de ecuaciones diferenciales está dada 
por el orden de la derivada de más alto orden que aparece en la ecuación. 

Definición 1.4. Se llama orden de una ecuación diferencial al orden más alto de la 
derivada que aparece en la ecuación. 

Las ecuaciones de los ejemplos 1 y 4 son ecuaciones de primer orden mientras 
que las ecuaciones 2 y 3 son de segundo orden. 

Una ecuación diferencial ordinaria de n-ésimo orden se puede escribir de la forma: 

F(x,y,y',y",...,yW)=0. (1.8) 


La noción de solución de una ecuación diferencial es, claramente, la misma que 
se tiene para una solución de cualquier ecuación. La definición exacta de solución de 
una ecuación diferencial ordinaria es la siguiente: 


Definición 1.5. Se llama solución explícita de la ecuación 1.8 en un intervalo I a 
cualquier función f definida en I, que admite hasta la n-ésima derivada para todo x 
( df d 2 f d n f\ 

de I y tal que F ( x, /, —, —-,...,-) es idénticamente nulo para cada x £ /. 

V dx dx ¿ dx n 1 


Por ejemplo, la función / definida para todo real x por f{x) = 2 sen x + 3 eos a: 

d 2 y 

es una solución explícita de la ecuación diferencial —- + y = 0 en todos los reales. 

dx ¿ 

En efecto, para comprobar que / es una solución de la ecuación diferencial tenemos 
que calcular la segunda derivada y reemplazar este valor en la ecuación diferencial. 


df 9 

— = 2 eos x 

dx 


3 sena:. 


dx 2 


= — 2 sena: — 3 eos x. 
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Reemplazando se tiene 


cPf 
dx 2 


+ / 


—2 sen x — 3 eos x + 2 sen x + 3 eos x, 


que es idénticamente igual a cero en todos los reales 


Observación 1.5. En la mayoría de los casos, lo más que se puede esperar al 
aplicar algún método de resolución de una ecuación diferencial es hallar una relación 
g(x,y,c) = 0 que satisfaga la ecuación diferencial. Una relación de este tipo, es decir, 
una relación que satisfaga la ecuación diferencial, se llama solución implícita si, a 
partir de esta, se puede definir una solución explícita en un intervalo I. 


Por ejemplo, la relación x 2 + y 2 — 25 = 0 satisface la ecuación diferencial 

dy 

x + y— = 0 (comprobarlo). Además, esta relación define la función f(x) = \/25 — x 2 
que es una solución explícita en el intervalo I = [—5,5] (comprobarlo). Luego la 
relación x 2 + y 2 — 25 = 0 es una solución implícita. 


Observación 1.6. Una relación g(x,y,c) = 0 que satisfaga una ecuación diferencial 

pero que, a partir de esta, no sea posible hallar una solución explícita de la ecuación 

diferencial en algún intervalo I, se llama solución formal de la ecuación. Por ejemplo, 

dy 

la relación x 2 + y 2 + 25 = 0 es una relación que satisface la ecuación x + y— = 0, 

dx 

pero claramente se ve que esta relación no define ninguna función que sea solución 
explícita de la ecuación. Luego, la relación x 2 + y 2 + 25 = 0 es una solución formal. 


De ahora en adelante, a menos que se indique lo contrario, solo diremos solución 
de la ecuación diferencial sin indicar si se está tratando con una solución explícita o 
implícita. 

Una ecuación diferencial ordinaria de primer orden es aquella ecuación diferencial 
que tiene o se puede reducir a la forma 

t = nx ’ y) - (L9) 

Proposición 1.1. Toda ecuación diferencial de primer orden se puede escribir de la 
forma M(x, y)dx + N(x, y)dy = 0. 


30 





Leónidas Cerda y Janneth Morocho 


Demostración. Ya que una ecuación de primer orden tiene la forma 1.9, ponemos 
M(x, y) = f(x, y) y N(x, y) = -1. □ 


Una ecuación diferencial de primer orden escrita como M(x, y)dx + N(x, y)dy = 0 
se dice que está escrita en forma diferencial. La forma diferencial de una ecuación de 
primer orden es importante porque a partir de ésta se puede decidir con que tipo de 
ecuación diferencial se está tratando. 


A continuación pedimos hallar la forma diferencial de una ecuación diferencial 
particular. 

x 2 

Ejemplo 1.6. Hallar la forma diferencial de la ecuación y' ---= 0 

sen(x — y) 


dy ¿re¬ 


solución. Escribimos la ecuación como — = - 

dx sen(x — y) 


. Luego, multiplicando y 


trasponiendo términos, tenemos que la forma diferencial de la ecuación es 


x 2 dx — sen(x — y)dy = 0. 


Sea g(x,y,c ) = 0 la solución de la ecuación M(x,y)dx + N(x,y)dy = 0. La 

relación g(x,y,c ) = 0 define una familia uniparamétrica de curvas llamada curvas 

integrales de la ecuación diferencial. Por ejemplo, las curvas integrales de la ecuación 
dy 

diferencial x + y— = 0 es la familia uniparamétrica de curvas x 2 + y 2 — c = 0 donde 
dx 

c> 0. 

En el siguiente gráfico, se pueden apreciar algunas de estas curvas integrales. 



Figura 1.10: Curvas integrales de la ecuación x + y 


dy_ 

dx 


0 
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En muchos problemas es común que aparezca una condición adicional que debe 

satisfacer la solución de la ecuación diferencial. Por ejemplo, se puede pedir hallar 

dy 

la solución de la ecuación — = f(x,y ) de tal forma que la solución, digamos y(x), 
satisfaga la siguiente condición: y = y 0 cuando x = Xq. Este tipo de problemas se 
les conoce como problemas con condiciones iniciales cuando io = Oy problemas con 
condiciones en la frontera en caso contrario. 

La forma habitual de escribir este tipo de problemas es como sigue: 

Resolver el problema: 

J t = 

[í/(®o) = yo- 

El sentido geométrico de un problema con condiciones iniciales (con condiciones 
en la frontera) consiste en elegir del conjunto de curvas integrales justo aquella que 
satisface la condicón adicional. 


El siguiente gráfico muestra algunas curvas integrales y cuál es la curva que está 
eligiendo la condición y{x o) = yo- 



Figura 1.11: Solución al problema 


dy 

dx 


f(x,y),y {x 0 ) = yo 


La noción de problema con condiciones iniciales y problema con condiciones en 
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la frontera se puede extender fácilmente a ecuaciones diferenciales de n-ésimo orden 
de la siguiente manera: 

Resolver: 

F (x,y,y',y",...,y (n) ) = 0 , 

V Oo) = Vo, 

y'(x i) = yi, 

< 

y"(x 2 ) = y 2 , 

í/ (n-1) (®n-l) = Vn-1- 

Cuando Xq — X\ — ... — x n _i, el problema se llama problema con condiciones 
iniciales, en caso contrario, se llama problema con condiciones en la frontera. 

Observación 1.7. Cuando resolvemos un problema con condiciones iniciales (o con 
condiciones en la frontera), es importante de que estemos seguros que tal solución 
existe, sino podríamos pasar toda la vida buscando una solución que podría no existir. 

El siguiente teorema, cuya demostración cae fuera del alcance de este texto, da 
una condición suficiente para la existencia y unicidad de soluciones de una ecuación 
diferencial de primer orden con una condición adicional del tipo y(x o) = yo- 

Sea 1Z = (a, b) x (c, d) un rectángulo del plano xy, (xq, yo) G 1Z. 

Teorema 1.1. Si la función definida por f(x,y ) y su derivada parcial especio a y 
son ambas funciones continuas en 7 Z, entonces el problema 

JÉ = 

[yfao) = y 0 - 

tiene solución única. 

Ahora retomamos nuestro estudio de las ecuaciones diferenciales de primer orden. 
Recordamos que la proposición 1.1 nos dice que toda ecuación diferencial de primer 
orden se puede escribir de la forma M(x,y)dx + N(x,y)dy = 0. 
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Definición 1.6 (Ecuación separable). La ecuación diferencial de primer orden 
M(x,y)dx + N(x,y)dy = 0 se llama ecuación diferencial en variables separables o 
simplemente ecuación separable si M(x,y ) = fi(x)g 2 (y) y N(x,y) = / 2 (x)< 7 i(?/). 

Notamos que, si M(x,y)dx + N(x,y)dy = 0 es una ecuación separable, entonces 
se la puede escribir como 



( 1 . 10 ) 


La solución de una ecuación separable se obtiene integrando, es decir, la solución 
de una ecuación separable es 



Ejemplo 1.7. Resolver la siguiente ecuación diferencial ( y 2 + xy 2 )y' + x 2 — yx 2 = 0. 


Solución. Primero escribimos la ecuación en su forma diferencial. Realizando 
operaciones, vemos que esta ecuación la podemos escribir como 


(. x 2 — yx 2 )dx + ( y 2 + xy 2 )dy = 0, 


por tanto M(x, y) = x 2 — yx 2 = x 2 (l — y) y N(x, y) = y 2 +xy 2 = y 2 ( 1+x). Dividiendo 
toda la ecuación para (1 — y)(l + x) la ecuación queda de la forma 


1 + x 


dx + - 1 -—dy = 0. 


f x 2 

Luego, la solución es / - dx + 

J 1 +x 

de la ecuación diferencial es 


/ 


/ 


——dy = c. Integrando, se tiene que la solución 

i -y 


x 


.2 


y 

x + ln(x + 1) —-— y — ln(l — y) = c. 


2 


Ejemplo 1.8. Resolver el siguiente problema con condiciones iniciales: 



Solución. Primero resolvemos la ecuación xyÓ — y 2 dx+y\/l — x 2 dy = 0. Separando 
variables e integrando se tiene que 


Vi ~ x 2 - Vi - y 


( 1 . 11 ) 
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es la solución de la ecuación. La condición y( 0) = 1 nos permite hallar el valor de la 
constante c. Como la condición inicial nos dice que x — 0, y — 1, reemplazamos estos 
valores en la ecuacuón 1.11 y resulta que c = 1. Luego, la solución del problema con 
condiciones iniciales es \/l — x 2 — \Jl — y 2 = 1 . 


El siguiente ejemplo muestra una ecuación diferencial que no es separable pero, 
con un cambio de variable adecuado, se la transforma en una ecuación diferencial 
separable. 


Ejemplo 1.9. Probar que la ecuación (x + y)dx 
transformar en una ecuación diferencial separable. 



0 se puede 


Solución. Afirmamos que la sustitución x = uy transforma la ecuación diferencial 


x 


(x + y)dx — — ) dy = 0 en un acuación diferencial separable. En efecto, tenemos 

\y J 

que dx = ydu + udy. Reemplazando estos valores en la ecuación diferencial tenemos 


(uy + y)(ydu + udy) 


u 2 y 2 


dy = 0. 


Multiplicando y agrupando términos tenemos (se pide al lector que realice las 
cuentas) y 2 (u + 1 )du + uydy = 0. La última ecuación es claramente una ecuación 
diferencial en variables separables. 

El ejemplo 1.9 muestra que existen ecuaciones diferenciales que son susceptibles de 
transformarlas en ecuaciones separables. La siguiente sección estudia casos particulares 
de estas ecuaciones. 


Ej ercicios 


El siguiente bloque de ejercicios ha sido tomado de Boyce y DiPrima [BD], 


1 . 


Demostar que la ecuación 


dy_ 

dx 


x 2 

i -y 2 ' 


2. Resolver el problema con valor inicial 


dy 3x 2 + Ax + 2 

dx 2 (y — 1) 

2 /( 0 ) = -1 
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3. Encontrar una solución al problema con valor inicial 


dy y eos x 
dx 1 + 2 y 2 

2 /( 0 ) = 1 


Hallar la solución general de las siguientes ecuaciones diferenciales: 

4. y' = (eos 2 x)(eos 2 2 y). 

5. xy' = (1 — y 2 ) 1 / 2 - 

, x~e~ x 

6 . y = -. 

y + yv 

dr r 2 

7 ' dé = l' 

8. sen 2x dx + eos 3 y dy = 0. 

9. y 2 (l — x 2 ) l ¡ 2 dy = are cosa: dx. 

1.3 Ecuaciones diferenciales ordinarias que se 
pueden reducir a ecuaciones separables 


En esta sección, estudiamos algunos tipos de ecuaciones que se pueden reducir a 
ecuaciones separables a través de sustituciones. 

Las ecuaciones diferenciales de primer orden homogéneas (ver definición 1.8) 
son una clase típica de ecuaciones diferenciales que pueden reducirse a ecuaciones 
separables. Antes de dar la definición de ecuación diferencial homogénea tenemos 
que dar la siguiente definición. 

Definición 1.7. Sea f una función de m variables. Se dice que f es una función 
homogénea de grado n si para todo X G M su cumple que 

/( Xxi, Xx 2 , • • •, Xx m ) = X n f(xi, x 2 ,..., x m ). 

En el caso particular de una función de dos variables x, y se tiene que / es una 
función homogénea de grado n si para todo número real A se cumple que 

f(Xx,\y) = X n f(x,y). 
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Son ejemplos de funciones homogéneas las funciones definidas como: 


1 . 

2 . 

3. 

4. 


f{x,y) = x 2 + xy. 

f(x, y) = \jx 2 + y 2 + 3x - 2 y. 

/(x, y) = x 3 eos ( — ) — x 2 ye*. 


f(x,y) 


\y. 

x 2 + 3 xy 
5 y 2 + x 2 


La función del ejemplo 1 es una función homogénea de grado 2, las funciones de los 
ejemplos 2, 3 y 4 son funciones homogéneas de grado 1, 3 y 0 respectivamente. Un 
ejemplo de función que no es homogénea es la función definida por /(x, y) = x 2 +2 x—y. 

Ahora ya podemos definir el concepto de ecuación diferencial homogénea. 


Definición 1.8. La ecuación diferencial M(x,y)dx+N(x,y)dy = 0 se llama ecuación 
diferencial homogénea si las funciones M y N son ambas funciones homogéneas del 
mismo grado. 


La ecuación diferencial (x — y)dx + xdy = 0 es una ecuación homogénea, pues 


M (Ax, Xy) = Ax — A y 
= A(x - y) 
= A M(x,y) 


y 


N(Xx, Xy) = Ax 

= A N(x,y). 

Es decir, M y N son ambas funciones homogéneas de grado 1. 


El siguiente teorema nos dice que toda ecuación diferencial homogénea se 
puede transformar en separable. Su demostración nos indica cómo realizar esta 
transformación. 


Teorema 1.2. Si M(x,y)dx + N(x,y)dy = 0 es una ecuación diferencial homogénea, 
entonces un cambio de variable la transforma en separable. 
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Demostración. Sea x = uy. Derivando implícitamente tenemos dx = ydu + udy , 
reemplazamos estos valores en la ecuación diferencial. Ya que las funciones M y N 
son ambas funciones homogéneas del mismo grado tenemos 


y n M(u, 1 )(ydu + udy ) + y n N(u, 1 )dy = 0. 


Luego yM(u, 1 )du + ( uM(u , 1) + N(u, 1 ))dy = 0. La última ecuación es claramente 
una ecuación diferencial en variables separables. □ 


Observación 1.8. Utilizando la sustitución y — vx también se consigue una ecuación 
diferencial en variables separables (realizar los cálculos). 

fx 2 \ 

Ejemplo 1.10. Resolver la ecuación (x + y)dx — ( — ) dy — 0. 

V V ) 


Solución. En este caso, se tiene que M(x, y) = x + y, N(x, y) 
funciones M y N son ambas funciones homogéneas: 


-. Veamos si las 

y 


M(Xx, Xy) = Xx + Xy 
= X (x + y) 
= X M{x,y) 

Para la función N tenemos 

Ar/ . . (Xx) 2 

N(Xx, Xy) = --— 


X 2 x 2 



= X N(x.y). 


Puesto que las funciones M y N son ambas funciones homogéneas de grado 1, la 
ecuación diferencial es homogénea. Por el teorema 1.2 la sustitución x = uy transforma 

f x 2 \ 

la ecuación (x + y)dx — ( — ) dy = 0 en una ecuación diferencial separable. En efecto, 

\V J 

la ecuación se transforma en y 2 (u + 1 )du + uydy = 0 (ver ejemplo 1.9 de la pag. 35). 
Resolviendo esta ecuación separable (realizar las cuentas) se tiene que ln (uy) + u = c 

X 

es su solución. Realizando la sustitución inversa se tiene que la relación lnx H— = c 

y 

es la solución de la ecuación diferencial original. 
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Ejemplo 1.11. Resolver la ecuación xy' = a Jy 1 — x 2 . 


Solución. Escribimos la ecuación en forma diferencial: 

\/ y 2 — x 2 dx — xdy = 0 

Como 

M(Xx, Xy) = s]\ 2 y 2 - \ 2 x 2 

= A \Jy 2 — x 2 

= X M(x,y). 

y 

N(Xx, Xy) = —(Xx) 

= A N(x,y), 

la ecuación es homogénea. Ponemos x = uy. Luego la ecuación se transforma en 


y 


Vi — u 2 du + yu(yj 1 — u 2 — 1 )dy = 0. 


Esta ecuación diferencial es una ecuación en variables separables. El lector puede 
comprobar que la solución de esta ecuación es 

^ + ln y — c. 

X 

Para completar la resolución del ejemplo reemplazamos u — — en la última ecuación. 

y 

Observación 1.9. La ecuación (a\X + b\ r y + ci)dx + (a 2 x + b 2 y + c 2 )dy = 0 no es 
una ecuación homogénea cuando por lo menos uno de los coeficientes ci o c 2 son 
diferentes de cero. 



Este tipo de ecuaciones se pueden transformar en ecuaciones homogéneas cuando 
el sistema 


{ arx + biy + ci = 0, 
a 2 x + b 2 y + c 2 = 0 


( 1 . 12 ) 


tiene sulución única. 


Teorema 1.3. Si el sistema 1.12 tiene solución única, entonces la ecuación diferencial 
(a\x + b\y + cfidx + (a 2 x + b 2 y + c 2 )dy = 0 se puede transformar en una ecuación 
diferencial homogénea. 
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Demostración. Sea (xo ; 2/o) I a solución del sistema 1.12. Ponemos 

y = u + yo- 

Tenemos dx = dy y dy = du. Sustituyendo en la ecuación diferencial resulta que 

(ai y + a±x 0 + biu + biyo + c±)dy + (a 2 y + a 2 x 0 + b 2 u + b 2 y 0 + c 2 )duj = 0. (1.13) 

Ya que (x 0 ,y 0 ) es la solución del sistema 1.12 tenemos que a x x 0 + biy Q + C\ = 0 y 
a 2 xo + b 2 yo + c 2 = 0. Luego, la ecuación 1.13 nos queda como 

(aifj, + b\u)dy + (« 2 /i + b 2 (jj)d(jj = 0, 


que es claramente una ecuación diferencial homogénea. 

Ejemplo 1.12. Resolver la ecuación (x + y)dx + (x — y — 1 )dy = 0. 


□ 


Solución. La solución del sistema 



es el punto ^ . Ponemos x = y+ -yy=iü — Reemplazando estos valores 

y simplificando nos queda (/i + u)d¡i + (/i — co)doj = 0 (realizar los cálculos) que es 
una ecuación diferencial homogénea. Ahora utilizamos el cambio de variable que se 
indica en la demostración del teorema 1.2. Es decir, ponemos y = v lo, derivando 
implícitamente la última igualdad obtenemos dy — u dv + v du. Reemplazamos 
estos valores en la ecuación diferencial (y + cu)dy + (y — co)du = 0; simplificando y 
agrupando términos obtenemos la ecuación u 2 (v + 1 )dv + c o[y 2 + 2v — l)du = 0 que 
es una ecuación separable. La última ecuación la podemos escribir como 


v + 1 1 

-dv H —du = (J. 


v 2 + 2v — 1 


u 


(1,14) 


v + i . Í1 

—- dv + / — du = c. 

v 2 + 2v — 1 I u 


La solución de la ecuación 1.14 está dada por J 
Es decir, la solución de la ecuación 1.14 es - ln(u 2 + 2v — 1) + lncu = c; luego, 


u 


yjv 2 T 2v — 1 = c. 
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IJ- 


Puesto que v — —, la solución de la ecuación (¡a + u:)dfi + (¡i — uj)du¡ = 0 esta 


ÜJ 


dada por la relación ¡i 2 + 2/rcu — oj 2 = c (hacer las cuentas). Finalmente, ya que 
H — x — - jcj — y+-, la solución de la ecuación diferencial 


(x + y)dx + (x — y — l)dy = 0 (la ecuación de partida) 


está dada por la relación [x — - j + 2 (x — 2 ) \ + 2 


y + l 


Ejemplo 1.13. Resolver la ecuación ( 3x + 4 y — 11 )dx — ( 2x + 5 y — 12 )dy = 0. 


Solución. La solución del sistema 

{ 3a; + Ay — 11 = 0 
-2x - 5y + 12 = 0 

es el punto (1,2). Ponemos i = /i + ly|/ = w + 2. Reemplazando estos valores en la 
ecuación diferencial nos queda la euación homogénea (3yU + 4 u)d¡i — (2/i + 5 u>)du> = 0. 
Para resolver esta ecuación utilizamos el mismo procedimiento empleado en los 
ejemplos 1.11 y 1.10 para hallar la solución de una ecuación homogénea. La solución 
de la ecuación original se la consigue con un proceso análogo al utilizado en el ejemplo 
1 . 12 . 


Ej ercicios 


El siguiente bloque de ejercicios ha sido tomado de O’Neil [OP]. 


1. Comprobar si la función definida por f(x, y) 


y y 

— H— es homogénea. 
x z x 


3 d V 


2. Comprobar si la ecuación diferencial x á — = x z y — 2y á es homogénea. 

dx 


Resolver las siguientes ecuaciones diferenciales: 


dy_ 

dx 




2 


4. 


dy_ 

dx 


Í 2x + y-l \ 2 
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En cada uno de los siguientes ejercicios determine si la ecuación es homogénea o 
reducible a homogénea. Si lo es, hallar su solución general. Si no lo es, no intente 
resolverla en este momento. 


c d v 2 

5. x— = y . 

dx 

/ . dy 

6. (* + y) s = y- 

7. (3x — y — 9 )dx — (x + y + 1 )dy = 0. 


. x — y + 6 
V = 


3x — 3y + 4 


n 3 dy 2 3 

9. x — = x y — y . 
dx 


10. = y + AJxy. 

dx 


1.4 Ecuaciones diferenciales exactas y factor de 

integración 


Antes de definir una ecuación diferencial exacta recordamos que si F es una función 
de dos variable que tiene derivadas parciales de primer orden continuas, el diferencial 
total de F está definido como 


dF 


dF 

dx 


dx + 


dF 

dy 


dy. 


Definición 1.9. La ecuación diferencial M(x, y)dx + N(x, y)dy = 0 se dice que 
es una ecuación diferencial exacta si existe una función de dos variables F tal que 
dF = M(x, y)dx + N(x, y)dy. 


El siguiente teorema nos dice cómo saber si una ecuación diferencial es exacta. 

Teorema 1.4. Si las funciones M y N tienen derivadas de primer orden continuas, 

entonces la ecuación diferencial M(x,y)dx + N(x,y)dy = 0 es exacta si y solo si 
dM _ dN 
dy dx 
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Demostración. Supongamos que M(x,y)dx + N(x,y)dy = 0 es una ecuación 
diferencial exacta. Por definición se tiene que existe una función de dos variable F 
tal que dF = M(x, y)dx + N(x, y)dy , luego 


dF dF Ar 

~ñ~ = M y = N - 
dx ay 

Ya que las funciones M y N tienen derivada de primer orden continua se tiene, por 
el Teorema de Schwarz, 

d 2 F d 2 F 


dydx dxdy ’ 


dM dN 
luego —— = —. 
dy dx 


dM dN 

Ahora supongamos que —— = ——. Tenemos que hallar una función de dos 

dy dx 


variables F tal que 


dF 

Hx ~ 
ay 


(1.15) 

(1.16) 


Afirmamos que la función F definida por: 


F(x,y ) = / M(x,y)dx + g(y) con g(y) = \ ( N(x,y) - ^ J M(x,y)dx ) dy. 


d 


cumple las condiciones que se piden en las ecuaciones 1.15 y 1.16. 


Antes de demostrar que esta afirmación es correcta, comprobemos que la función 
g es efectivamente una función que depende únicamente de la variable y. Es suficiente 
ver que 


(N(x, y) — [ M(x,y)dx ) = 0 


dx 


dy 


Resulta que 


d 


dx 


d 


— N(x,y ) - tt / M(x,y)dx = 


dy 


dN(x, y) 

- 32 [a 

dx 

dxdy J 

dN{x , y) 

- 32 [t 

dx 

dydx J 

dN(x, y) 

d ( d 

dx 

dy \dx J 

dN{x , y) 

dM(x , y) 

dx 

dy 


= 0 . 


43 














INTRODUCCIÓN A LAS ECUACIONES DIFERENCIALES ORDINARIAS 


Ahora veamos que la función F en verdad satisface las ecuaciones 1.15 y 1.16. 
La ecuación 1.15 es inmediata. Por otro lado tenemos que 


of d r . d 

7r = 7T M(x,y)dx + — 
dy dy J dy 


N(x, y)dy — / ( ^ j M(x,y)dx)dy 


di d í 

= — y M(x,y)dx + N(x,y) « — J M(x,y)dx 

= N(x,y). 


□ 


Observación 1.10. En la demostración del teorema l.f también podemos utilizar 
la función F definida como 

F(x,y) — j N(x,y)dy + h{x) con h(x) — j ^M(x,y) — — j N(x,y)dy^j dx. 

Si M(x,y)dx + N(x,y)dy = 0 es una ecuación diferencial exacta, entonces su 
solución está dada por F(x, y) = c. La verificación de que esta es efectivamente la 
solución de la ecuación diferencial exacta se deja como ejercicio para el lector. 

El siguiente ejemplo muestra el procedimiento estándar para hallar la solución 
de una ecuación diferencial exacta. 

Ejemplo 1.14. Resolver la siguiente ecuación diferencial: 

( x 3 + xy 2 )dx + ( x 2 y + y 3 )dy = 0. 


Solución. Primero comprobamos que la ecuación es exacta: 

dM d(x 3 + xy 2 ) ^ d(x 2 y + y 3 ) dN 
dy dy dx dx 

Ya que la ecuación es exacta, existe una función de dos variables F tal que 

(dF „ r , 2 

—— = M = x 3 + xy~, 
dx 

dF AT 2 3 

= N = x 2 y + y 3 . 

{ dy 

Integrando respecto ax se tiene que 

F{x, y) = I(x 3 + xy 2 )dx + g(y) 

4 2 2 

x x y , \ 

= ^ + ^- + g(y) 
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Para hallar la función g derivamos F con respecto a ye igualamos este resultado a 

N. 

dF 2 , dg 2,3 

— = xy + — = x y + y . 
óy dy 

dg y ^ 

Esta última igualdad se cumple solamente cuando — = y 3 . Luego g(y) = —, 

dy 4 

reemplazando en la expresión para F(x, y) se tiene que la solución de la ecuación 
diferencial está dada por 

4 2 2 4 

x x y y 

— H---h — = c. 

4 2 4 

Ejemplo 1.15. Hallar la solución de la ecuación diferencial 


( e 2y — y cos(xy))dx + (2xe ¿y — xcos(xy) + 2 y)dy = 0. 

Solución. Ya que = 2e 2y — eos (xy) + xysen(xy) = — la ecuación es exacta. 

dy ox 

Luego existe una función F en las variable x, y tal que 

' dF 

—— = M = e 2y — y eos (xy) 
dx 

dF 

= N = 2xe 2y — x cos(xy) + 2y. 

I dy \ y 

Integrando respecto ax se tiene que 

F(x, y) = I(e 2y - y eos {xy))dx + g(y) 

= xe 2y — sen (xy) + g(y)- 


Ahora derivamos respecto a y la expresión hallada para F(x,y). Se tiene que 


dF 

dy 


2xe 2y — íccos (xy) + g'(y)- 


Como — — = N, tenemos 2xe 2y — íccos (xy) + g'(y ) = 2xe 2y — xcos(xy) + 2y, luego 
dy 

g'(y ) = 2y, por tanto g(y) = y 2 . Así, la solución de la ecuación diferencial de este 
ejemplo está dada por 

xe 2y — sen (xy) + y 2 = c. 


Observación 1.11. Consideremos la ecuación diferencial 

2 xy ln ydx + (íc 2 + y 2 \Jy 2 + 1 )dy = 0. 
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En este caso se tiene que M[x,y ) = 2xy ln y mientras N(x,y) = x 2 + y 2 \J y 2 + 1, de 
dM dN 

esta forma —— ^ ——. Luego, la ecuación diferencial no es exacta, 
ay dx 


Multiplicando la ecuación diferencial por el factor y = - tenemos 

y 


M(x,y ) = 


2 xy ln y 

y 


N(x, y) = 


x 2 + y 2 \/y 2 + 1 


y 


rr , , dM dN T 

Haciendo las cuentas respectivas vemos que —— = ——. Luego, la ecuación 

ay ox 

diferencial es exacta. 


El factor y ( 1/y en la observación 1.11) se llama factor de integración. La 
definición de factor de integración de una ecuación diferencial es la siguiente: 

Definición 1.10. Una función y(x, y) se llama factor de integración para la ecuación 
diferencial M(x,y)dx + N(x,y)dy = 0 si la ecuación 

y(x, y)M(x , y)dx + y(x, y)N(x , y)dy = 0 

es una ecuación diferencial exacta. 


El procedimiento para hallar el factor de integración, siempre que sea posible, es 
el siguiente: 

A partir de la exactitud de la ecuación y(x, y)M(x, y)dx + y(x, y)N(x, y)dy = 0 
tenemos 

d d 

—y(x,y)M(x,y) = —y(x,y)N(x,y). 

Derivando parcialmente tenemos 

9 í \ i\/r/ x w/ ,dy{x,y) dM(x,y) 

y) M y) = y) —^— + Mu y) —^— 


y 


Luego 




m(x. y) d ->M^ + Kx ,y) = N{XtV) Mm> +Kx¡v) ^ y ) 


dy 


dy 


dx 


dx 
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por tanto 
H(x,y) 


dM(x,y) dN(x,y) 


dy 


dx 


= N(x,y) 


dy{x,y) 

dx 


M(x,y ) 


d»{x,y) 

dy 


(1.17) 


La ecuación 1.17 es una ecuación diferencial en derivadas parciales muy difícil 
de resolver si no se tiene algún conocimiento adicional sobre fj,(x,y)\ sin embargo, 
suponer que y(x, y) depende de una sola variable la vuelve muy fácil de resolver. 


Supongamos que y(x,y) = f(x). En este caso la 1.17 se reduce a 

Arreglando términos se tiene 

df(x) = 1 / dM(x,y ) _ dN(x,y) \ 

f(x) N(x,y ) V d y dx ) X ’ 

integrando y despejando f{x) tenemos 

f, ^ ( dM n X ’ V) ~ ) dx 

t(x) — c J N ( x ’y)\ °y °x ) 


(1.18) 


Observamos que una condición necesaria y suficiente para que el factor de 

integración dependa solo de la variable x es que la expresión 

1 ídM(x,y) dN(x,y)\ 

N(x,y) \ dy dx ) 

dependa solo de la variable x. 


Ahora supongamos que ¡i(x,y) = g(y). En este caso, la ecuación 1.17 se reduce a 


fdM(x,y ) dN(x,y)\ 
g(y) {—dy - dx J 


-M(x,y) 


dg(y) 

dy 


Arreglando términos, se tiene 


dg(y) _ 1 / dM(x,y) _ dN(x,y) \ 

g{y) M(x, y) V dy dx ) V ’ 

integrando y despejando g(y) tenemos 


g(y) = e f 


1 / dM(x,y) 

M(x,y ) V dy 


dN(x,y) 

dx 


)dy 


(1.19) 


Observamos que una condición necesaria y suficiente para que el factor de 
integración dependa solo de la variable y es que la expresión 

1 ídM(x,y) dN(x,y)\ 

M(x,y) \ dy dx J 

dependa solo de la variable y. 
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Ejemplo 1.16. Resolver la ecuación (2 xy 2 — 3 y 3 )dx + (7 — 3 xy 2 )dy = 0. 


Solución. Como 


<9(2 xy 2 — 3 y 3 ) 2 0 2 <9(7 — 3 xy 2 ) 


dy 


= 4 xy - 9y ± -3 y = 


dx 


la ecuación 110 es exacta. Se puede comprobar (hacerlo) que el factor de integración 
no depende solo de la variable x, veamos si este depende solo de la variable y. 

U- ({Axy - 9i- + 3yi) 


M (x, y) V dy 


2 xy 2 — 3 y 3 
2y(2x - 3 y) 
y 2 [2x - 3 y) 
2 

y' 


Luego, según la ecuación 1.19 el factor de integración está dado por: 


9{y) = e f2ydy 


= e~ 2lny 


= V 


-2 


Multiplicando la ecuación diferencial por este factor de integración obtenemos la 
(nueva) ecuación 

(2x — 3 y)dx + (7 y~ 2 — 3 x)dy = 0. 

Para esta ecuación vemos que se satisface la condición de exactitud. En efecto, 
tenemos 

d(2x — 3 y) <9(7 y~ 2 — 3x) 

dy dx 

A partir de este punto resolvemos la ecuación diferencial con el procedimiento que 
vimos anteriormente. 

Ejemplo 1.17. Hallar un factor de integración para la ecuación fx 2 + y)dx — xdy = 0. 


Solución. Como 


<9(-x) = _ x , l = d(x 2 + y) 


dx 


dy 
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la ecuación no es excata. Vea mos si el factor de integración depende únicamente de 
la variable x. 

1_ / dM(x,y) _ dN(x,y) \ _ J_ 

x , y) \ dy dx J —x y 

2 
x 

Luego, según la ecuación 1.18, el factor de integración está dado por: 

m = 

_ g— 2\nx 

= x~ 2 . 

El lector puede comprobar sin dificultad que la función f(x) = x~ 2 es un factor de 
integración de la ecuación ( x 2 + y) dx — xdy = 0. 

Observación 1.12. Como se indicó anteriormente, hallar un factor de integración 
es una tarea relativamente sencilla si este factor depende únicamente de una sola 
variable. Sin embargo, cuando el factor de integración depende de las dos variables y 
además tenemos alguna condición adicional sobre este factor, es posible calcularlo. 

Ejemplo 1.18. Resolver la ecuación (3 y 2 — x)dx + (2 y 3 — 6 xy)dy = 0 si se sabe que 
un factor de integración y(x,y) depende de las variables x,y según la ley ip(x + y 2 ), 
es decir ¡i = tp(x + y 2 ). 



Solución. Se deja como ejercicio comprobar que la ecuación no es exacta y 
además comprobar que el factor de integración no depende de una sola variable. 
El procedimiento que se sigue en este caso es el siguiente: Ponemos z = x + y 2 y 


calculamos las derivadas parciales de z respecto a las variables x, y. Tenemos 


{ dz 
dx 
dz 

dy = V ' 


Aplicando la regla de la cadena tenemos 


dji 
dx 
dy 
, dy 


dy dz 
dz dx 
dy dz 
dz dy 


dy 
dz ’ 


2 y 


dy 

dz 
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Notamos que —— = — está justificada por la condición u = (p(x + y 2 ) = <f(z). 
dz dz 

Reemplazando estos valores en la ecuación 1.17 tenemos 

. . (dM(x,y) dN(x,y)\ ■.da(z) „ „ r/ . du(z) 

Mu> -- = N(x,y)^l - 2 yM(x,y)-^ 

\ ay ox J az az 


Lnego 


l*(z. 


-d/a(z) = 


_1_ í dM{x, y) _ dN{x,y) \ ^ 

N(x, y) — 2yM(x, y) \ dy dx 1 


12 y 


-4y(x + y 2 ) 

= ~-dz. 
z 


dz 


Integrando y despejando el factor de integración se tiene que /a(x,y) — (x + y 2 ) 3 . 
Multiplicando la ecuación diferencial por este factor y(x, y) tenemos la nueva ecuación 


3 y 2 - x 2 y 3 - 6 xy 

;dx + —-;- IWZ-dy = 0 . 


( 1 . 20 ) 


(x + y 2 ) 3 (x + y 2 ) 3 

Comprobemos que esta ecuación diferencial es exacta. 

La derivada parcial de M(x, y) con respecto a la variable y es: 

dM(x,y) d / 3 y 2 — x \ 

dy dy \(x + y 2 ) 3 ) 

_ 12 y(x-y 2 ) 

(x + y 2 ) 4 

La derivada parcial de N(x,y ) con respecto a la variable x es: 

dN(x,y) d Í2y 3 — 6xy\ 
dx dx \ (x + y 2 ) 3 J 
_ 12 y(x - y 2 ) 

(x + y 2 ) 4 

Como la ecuación 1.20 es excata, existe nna función F en las variables x,y tal que 


' dF 

3 y 2 — x 

dx 

(x + y 2 ) 3 

dF 

2 y 3 — 6 xy 

. dy 

(x + y 2 ) 3 


Integrando respecto a la variable x. tenemos F(x, y) = 


3 y 2 


x 


(x + y 2 ) 3 


dx. Puesto que 
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x 


1 


y 


(x + y 2 ) 3 (x + y 2 ) 2 (x + y 2 ) 3 

3 y 2 — x 


, tenemos 


F(x,y) = 


(x + y 2 ) 3 


dx 


= 3 y 2 


dx 


x 


(x + y 2 ) 3 


(x + y 2 ) 3 


3 y 


2 1 
+ 


2 {x + y 2 ) 2 x + y 2 2(x + y 2 ) 2 

2 » 2 + -C + %). 


dx + h(y) 


+ %) 


(x + y 2 ) 2 x + y 2 

3 F — Gxy 

Derivando esta expresión respecto a la variable y tenemos —— = —- 77— + h'(y). 

ay (x + y 2 ) 3 

Luego h'{y) = 0, por tanto h(y) = k, donde k es una constante. Así, la solución de la 


ecuación diferencial es 


2 y 2 , 1 

H- 77 = c. 


(x + y 2 ) 2 x + y 2 

Ejemplo 1.19. Hallar un factor de integración para la ecuación diferencial de primer 
orden ( x 2 + y 2 + 1 )dx — 2 xy dy = 0 si se sabe que este depende de las variables x, y 
según la ley <p(x, y) = y 2 — x 2 . 


Solución. Ponemos z = y 2 — x 2 . Aplicando la regla de la cadena se tiene 

dy dy 

dx dz ’ 

dy dy 
~dy = 2V Tz 

„ dM dN , 

Como —— = 2y y —— = —2 y, reemplazando en la cuación 1.17 tenemos 
dy dx 

2 y (y 2 — x 2 + l)-j- = —4 y y, luego (z + 1)-/- = —2 y. Integrando, se tiene que 


dz 


dz 


= {y 2 -x 2 + 1 ) 2 . 

Se deja como ejercicio para el lector la comprobación de que ( y 2 — x 2 + 1)~ 2 es 
efectivamente un factor de integración de la ecuación ( x 2 + y 2 + 1 )dx — 2xy dy = 0. 

Ej ercicios 

Resolver las siguientes ecuaciones diferenciales: 


1 . 2 xy dx + ( x 2 — 1) dy = 0. 
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dy x — y 2 
dx 2 xy + y 



dy — 0 


En cada una de las siguientes ecuaciones hallar un factor de integración y. 

4. (2 xy + y 4 ) dx{3x 2 + 6 xy 3 ) dy = 0, y = <p(y). 

5. (1 + y) dx + (1 — x) dy = 0, y = <p(x). 

6 . (3 y 2 — x) dx + (2 y 3 — 6 xy) dy = 0, y = ip(x + y 2 ). 

Ecuaciones lineales de primer orden 

Un tipo particularmente importante de ecuaciones diferenciales de primer orden son 
las ecuaciones lineales. Su forma es la siguiente: 



( 1 . 21 ) 


donde p(x)q(x) ^ 0. Cuando p{x) = 0 o q(x) = 0, la ecuación 1.21 se reduce a una 


ecuación diferencial separable. 

Para hallar la solución de la ecuación 1.21 la escribimos en su forma diferencial 


(q{x) — p(x)y)dx — dy = 0. 


Se puede comprobar fácilmente que esta ecuación no es exacta. Veamos que esta 
ecuación tiene un factor de integración que depende únicamente de la variable x. En 
efecto, 



Se tiene que un factor de integración de la ecuación es e / pUU*_ Luego, la ecuación 


e -T p ( x)dx { q (x) - p(x)y)dx - e^^ x)dx dy = 0 


es una ecuación diferencial exacta, por lo tanto existe una función de dos variables F 
tal que 


' dF 
dx 


fp(x)dx (g(x) -p(x)y), 


dF 
, dy 


— _ e Spi x ) dx 
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De la segunda ecuación tenemos 

F{x,y) = -ye$ p{ - x)dx + h(x). 

Puesto que la función F tiene que satisfacer la primera ecuación tenemos 

f = _ vp(x)e t + 

dx dx 

= ef p W dx (q(x) -p(x)y). 

dh(x) 

Luego —-— = q(x)e 3 p ( x ) dx . Esta última ecuación nos dice que el valor para h(x) 
dx 

está dado por: 


h(x) = / q{x)ef p ( x)dx dx. 


La solución de la ecuación 1.21 es —y(J pí ' x ' )dx + J q(x)e$ p ^ dx dx = c. Despejando 
y nos queda 

V = eS^(J^^ + o). (, 22 , 

Ejemplo 1.20. Resolver la ecuación (3 x 2 y — x 2 )dx + dy = 0. 

dy 

Solución. Primero escribimos la ecuación en la forma ——f- p(x)y = q (x). Dividiendo 

dx 

la ecuación por dx y trasponiendo términos se tiene 

dy | q 2 2 

-—b ¿x y = x . 
dx 

La ecuación tiene un factor de integración de la forma ¡i{x) = ef 3x2dx = e x3 . Luego, 
la ecuación x 2 e x3 (3y — 1 )dx + e x3 dy = 0 es exacta, por tanto existe una función F en 
las variables x, y tal que 

9F i \ 

d^ =xe 
dF x3 
dy ~ £ ' 

Integrando la segunda ecuación respecto a la variable y tenemos F(x, y) = e x y + h(x). 

dF 3 

Ahora, derivando F respecto a x tenemos = 3x 2 e x y + h'(x). Igualando esta 

dx 

-r-3 


expresión con la primera ecuación se tiene h’(x) = —x 2 e x . Así, h(x) = -. La 

X 3 

3 6 

solución de la ecuación diferencial esta dada por e x y -— = c. Despejando y se 

tiene 

y = e '* 3 í V + c I • 


53 





INTRODUCCIÓN A LAS ECUACIONES DIFERENCIALES ORDINARIAS 


Ejemplo 1.21. Resolver la siguiente ecuación diferencial: 


dy 

— + 2 xy = 2xe 
dx 


Solución. Para utilizar la fórmula 1.22, primero identificamos las funciones p y q. En 
este ejemplo, se tiene que p(x) = 2x y q{x) = 2xe~ x2 . Ahora calculamos e~ ^ p ^dx. 

e -fp^dx = e~f 2xdx 


Luego 


y — e -fp(?)<te 


q(x)eS p( - x)dx dx + c 


= e 


2xe x e x "dx + c 


I 2xdx+C 

= e~ x2 (x 2 + c). 

Observación 1.13. La ecuación 1.21 se dice que es lineal en la variable y. Una 
ecuación lineal, por ejemplo, en la variable x tiene la forma 

C1T 

-+p(y)x = q(y). 



Ej ercicios 

Los siguientes ejercicios han sido tomados de O’Neil [OP]. 
Resolver las siguientes ecuaciones diferenciales: 

1. y' + y = sena:. 

2. y H— y = 3x 2 . 

x 

3. sen(2x)y' + 2 sen 2 (x)y = 2 sen(a:). 

4. y' + xy = cosx. 

5. (x 2 — x — 2)y' + 3a :y = x 2 — Ax + 4. 
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6. (1 - 2xe 2y ) ^ - e 2y = 0. 

Ecuaciones de Bernoulli 

Se llama ecuación diferencial de Bernoulli a una ecuación diferencial de la forma 

^ +p(x)y = q(x)y n , (1.23) 

donde n es un número real diferente de 0 y 1. 

Observación 1.14. Para resolver la ecuación 1.23 primero la multiplicamos por 
y~ n y luego la reducimos a una ecuación diferencial lineal de primer orden con la 
sustitución z = y x ~ n . 


En efecto, multiplicando la Ecuación 1.23 por y n nos queda la nueva ecuación 


V n< Jx +P ^ yl " = Q 


dz 


~n d V 


Ahora usamos la sustitución z = y 1 n . Tenemos que — = (1 — n)y n —. Luego 


n d V 


1 dz 


y dx 1 — ndx 
podemos escribirlo como 


. Sustituyendo estos valores tenemos 


dx 
1 dz 
1 — n dx 


dx 

+ p(x)z = q(x) que 


-h (1 — n)p(x)z = (1 — n)q(x). 

dx 

Esta última ecuación es, claramente, una ecuación diferencial lineal en la variable z. 
Ejemplo 1.22. Resolver la siguiente ecuación diferencial: 

dy . 2 , 

x ——b y = y Inx. 
dx 

Solución. Primero dividimos la ecuación por x para que esta tenga la forma estándar 

de una ecuación de Bernoulli. El resultado es la siguiente ecuación: 

dy y y 2 ln x 

dx x x 

Ahora multiplicando esta última ecuación por y~ 2 se llega a la siguiente ecuación: 

ndy y -1 lnx 

y I — • 

dx x x 
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dz _ 2 dy 

Ponemos z = y , luego — = —y —. 

dx dx 

multiplicando por —1 se obtiene la ecuación 


Reemplazando estos valores y 


dz z ln x 

dx x x 


Esta ecuación es una ecuación de primer orden lineal cuya solución es 


z = x 




Luego, la solución de la ecuación de Bernoulli está dada por: 


y 1 = x 


ln x 1 
-1-he 

X X 


Ej ercicios 

Los siguientes ejercicios han sido tomados de O’Neil [OP]. 
Resolver las siguientes ecuaciones diferenciales: 


-i dy 4 

i- ~r + y = y ■ 

dx 

2. x 2< ^~ + xy = — y -3 / 2 . 
dx 

o 2 d y , r, 4 

ó. y - — h 2xy = x . 
dx 

a 2 , . d y 

4. xy +y+ —. 

dx 

5. x 3 ^- + x 2 y = 2y~ 4 / 3 . 

dx 


1.5 Ejercicios del Capítulo 1 

El siguiente bloque de ejercicios ha sido tomado de Zill [ZD, Cap. 3] y Makerenko, 
Kiseliov y Krasnov [MK, pag. 34-60], a excepción de los ejercicios 14 y 15. 

1. Resolver las siguientes ecuaciones diferenciales 
a) (1 + y 2 )dx + (1 + x 2 )dy = 0. 
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b) 

c) 

d) 

e) 
/) 
9) 
h ) 
0 

j) 

k) 

o 


y + sen 


x + j/\ ( x — y 

= sen 


ydx + (2y/xy — x)dy = 0. 
x — y eos j j dx + x eos j dy = 0. 


(2x — %)<ix + (x + y — 3 )dy = 0. 


x 


3xy' - 2y = —. 

y 

2 sen (x)y' + ycos(x) = y 3 (xcos(x ) — sen(x)). 
2x + 2y — 1 + y'{x + y — 2) = 0. 

(x — y + 3)dx + (3x + y + l)d?/ = 0 

y' — 2 xy = 2xe x2 . 

. 1 

8xy — y = - . 

y 3 \Jx + 1 

2?/' sen x + y eos x = y 3 (x eos x — sen x). 


2. Hallar un factor de integración de las siguientes ecuaciones 

a) ( x 2 + y)dx — xdy = 0. 

b ) {x + sen x + sen y)dx + eos y dy — 0. 

c) (2x?/ 2 — 3y 3 )dx + (7 — 3 xy 2 )dy = 0. 

d) xdx + ydy + x(xdy — ydx) = 0, y = <¿?(:r 2 + y 2 ). 

e) 3y<¿r + 4xch/ = 0. 

/) Aydx — xdy = 0. 

<?) 3x 2 ydx + ydy = 0. 

h ) x~ 2 y~ 5 dx + x~ 3 y~ A dy = 0, ¡1 = tp(x 3 y 5 ). 


3. Un cultivo tiene una cantidad inicial N 0 de bacterias. Cuando t — 1 hora, la 

3 

cantidad medida de bacterias es -iVo. Si la razón de reproducción es proporcional 
a la cantidad de bacterias presentes, calcule el tiempo necesario para triplicar 
la cantidad inicial de los microorganismos. 
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4. Un reactor de cría convierte al uranio 238, relativamente estable, en plutonio 
239, un isótopo radiactivo. Al cabo de 15 años, se ha desintegrado el 0.043% de 
la cantidad inicial, A 0 , de una muestra de plutonio. Calcule el período medio de 
ese isótopo, si la razón de desintegración es proporcional a la cantidad presente. 

5. Al sacar un pastel del horno, su temperatura es 300 °F. Después de tres minutos, 
200 °F. ¿En cuanto tiempo se enfriará hasta la temperatura ambiente de 70 °F1 

6 . Un acumulador de 12 volts se conecta a un circuito en serie LR, con una 
inductancia de - henry y una resistencia de 10 ohms. Determinar la corriente 
i, si la corriente inicial es cero. 

7. Se sabe que la población de cierta comunidad aumenta con una razón 
proporcional a la cantidad de personas que tiene en cualquier momento. 
Si la población se duplicó en cinco años, ¿en cuanto tiempo se triplicará 
y cuadruplicará? 

8 . Cuando t — 0, había 100 miligramos de una sustancia radiactiva. Al cabo de 
seis horas, esa cantidad disminuyó el 3%. Si la razón de desintegración, en 
cualquier momento, es proporcional a la cantidad de la sustancia presente, 
calcule la cantidad que queda después de dos horas. 

9. Cuando pasa un rayo vertical de luz por una sustancia transparente, la razón 
con que decrece su intensidad / es proporcional a I(t), donde t representa el 
espesor, en pies, del medio. En agua de mar clara, la intensidad, a tres pies 
bajo la superficie, es el 25 % de la intensidad inicial % del rayo incidente. ¿Cuál 
es la intensidad del rayo a 15 pies bajo la superficie? 

10. Un termómetro se lleva de un recinto interior hasta el ambiente exterior, donde 
la temperatura del aire es 5 °F. Después de un minuto, el termómetro indica 
55 °F, y después de cinco marca 30 °F. ¿Cuál era la temperatura del recinto 
interior? 

11. Se aplica una fuerza electromotriz de 30 voltios a un circuito en serie LR con 
0.1 henry de inductancia y 50 ohmio de resistencia. Determine la corriente i(t), 
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s í( 0) = 0. Halle la corriente cuando t —> oo. 


12. Un tanque está parcialmente lleno con 100 galones de salmuera, con 10 libras 
de sal disuelta. Le entra salmuera con - libra de sal por galón a un flujo de 6 
gal/min. El contenido del tanque está bien mezclado y de él sale un flujo de 4 
gal/min de solución. Calcule la cantidad de libras de sal que hay en el tanque 
a los 30 minutos. 


13. Una ecuación diferencial que describe la velocidad v de una masa m en caída 
sujeta a una resistencia del aire proporcional a la velocidad instantánea es 

dv 

m — = mg — kv, 
dt 

donde k es una constante de proporcionalidad positiva. 

a) Resuelva la ecuación, sujeta a la condición inicial u(0) = vq. 

b ) Calcule la velocidad límite (o terminal) de la masa. 

c) Si la distancia s se relaciona con la velocidad por medio de la ecuación 

— = v, deduzca una ecuación explícita para s, si también se sabe que 
dt 

s(0) = s 0 . 

14. Inicialmente había 100 miligramos de una sustancia radiactiva. Después de seis 
horas su masa disminuyó en un 3 %. Si en un instante cualquiera la rapidez de 
desintegración es proporcional a la cantidad de sustancia presente, determinar 
la cantidad que queda después de 24 horas. 


15. La rapidez con que cierto medicamento se disemina en el flujo sanguíneo se 

rige por el siguiente problema con condiciones iniciales 

dx . 

— = A — Bx 
dt 

x(0) = 0, 


donde Ay B son constantes positivas. La cantidad x(t) describe la concentración 
del medicamento en el flujo sanguíneo en un instante cualquiera t. Encontrar 
el valor límite de x cuando t —> oo. ¿Cuánto tarda la concentración en alcanzar 
la mitad de este valor límite? 
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Ecuaciones Diferenciales lineales de orden superior 

En el capítulo 1, se estudió las ecuaciones diferenciales de primer orden y se 
indicó como formular modelos matemáticos para hallar la solución de un problema 
particular (ver sección 1.1). Naturalmente los modelos matemáticos ahí planteados se 
presentaban con el apelativo de problemas con condiciones iniciales donde la ecuación 
diferencial es una ecuación diferencial de primer orden. 

Los modelos matemáticos que se construyen con ecuaciones diferenciales de 
orden superior cubren una amplia gama de situaciones que van desde el problema 
de determinar el movimiento vibratorio de sistemas mecánicos hasta problemas 
relacionados con la cardiología. 

En este capítulo, estudiamos técnicas resolutivas de una ecuación diferencial 
lineal de orden superior y damos algunos ejemplos de aplicación de estas. En la 
sección 2.1, describimos las propiedades básicas de las ecuaciones diferenciales lineales 
de orden superior. En la sección 2.2, indicamos cómo hallar las soluciones de una 
ecuación diferencial lineal de orden superior homogénea con coeficientes constantes. 
En la secció 2.3, presentamos los métodos de coeficientes indeterminados y variación 
de parámetros para hallar una solución particular de la ecuación diferencial lineal de 
orden superior completa con coeficientes constantes. La sección 2.4 (que puede ser 
omitida en un curso introductorio a las ecuaciones diferenciales ordinarias) estudia 
la forma de hallar las soluciones en series de potencias de una ecuación diferencial de 
segundo orden con coeficientes variables. Finalmente, en la sección 2.5, mostramos 
algunos modelos matemáticos con ecuaciones diferenciales lineales de segundo orden. 

2.1 Teoría básica de las ecuaciones diferenciales 

lineales 

Definición 2.1. Una ecuación diferencial lineal de orden n en la variable 
independiente y, y variable dependiente x es una ecuación diferencial que tiene, 
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o se puede escribir, la forma 

a 0 (x)y + a i( x )^, + • • • + a ^ x )~¿ = F ( x )> (2-1) 

donde las funciones üq, ai ,..., a n y F(x) son funciones continuas en un intervalo 
[a, 6] con a n (x) ^ 0 para todo x en el intervalo [a, b\. 


Si, en la ecuación 2.1, la función F es idénticamente igual a cero para todo x en 
[a, b], la ecuación se llama homogénea. En caso contrario, se llama completa. Además, 
en el caso de una ecuación diferencial lineal de n-ésimo orden completa, la función F 
se llama término no homogéneo. 

Observación 2.1. Si, en la ecuación 2.1, las funciones üq, ai, ..., a n son constantes 
en [a, b], la ecuación se llama ecuación diferencial lineal de orden n con coeficientes 
constantes. La ecuación 2.1 se dice que es una ecuación diferencial lineal con 
coeficientes variables cuando al menos una función ai no es idénticamente igual 
a una constante en [a, b]. 


Las siguientes ecuaciones son ejemplos de ecuaciones diferenciales lineales de 
orden superior. 


d 3 y d 2 y 2 x 

5 ~ Td + y = x e 

dx ó dx 2 


es una ecuación diferencia lineal de tercer orden 
completa con coeficientes constantes. 


dfy dy 

2. x 2 -— — 4sena;-—K (. x 3 + 4a; — 1 )y = e x es una ecuación diferencial lineal de 

dx 2 dx 

segundo orden completa con coeficientes variable. 

3. sena;—^ + cosa;— 1 - — y = 0 es una ecuación diferencial lineal de tercer orden 

dx A dx 

homogénea con coeficientes variables. 

4. 2—^r - —\ + 7— — 3a/ = 0 es una ecuación diferencial lineal de cuarto orden 

aa; 4 dx ¿ dx 

homogénea con coeficientes constantes. 


El siguiente teorema es una generalización del teorema 1.1 para el caso de 
ecuaciones diferenciales lineales de n-ésimo orden. Este teorema nos da una condición 
suficiente para la existencia y unicidad de soluciones de un problema con condiciones 
iniciales donde la ecuación diferencial es una ecuación lineal de orden n. 
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Teorema 2.1. Si las funciones üq, a\,..., a n , F son continuas en un intervalo [a, b], 
a n (x) 7 ^ 0 para todo x de [a, 6 ] y xo está en [a,b], entonces existe una solución única 
para el problema con condiciones iniciales 

a 0 (x)y + afix)^- + ... + a n (x)^ = F(x), 
i dx dx n 

y{x o) = y 0 , 

y'(x o) = 2 / 1 , 

< 

y"(x o) = 2 / 2 , 


2/ (n ^(xo) = Vn—l- 

Además, esta solución está definida para todo x en [a, b]. 

Igual que en el capítulo 1, no damos la demostración de este teorema pues esta 
utiliza técnicas que caen fuera del alcance de este texto. 

Los siguientes ejemplos muestran la utilidad de este teorema. 

Ejemplo 2.1. Decidir si el siguiente problema con condiciones iniciales: 

'. o . d 3 y _ dy , 

(x + 1 )-— — 5cosa;-—f- [x — 1 ) y = senx 
dx á dx 

2/(3) = 2 

< 

2/(3) = -1 
/'( 3 ) = 0 

tiene solución única en algún intervalo. 

Solución. Ya que las funciones (x 2 + 1), —5 eos x, (x — l ) 3 y senx son continuas para 

todo x tal que —00 < x < + 00 , x 2 + 1 7 ^ 0 para todo x tal que — 00 < x < 00 y, 

Xo = 3 está en el intervalo (— 00 , + 00 ). Por el teorema 2.1, se puede asegurar que 

existe una única función real / definida en (— 00 , + 00 ) tal que 

2 d 3 f(x) df(x) 3 

x ——-5 cosx— -h (x — 1 ) j(x) = senx 

dx á dx 

para todo x real. Además, esta función satisface las condiciones iniciales, es decir 
/(3) = 2 , / 7 (3) = —1 y /"(3) = 0 . 
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Ejemplo 2.2. Decidir si el problema 


x- 


_1_ <Py_ 

— 1 dx 2 


+ Vx 2 — 25— + (x — 3)y = e 
dx 


2z+3 


2/(5/2) = -3 
2/'(5/2) = 4 
í/"(5/2) = 0 


¿¿ene solución. 


Solución. Se puede demostrar (hacerlo) que las funciones —- , \Jx 2 — 25, x — 3 

x 2 — 1 

y e 2a:+3 son continuas en el intervalo [2,4], Además, —--- 0 para cada x G [2,4], 

5 

Por el teorema 2.1, ya que - G [2,4], se puede asegurar que existe una única 
función real / definida en [2,4] tal que 


1 d 2 f{x ) 

x 2 — 1 dx 2 


+ \Jx 2 - 25 


dx 


+ (x 


3)/(x) = e 2x+3 


y /(5/2) = -3, /'(5/2) = 4 y /"(5/2) = 0. 


Como muestran los ejemplos 2.1 y 2 . 2 , el teorema 2.1 no nos dice como hallar la 
solución a un problema con condiciones iniciales pero, nos dice cuando tiene sentido 
buscar una solución. 


Observación 2.2. Como se indicó en el capitulo 1 (ver observación 1.5), la solución 
que se pueda obtener aplicando algún método resolutivo puede ser una solución 
implícita. A partir del teorema 2.1, podemos estar seguros que esta solución implícita 
determina una función f (solución explícita) definida en algún intervalo I. 


El siguiente corolario es de mucha utilidad, pues nos dice bajo qué condiciones la 
solución de una ecuación diferencial lineal de n-ésimo orden se reduce a la solución 
trivial. 

Corolario 2.1. Sea f una solución en el intervalo cerrado I de la ecuación diferencial 
lineal homogénea de n-ésimo orden 

a 0 (x)y + °i(x)^ + • • • + = 
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Xo un punto de I tal que f(x o) = f'(x o) = f"(x o) = . .. = /^ n_1 ^(xo) = 0. Si las 

funciones ao,ai, ... ,a n son continuas en I, a n (x) 7 ^ 0 para todo x G I y xq G /, 

entonces f(x) = 0 para todo x de I. 

Por ejemplo, la única solución de la ecuación diferencial homogénea de tercer 
orden 

d 3 y 0 du , 

cosa:—- — 5ar-—h 3ary = 0, 
dx ¿ dx 

tal que y(l) = y'(l) — y"( 1 ) = 0 es la función / definida por f{x) =0 para todo x 
real. 


La ecuación lineal homogénea de n-ésimo orden 


Ahora consideramos resultados fundamentales sobre la ecuación diferencial homogénea 
de n-ésimo orden 

a 0 (x)y + a Á x )^ + • • • + = °» ( 2 - 2 ) 

donde las funciones a 0 , ai,..., a n son funciones continuas en un intervalo [a, 6 ] con 
a n (x) 0 para todo x en el intervalo [a, b\. 

El primer resultado básico (cuya demostración es trivial) es el siguiente teorema: 

Teorema 2.2. Si f \, f '2 ■,■■■■, fm son m soluciones arbitrarias de la ecuación 
diferencial homogénea 2.2 y ci,C 2 ,... ,c m son constantes cualesquiera, entonces la 
función C 1/1 + C 2/2 + • • • + c m f m también es una solución de la ecuación 2.2. 


El Teorema 2.2 nos dice cómo crear “nuevas soluciones” de una ecuación diferencial 
homogénea de n-ésimo orden cuando se conocen algunas soluciones de la ecuación. 
Por ejemplo, las funciones sena; y cosa; son soluciones (comprobarlo) de la ecuación 
diferencial homogénea 


d 2 y 
dx 2 


+ y 


0. 


Por el teorema 2.2, la función Ciseno; + C 2 eos x (claramente diferente a las 
funciones seno; y coso;) también es solución de la ecuación diferencial donde, c\ y C 2 
son constantes arbitrarias. 


64 




Leónidas Cerda y Janneth Morocho 


Observación 2.3. Si /j, f 2 , ■ ■ ■, fm son m funciones dadas y c±, C 2 ,..., c m son m 
constantes, entonces la expresión 

Ci/i + C0/2 H-h Cmfm 

se llama una combinación lineal de f\, f 2 ,..., f m . 

En términos de combinaciones lineales, el teorema 2.2 expresa el hecho de que 
cualquier combinación lineal de soluciones de la ecuación diferencial homogénea 2.2 
sigue siendo una solución. 

Nuestro siguiente paso es entender el concepto de solución general de la ecuación 
diferencial homogénea 2.2. Para comprender completamente este concepto debemos 
primero definir los conceptos de dependencia e independencia lineal. 

Definición 2.2. Las funciones f \, / 2 ,..., f m se dicen que son linealmente 
dependientes en el intervalo [a, b] si existen constantes ci, C 2 , ..., c m no todas iguales 
a cero tal que 

cifi(x) + c 2 f 2 (x) + ... + Cmf m {x) = 0 

para todo x en [a,b\. 

Observación 2.4. En la definición 2.2, pedimos que todas las funciones fi, para 
1 < i < m, sean diferentes de la función nula en [a,b\, pues, si alguna de ellas es 
nula, entonces estas son linealmente dependientes. 

Un ejemplo de funciones linealmente dependientes en el intervalo [0,1] son las 
funciones fi y / 2 definidas por fi(x) = x y f 2 (x) = 2x pues existen las constantes 
Ci = 2 y C 2 = —1 tal que 

C 1 / 1 O) + c 2 f 2 (x ) = 2x - 2x = 0 
para todo x en el intervalo [ 0 , 1 ]. 

Otro ejemplo de funciones linealmente dependientes son las funciones sena;, 
3 sen x y — sen x que resultan ser linealmente dependientes en cualquier intervalo 
[a, b] pues existen las constantes ci = 1, C 2 = 1 y C 3 = 4 tal que 

Ci seno; + c 2 (3sena;) + c 3 (— sena;) = 0 
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para todo valor de x en [a, b\. 

Observación 2.5. Desde los ejemplos anteriores, se podría pensar que si las 
funciones fi, f 2 , ■ ■ ■, fm son linealmente dependientes, entonces una de ellas se puede 
expresar como combinación lineal de las otras. Dejamos como ejercicio para el lector 
demostrar que esta intuición es correcta. 

Definición 2.3. Las funciones fi, f 2 , ■ ■ ■, fm se dicen que son linealmente 
independientes en el intervalo [a, b] si no son linealmentes dependientes en el intervalo 



Notamos que las funciones fi, f 2 , ■ ■ ■ , fm son linealmente independientes en el 
intervalo [a, b] cuando se cumple la siguiente condición: Si la ecuación 

cifi(x) + c 2 f 2 (x ) + ... + Cmf m {x) = 0 

se cumple para todo x en el intervalo [a,b\, entonces C\ — C2 — ■ ■ ■ — c m — 0. En 
otras palabras, la única combinación lineal de _/j, f 2 , ■ ■ ■, fm que es idénticamente 
igual a cero en el intervalo [a, b] es la combinación lineal trivial 


0 • fi(x) + 0 ■ f 2 {x) + ... + 0 ■ f m (x) = 0. 


Ejemplo 2.3. Averiguar si las funciones x y x 2 son linealmente independientes en 
el intervalo [0,1]. 

Solución. Para ver que estas funciones son linealmente independientes en el intervalo 
[0,1] debemos probar lo siguiente: Si cpx + c 2 x 2 = 0 para todo x en el intervalo [0,1], 


entonces C\ = c 2 = 0. En efecto, a partir de cix-j-c 2 x 2 = 0 para todo x en el intervalo 


[0,1] tenemos (derivando ambos lado de la ecuación) c± + 2c 2 x = 0 para todo x en el 
intervalo [0,1], luego C\X + 2c 2 x 2 = 0 para todo x en el intervalo [0,1]. Así, tenemos 
el siguiente sistema de ecuaciones 
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que debe satisfacerse para todo x en el intervalo [0,1]. Restando la segunda ecuación 
de la primera se tiene que C 2 X 2 = 0 todo x en el intervalo [0,1], luego C 2 = 0. De 
manera similar se tiene que C\ = 0. 

El siguiente teorema asegura la existencia de un conjunto de soluciones linealmente 
independientes de la ecuación diferencial homogénea 2.2 y cual es el significado de 
tales conjuntos linealmente independientes. 

Teorema 2.3. La ecuación diferencial homogénea de n-ésimo orden 2.2 siempre tiene 
n soluciones linealmente independientes. Además, si f \, / 2 ,..., f n son n soluciones 
linealmente independientes de la ecuacón diferencial homogénea de n-ésimo orden 
2.2, entonces cualquier solución f de la ecuación 2.2 puede expresarse como una 
combinación lineal c\fi + C 2/2 + ... + c n f n de estas n soluciones linealmente 
independientes. 

Para ver cómo se puede aplicar el teorema 2.3 consideremos las funciones definidas 
por fi(x) = senx y /^(x) = eos x. Como vimos anteriormente (ver pag. 64) estas 
funciones son soluciones de la ecuación diferencial homogénea de segundo orden 



(2.3) 


para todo x tal que — 00 < x < + 00 . Además, estas funciones son linea l m ente 
independientes para todo x en el intervalo (— 00 , + 00 ) (demostrarlo). Luego, si una 
función / es solución de 1a. ecuación diferencial 2.3, entonces, por el teorema 2.3, / 
puede ser expresada como combinación lineal de las funciones sen x y eos x, esto es 
f{x) = C\ sen x + c 2 cosx para valores adecuados de las constantes Ci y c 2 . Además, 
esta igualdad se cumple para todo x en el intervalo (— 00 , + 00 ). 

De esta forma, por ejemplo, 1a. función definida por 


sen (x + — J 


7r 

sen x ■ eos —h eos x ■ sen 
6 


7T 


6 / 


6 


es solución de la ecuación 2.3. 

Observación 2.6. Sean f \, f-¿,, f n un conjunto de n soluciones linealmente 
independientes de la ecuación diferencial 2.2. Por un lado, desde el teorema 2.2 
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conocemos que la conbinación lineal 

c ifi + C 2/2 + • • • + c n f n , (2-4) 

donde Ci,c 2 , ■ ■ ■ ,c n son constantes arbitrarias, es una solución de la ecuación 2.2. 
Por otro lado, por el teorema 2.3 conocemos que si f es cualquier solución de la 
ecuación diferencial 2.2, entonces esta puede ser expresada como la combinación 
lineal 2-4 de las n soluciones linealmente independientes f \, / 2 ,..., f n para una 
adecuada elección de las constantes C\, c 2 ,..., c n . Así, una combinación lineal de las 
n soluciones linealmente independientes / 1 , / 2 , ■ ■ •, fn con c\, C 2 ,..., c n constantes 
arbitrarias debe contener todas las soluciones de la ecuación diferencial 2.2. Por esta 
razón, nos referimos al conjunto de n soluciones linealmente independientes como 
“un conjunto fundamental” de soluciones de la ecuación diferencial 2.2 y llamamos 
solución “general” de la ecuación diferencial 2.2 a una combinación lineal de las n 
soluciones linealmente independientes. 

A continuación se define rigurosamente lo que entendemos como conjunto 
fundamental de soluciones de una ecuación diferencial homogénea de orden n y 
cual es su solución general. 

Definición 2.4. Sean f\, f 2 ,..., f m soluciones de la ecuación diferencial homogénea 
de n-ésimo orden 2.2 en el intervalo [a, b\. El conjunto A = {/ 1; / 2 ,..., f m } se llama 
conjunto fundamental de soluciones de la ecuación 2.2 en el intervalo [a,b\ si y solo 
si se cumplen las siguientes condiciones: 

1. A es un conjunto maximal de funciones linealmente independientes. 

2. Cualquier solución de la ecuación 2.2 se puede expresar como 

Clfl + c 2 f 2 + • • • + Cmfmi 


donde C\, c 2 ,..., c m son constantes arbitrarias. 

En la definición 2.4, se entiende como conjunto maximal de funciones linealmente 
independientes a aquel conjunto que contiene el máximo número de funciones 
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linealmente independientes, es decir, si B es un conjunto de funciones linealmente 
independientes, entones B C A. Además, la función definida como la combinación 
lineal de las funciones /i,...,/ m de este conjunto maximal A se llama solución 
general de la ecuación diferencial 2.2 en el intervalo [a,b\. 

Observación 2.7. Por el teorema 2.2, la solución general de la ecuación diferencial 
2.2 es efectivamente una solución. Además, por el teorema 2.3 y la observación 2.5, 
un conjunto maximal de soluciones linealmente independientes de la ecuación 2.2 
tiene exactamente n elementos. 


Es fácil verificar que las funciones e x ,e x ,e 2x son soluciones de la ecuación 


diferencial 


d 3 y d 2 y dy 

— - 2—-- + 2y = 0 

dx 3 dx 2 dx y 


(2.5) 


para todo x tal que —oo < x < +oo. Además, se puede demostrar que estas 


funciones son linealmente independientes en el intervalo (—oo, +oo). Luego, según la 
definición 2.4 y la observación 2.7, el conjunto {e x , e ~ x , e 2x } es el conjunto fundamental 
de soluciones de la ecuación 2.5 y su solución general es la función definida por 

c 1 e x + c 2 e~ x + c 3 e 2x . 


Observación 2.8. Saber que un conjunto de n funciones es linealmente independiente 
es de transcendental importancia, tanto para conocer cual es el conjunto fundamental 
de soluciones, como para determinar cual es la solución general de una ecuación 
diferencial. 


Los teorema 2.4 y 2.5 (ver más abajo) nos dan un criterio sencillo para 
averiguar cuando un conjunto de soluciones de una ecuación diferencial es linealmente 
independiente en un intervalo [a, b\. Primero necesitamos dar la siguiente definición. 


Definición 2.5. Sean /i, f 2 , ■ ■ ■, f n , n funciones reales que tienen derivadas hasta 
de orden n — 1 en un intervalo real [a, b]. El determinante 


W(f h f 2 ,...J n ) = 


fi h 

f fk 


fn 

fn 


yfii- 1 ) fpi-l) ff'- 1 '* 
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se llama Wronskiano de estas n funciones. 


Observamos que el Wronskiano de /i, / 2 , ..., f n es una función real en el intervalo 
[a, b], su valor lo denotamos por W(fi, / 2 , ■ ■ ■, fn)(x). 

Ejemplo 2.4. Hallar el Wronskiano de las funciones definidas por fi(x) = e x , 
h{x) = e~* V h(x) = e 2x . 


Solución. Desde la definición 2.5 tenemos 


W(e x , e~ x , e 2x ) = 


X „—x 


e e 


-2x 


e x —e x 2e 


2x 


x „—x 


e e 


4e : 


2x 


luego W(e x ,e x ,e 2x )(x ) = — 6e 2x . 


Ejemplo 2.5. Hallar el Wronskiano de las funciones definidas por f\(x) = sena; y 
/ 2 (a;) = 2 sena;. 


Solución. En este caso tenemos 


IU(sena;, 2 sen a:) 


sen x 2 sen x 
eos x 2 eos x 


luego VE (sen a;, 2 sen) (a;) = 0. 

Los siguientes teoremas nos permiten averiguar si un conjunto de soluciones 
de una ecuación diferencial de n-ésimo orden es linealmente independiente en un 
intervalo. 


Teorema 2.4. Sean f \, / 2 ,..., f n soluciones de la ecuación diferencial lineal 
homogénea de n-ésimo orden 2.2. Las soluciones f \, / 2 ,..., f n son linealmente 
independientes en el intervalo [a,b\ si y solo si VE(/j, / 2 ,..., f n )(x) es diferente 
de cero para algún x en [a, b]. 

Teorema 2.5. El Wronskiano de n soluciones f\, f 2 ,, f n de la ecuación diferencial 
2.2 es, o bien idénticamente igual a cero en [a, b] o bien diferente de cero en [a,b\. 
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Observación 2.9. Los teoremas 2-4 y 2.5 no tienen validez (al menos en todos los 
casos) cuando se elimina la hipótesis que pide que las funciones sean soluciones de la 
ecuación diferencial. 


Para una discución general de la relación entre dependencia lineal de un conjunto 
de funciones y el Wronskiano asociado a esas funciones ver, por ejemplo, Bóclier 
[BM] o Courant [CR, pag. 756-763]. 

Ejemplo 2.6. Decir si las funciones definidas por f\(x) = sen x y /^(x) = 2sena; 

d 2 y 

forman un conjunto fundamental de soluciones de la ecuación —- + y — 0. 

dx 2 

Solución. Primero vemos que las funciones f\ y / 2 son soluciones de la ecuación 
diferencial. En efecto: 


/l(z) 

fí(x) 

m*) 

hix) 

/ 2 O) 

/£(*) 

seno; 

cosx 

— seno; 

2 sen x 

2 cosx 

—2 senx 


Reemplazando estos valores, vemos claramente que fi y /2 son soluciones de la 
ecuación diferencial. En el Ejemplo 2.5, se vio que lE(senx, 2 sen a;) (a;) = 0. Luego, 
por el teorema 2.4 y el teorema 2.5, estas funciones son linealmente dependientes. 
Por lo tanto, no constituyen un coniunto fundamental de soluciones de la ecuación 

diferenda , pL + y = 0 . 

dx ¿ 

Ejemplo 2.7. Conociendo que las funciones seno; y cosx son soluciones de la 
d 2 y 

ecuación diferencial + y — 0. Hallar la solución general de la ecuación diferencial. 


Solución. El Wronskiano de las funciones seno; y coso; es 


W (coso;, seno;) = 


sen x eos x 
eos x — sen x 


= — sen 2 x — eos 2 x 


= - 1 . 


Luego, por los teoremas 2.4 y 2.5 se tiene que seno;, cosx son linealmente 
independientes en todos los reales. Luego (ver observación 2.7) el conjunto 
{sen x, eos x} es el conjunto fundamental de soluciones de la ecuación diferencial 
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(j~ y 

+ y = 0 y su solución general es la función definida por f{x) 
Además, el intervalo solución es M. 


c\ sen x + C 2 eos x. 


Ejemplo 2.8. Hallar la solución general de la ecuación 

d 3 y d 2 y dy 

— - 2—-- + 2y = 0. 

dx 3 dx 2 dx y 


Solución. Derivando, se demuestra que las funciones e x , e~ x y e 2x son soluciones 
de la ecuación diferencial. Como W(e x , e~ x , e 2x )(x) = —6e 2x ^ 0 para todo x real. 
Tenemos que e x , e~ x y e 2x son funciones linealmente independientes para todo x tal 
que —oo < x < +oo. Luego, la solución general de la ecuación de este ejemplo está 
definida por 

f(x) = cxe x + c 2 e~ x + c 3 e 2:r . 


Reducción de orden 

Reducir el orden de una ecuación diferencial lineal homogénea de n-ésimo orden 
es una técnica que facilita mucho el trabajo de hallar soluciones de una ecuación 
diferencial. Su aplicación depende del conocimiento de una solución no trivial de la 
ecuación. 

El próximo teorema nos indica c’omo podemos reducir el orden de una ecuación 
diferencial. 

Teorema 2.6. Si f es una solución no trivial de la ecuación lineal homogénea de 

n-ésimo orden 2.2, entonces la transformación y = fv reduce la ecuación 2.2 a una 

dv 

ecuación diferencial lineal homogénea de (n — l)-ésimo orden en la variable w = —. 
Además, la función definida como f(x)v(x) también es solución de la ecuación 2.2. 

Para ejemplificar la utilidad del teorema 2.6, realizamos un estudio completo de 
la ecuación de segundo orden. 

Sea / una solución no trivial de la ecuación diferencial lineal homogénea de 
segundo orden 

a Q (x)y + ai(x)^- + +a 2 (x)^ = 0 . (2.6) 

dx dx z 
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Utilizamos la transformación 


y = 


(2.7) 


El objetivo es determinar la función v. Diferenciando obtenemos 

dy = + ^df_ 

dx ax dx' 

d 2 y d 2 v df dv df 2 

dx 2 dx 2 dx dx + 1 dx 2 ' 

Sustituyendo estos valores en la ecuación 2.6 obtenemos 

. . ( dv df\ . f d 2 v df dv df 2 \ 
a „(x)f v + a,(x) + ■” Tx j + “= M 

Ya que / es una solución de la Ecuación 2.6, la última ecuación se reduce a la ecuación 

a 2( x )fi x )^ + (a^fix) + 2a 2 {x)^j ^ = 0 . 

Si ponemos w = — , entonces la última ecuación se puede escribir como 
dx 


a 2 (x)f(x)P + ía¡(x)f(x) + 2a 2 (x)f- ) w = 0. 


( 2 . 8 ) 


La ecuación 2.8 es una ecuación diferencial en variables separables cuya solución es 


w = 


_ f 

ce J a 2 ( x ) 

I75F' 


Sin pérdida de generalidad, podemos hacer c = 1. Luego se tiene que 


v = 


_ r 

e J “ 2 (*> 

JfW 


dx. 


Finalmente, una solución de 1a. ecuación diferencial 2.6 está dada por 


y = /O) 


_ríiM dx 
£ J a 2\ x ) 

JfW 


dx. 


(2,9) 


Observación 2.10. Se puede demostrar (hacerlo) que W(f,g) 0, donde f es una 


solución no trivial de la ecuación 2.6 y g = f{x) 


_ r ttl(T) , 
g J a 2 (x) 


dx. Por tanto, la 


(fW 

solución general de la ecuación lineal homogénea de segundo orden esta dada por 

cif + c 2 g- 
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Ejemplo 2.9. Se puede comprobar fácilmente que la función definida por f(x) = e x 

(f 2 y 

es solución de la ecuación diferencial —— — y = 0. Hallar una solución que sea 

J dx 2 y H 

linealmente independiente con la solución f. 


Solución. Según la observación 2.10, reemplazando los valores de las funciones ai, Ü 2 

d 2 y 

y / en la ecuación 2.9, una solución de la ecuación diferencial —— — y = 0 que 

dx 1 

es linealmente independiente con la función definida por f{x) = e x está dada por 


g(x) = e 1 


e 2x dx. Luego g(x) = — -e 


Ejemplo 2.10. Se sabe que la función definida por f(x) = x 2 es solución de la 


2 d 2 y 2 n _dy 


ecuación x 2 — f-x 2 — 3x ——|- 4y = 0. Hallar otra solución de la ecuación diferencial, 
d dx y J 


Solución. En este caso se tiene 


ai(g) 

a 2 {x) 


dx= —dx 

J x 

= — 31nx. 


Luego 



= lnx. 


2 d 2 y 2 n ~dy 


De esta forma, una solución de la ecuación deferencial x - x — 3x— —h 4y = 0 es 

d dx 

la función definida por g(x) = x 2 lnx. 


La ecuación lineal no homogénea de n-ésimo orden 


Para terminar esta sección enunciamos algunos resultados básicos sobre la ecuación 
diferencial lineal de n-ésimo orden no homogénea o completa 

a 0 (x)y + afix) 1 ^- + ... + a n (x )^-| = F(x). (2.10) 

dx dx n 

E1 resultado básico es el teorema que relaciona la solución de la ecuación completa y 
su correspondiente ecuación homogénea asociada. Primero definimos el concepto de 
ecuación homogénea asociada. 
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Definición 2.6. Se llama ecuación homogénea asociada a la ecuación 2.10 a la 
siguiente ecuación: 

a 0 (x)y + a Á x )^ + • • • + a ^ x )~¿ = °- (2- 11 ) 

Ahora enunciamos el teorema que relaciona la solución de la ecuación 2.10 con la 
solución de la ecuación 2.11. 

Teorema 2.7. Si f, g son soluciones arbitrarias de las ecuaciones 2.10 y 2.11 
respectivamente, entonces la función f + g es solución de la ecuación 2.10. 


En efecto, calculando las derivadas hasta el orden n de la función / + g se tiene 


que 


a 0 (x)(f + g) + ai(x) d ^ + ^ + ... + a n (x) d ^ ^ = F(x ) + 0 


dx 


dx r 


= F{x). 

Por ejemplo, la función definida como f(x) = x es solución de la ecuación 

Cpy 

-j—^ + y = x (comprobarlo) y la función dada por g(x ) = sen x es una solución de la 

d 2 y 

ecuación homogénea asociada —— + y = 0 (comprobarlo). Luego, por el teorema 2.7, 

dx 1 

Cpy 

la función definida por x + sen x es solución de la ecuación —— + y = x. 

dx 2 

Ahora, aplicando el teorema 2.7 al caso especial de tener una solución sin 
parámetros de la ecuación completa y la solución general de su ecuación homogénea 
asociada tenemos el siguiente corolario: 


Corolario 2.2. Si y p es una solución particular de la ecuación 2.10, y c es la solución 
general de la ecuación 2.11, entonces cualquier solución (p de la ecuación 2.10 está 
dada por <p = y p + y c . 


El corolario 2.2 sugiere la siguiente definición: 

Definición 2.7. Consideremos la ecuación diferencial no homogénea (ver ecuación 
2.10) y su ecuación diferencial homogénea asociada (ver ecuación 2.11). 


1. La solución general de la ecuación homogénea asociada se llama función 
complementaria de la ecuación no homogénea. A esta solución la notaremos 
por y c . 


75 






INTRODUCCIÓN A LAS ECUACIONES DIFERENCIALES ORDINARIAS 


2. Cualquier solución de la ecuación no homogénea que no contenga parámetros 
se llama integral particular de la ecuación no homogénea. A esta solución la 
notaremos con y p . 

3. La solución y c + y p , donde y c es la función complementaria, y p es una integral 
particular de de la ecuación no homogénea, se llama solución general de la 
ecuación no homogénea. 

d 2 y 

Ejemplo 2.11. Hallar la solución general de la ecuación —- + y = x. 

dx 2 


Solución. La función y p = x es una integral particular de la ecuación. Además, 
la función y c = Ci sen x + C 2 eos x es la función complementaria de la ecuación 


diferencial, luego x + C\ sen x + c^ eos x es la solución general de la ecuación diferencial 
d 2 y 

—| + y = x. 
dx 2 


2.2 Ecuaciones lineales homogéneas con 

coeficientes constantes 


Hallar la solución general de una ecuación diferencial lineal homogénea de n-ésimo 
orden con coeficientes variables es una tarea muy complicada de realizar y, en la 
mayoría de los casos es simplemente imposible hallar tal solución. Sin embargo, si la 
ecuación tiene coeficientes constantes la tarea se vuelve bastante fácil. 

dy d n y 

En efecto, sea a^y + a\— + ... + a n - — = 0, donde a n ^ 0 una ecuación 

dx dx n 

diferencial homogénea de n-ésimo orden con coeficientes constantes. Mirando la 
ecuación diferencial nos damos cuenta de que una solución de esta, se podría pensar, 
es una función / que tiene la siguiente propiedad: “/ y sus derivadas son múltiplos 
de sí mismo”. Ahora, ¿conocemos alguna función que cumpla esta propiedad? La 
respuesta es SI, la función exponencial e mx , donde m es una constante, es tal que 

Jk ( p mx\ 
u 


dx k 


= m k e mx . De esta forma podríamos decir que la solución de la ecuación 


diferencial homogénea de n-ésimo orden con coeficientes constantes tiene la forma 
y = e mx donde la constante m se debe escoger de tal forma que se satisfaga la 
ecuación. 
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Ya que 


f(x) = e mx , 

df(x) mx 
= me , 
dx 

d 2 f(x ) 2 

—= jn e 
dx- 


d n f(x ) 
dx 11 


= m n e mx . 


Reemplazando estos valores en la ecuación diferencial se tiene que 


a 0 e mx + ai me mx + a 2 m 2 e mx + ... + a n m n e mx = 0 


o 

e mx (a 0 + a\m + a 2 m 2 + ... + a n m n ) = 0 . 

Como e mx 7 ^ 0 para todo valor de x se tiene, desde la última expresión, que 

a 0 + Gqm + a 2 m 2 + ... + a n m n = 0 . ( 2 - 12 ) 

Podemos concluir que: si la solución de la ecuación diferencial homogenéa de n-ésimo 
orden con coeficientes constantes tiene la forma f(x) = e mx , entonces la constante m 
debe satisfacer la Ecuación 2.12. 

El objetivo de esta sección es presentar un método explícito para hallar la solución 
general de la ecuación diferencial lineal homogénea de n-ésimo orden con coeficientes 
constantes: 

dy d n y 

&o V + + • • • + ®nr 0- (2-13) 

El método que vamos a exponer está fundamentado por lo expuesto al inicio de esta 
sección. Antes de exponer el método, necesitamos dar la siguiente difinición: 

Definición 2.8. Se llama ecuación característica asociada a la ecuación 2.13 a la 
siguiente ecuación algebraica: 

ao + a±7n + a 2 m 2 + ... + a n m n = 0. (2.14) 

Notamos que la potencia de la variable m en la ecuación característica 2.14 está 
determinada por el orden de la derivada de y en la ecuación diferencial 2.13, donde 
la derivada de orden cero de y coincide con la función y. 

77 







INTRODUCCIÓN A LAS ECUACIONES DIFERENCIALES ORDINARIAS 


Ejemplo 2.12. La ecuación diferencial 2y — 5 


dy d 3 y 
dx dx 3 


0 tiene ecuación 


característica asociada 2 — 5 m — m 3 = 0. 

Existe una correspondencia unívoca entre las ecuaciones diferenciales lineales con 
coeficientes constantes y las ecuaciones algebraicas. De esta forma, a cada ecuación 
diferencial de la forma 2.13 le corresponde una ecuación algebraica de la forma 2.14 
y viceversa, a cada ecuación algebraica de la forma 2.14 le corresponde una ecuación 
diferencial de la forma 2.13. 

Ejemplo 2.13. Hallar la ecuación diferencial lineal con coeficientes constantes 
asociada a la ecuación algebraica 3 — 5 m 2 + 2 m 3 + 4m 5 = 0. 

Solución. Lo único que se tiene que hacer es reemplazar la potencia de m por el 
orden de la derivada de la ecuación diferencial, recordando que la potencia “0” de m 
corresponde a la derivada de orden cero en la ecuación diferencial. Así, la ecuación 



diferencial que buscamos es 3y — 5 


dx 2 dx 3 dx 5 


Se puede demostrar que: la forma de las soluciones de la ecuación diferencial 
homogénea 2.13 dependen de la naturaleza de las soluciones de su ecuación 
característica asociada 2.14. 

Ya que una ecuación algebraica de n-ésimo grado tiene exactamente n raíces y el 
conjunto fundamental de soluciones de una ecuación diferencial lineal con coeficientes 
constantes tiene exactamente n elementos, podemos concluir que existe una solución 
fundamental de la ecuación 2.13 por cada raíz de la ecuación 2.14. Además, las raíces 
de cualquier ecuación algebraica pueden ser: 

a. reales simples, 

b. reales múltiples, 

c. complejas simples y, 

d. complejas múltiples. 

Así, se tiene cuatro casos, dependiendo de que tipo de solución se tenga para la 
ecuación característica asociada 2.14. 
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Caso raíces reales simples 


Si es una raíz real simple de la ecuación característica asociada 

a 0 + aim + a 2 m 2 + ... + a n m n = 0, 

entonces la función definida por e mi es una solución de la ecuación diferencial 

dy d n y 

a oV + a i ~,—!■■■■ + a n~¡—~ = 0. 

dx dx n 

Observación 2.11. Si son raíces reales simples y diferentes (para if^j) de 

la ecuación característica asociada, entonces las funciones definidas pore mi ,e m:i son 
soluciones linealmente independientes de la ecuación diferencial. 


Cuando la ecuación característica asociada tiene n raíces reales simples, digamos 
• • • ,rn n , la solución general de la ecuaión diferencial homogénea 2.13 está 
dada por: 

y(x) = Cl e mix + c 2 e m2X + ... + c n e m ^, 
donde ci, c 2 ,..., c n son constantes arbitrarias. 


Ejemplo 2.14. Hallar la solución general de la siguiente ecuación diferencial: 


„ d A y „ „ d 3 y 
6 d¡rí ~ 


29^ + 52^ + 20t/ 

dx z dx 


0 . 


Solución. La ecuación característica asociada a esta ecuación diferencial homogénea 
es: 

6777 4 — 13m 3 — 29?n 2 + 52 m + 20 = 0. 

Las raíces de la ecuación característica asociada son m — 2, m — —2, m = — 

’ ’ 3 

5 

y m = - (hacer los cálculos). Luego, las soluciones fundamentales de la ecuación 
diferencial son e 2x , e~ 2x , e~^ x y e^ x . Así, la solución general requerida es: 

y(x) = c\e 2x + c 2 e -2;c + c 3 e~^ x + c A e^ x . 


Ejemplo 2.15. Hallar la solución general de la ecuación 
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Solución. La ecuación 15m 4 — 67m 3 +51m 2 +43m+6 = 0 es la ecuación característica 

asociada a la ecuaión diferencial. Ya que m — 2, m — 3, m = — y rn — — son 

3 5 

las soluciones de la ecuación asociada, se tiene que las soluciones fundamentales de 
la ecuación diferencial son e 2x , e 3x , e~^ x y e~^ x . Luego, la solución general de la 
ecuación diferencial está dada por 


y(x ) = c\é ¿x + c 2 e 3x + c 3 e + c 4 e ¿ 


Caso raíces reales múltiples 


Si rrij es una raíz real con multiplicidad a de la ecuación característica asociada 2.14, 
entonces é mjX ,xé mjX ,... ,x a ~ 1 é mjX son a soluciones linealmente independientes de la 
ecuación diferencial 2.13. 


Ejemplo 2.16. Hallar la solución general de la siguiente ecuación diferencial: 


dfy „ dfu „ d 3 y _ d 2 u ^dy 

d^~ 2 d^~ 6 d^ + 2 °d^~ 19 ár + V 


0 . 


Solución. La ecuación característica asociada a la ecuaión diferencial es: 


m 5 — 2 m 4 — 6 m 3 + 20 m 2 — 19m + 6 = 0. 


Las raíces de esta ecuación son m — 2, m — —3 y m — 1 donde las dos primeras 
soluciones son simples (sin multiplicidad) y la tercera solución tiene multiplicidad 3 
(hacer los cálculos). 

Según el caso de raíces reales simples, se tiene que las funciones definidas por 
e 2x y e~ 3x son soluciones de la ecuación diferencial. Por otro lado, según el caso de 
soluciones reales múltiples se tiene que las funciones dadas por e x , xe 2 y x 2 e x son las 
tres soluciones asociadas a la solución múltiple m — 1. Luego la solución general de 
la ecuación diferencial está dada por: 

y(x ) = Cie 2x ’ + c 2 e~ 3x — c 3 e x + c 4 xe x + c^x 2 e x . 


Ejemplo 2.17. Hallar la solución general de la ecuación 


4Í|-20^ + 37^ 

dx a ax 4 dx d 


-30 


d 2 y 

dx 2 



0 . 
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Solución. La ecuación característica asociada a la ecuación diferencial es 

4m 5 — 20 m 4 + 37 m 3 — 30m 2 + 9 m = 0. 

3 

Las raíces de esta ecuación son m = 0, m = - y m = 1 donde, la primera solución 
es simple y las dos últimas soluciones tienen multiplicidad 2 cada una (realizar los 
cálculos). Luego, la solución general de la ecuación diferencial está dada por la función 
definida como 

3 3 

y(x ) — ci + c 2 e^ x + c 3 xez x + c 4 e x + c 5 xe x . 

Caso raíces complejas 

Si rrij = a + fii es un a raíz compleja simple de la ecuación característica asociada 
2.14, entonces las funciones definidas por e ax sen(/3x), e ax cos(fix) son soluciones 
linea l m ente independientes de la ecuación diferencial 2.13. 

Observación 2.12. Sirrij = a+fii es solución de la ecuación característica asociada, 
entonces rñj = a — /3i también es solución de la ecuación caracteréstica. Esta solución 
rrTj no genera nuevas soluciones linealmente independientes (¿por qué?), luego es 
suficiente tomar en cuenta la solución rrij para construir la solución general de la 
ecuación diferencial 2.13. 

Ejemplo 2.18. Hallar la solución general de la ecuación diferencial: 

f!-5^+ 3 j!|- 1 5^- 4 ? + 20 ! , = 0 . 

dx° ax 4 ax á dx z dx 

Solución. La ecuación característica asociada a esta ecuación diferencial es: 
m 5 — 5 m 4 + 3 m 3 — 15 m 2 — 4 m + 20 = 0. 

Las raíces de esta ecuación algebraica son rn — 1, m — 2 + i, m — 2 — i, m — 2i, 
m = — 2 i (comprobarlo). Además, todas las raíces son simples. 

Las funciones e x (asociada a la raíz m — 1), e 2x senx, e 2x cosx (asociadas a la 
raíz m = 2 + i), sen( 2 x), cos( 2 x) (asociadas a la raíz m = 2 i) son las soluciones 
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linealmente independientes de la ecuación diferencial. Así, la solución general de la 
ecuación diferencial está dada por: 

y(x) = c\e x + C 2 e 2x sena; + c 3 e 2x cosa; + C 4 sen( 2 a;) + C 5 cos( 2 a;). 

Ejemplo 2.19. Hallar la solución general de la ecuación diferencial: 

dfy d 3 y d 2 y 

dx 5 dx 3 dx 2 ^ 

Solución. La ecuación característica asociada a esta ecuación diferencial es: 

m 5 — m 3 — m 2 + 1 = 0 . 

1 \/3 

Las raíces de esta ecuación algebraica son m = 1, m = —1, m — - \——i y 

1 a/3 

m = — - - —i (comprobarlo). La primer raíz tiene multiplicidad 2 y las restantes 

raíces son simples. 

Las funciones 


e x , xe x asociadas a la raíz múltiple m — 1 , 


e x asociada la la raíz simple m = — 1 , 

e~i x sen ; e ~* x cos I asociadas a la raíz compleja 

1 , a/3. 

m = -- % 

2 2 


son las soluciones linealmente independientes de la ecuación diferencial. Así, la 
solución general de la ecuación diferencial está dada por: 

y{x) = cie x + c 2 xe x + c 3 e~ x + c 4 e~^ x sen j ^ x J + c 5 e~¿ x cos | ~^~ x I • 


Caso raíces complejas múltiples 

Si nij = a + (3i es una raíz compleja con multiplicidad 7 de la ecuación característica 
asociada 2.14, entonces las 2y funciones definidas como e ax sen(/3x), e ax cos(/3x), 
xe ax sen (/3x), xe ax cos (/3x ),..., x 7_1 sen(,da;), x 7 ^ 1 cos (/3x) son soluciones linealmente 
independientes de la ecuación diferencial 2.13. 
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Observación 2.13. Si la ecuación característica asociada tiene una raíz compleja 
mj = a + /3i con multiplicidad 7 , entonces los 7 números complejos mj = a — f3i 
también son raíces de la ecuación característica asociada. Estas soluciones no generan 
nuevas soluciones linealmente independientes (ver observación 2.12), por tanto, es 
suficiente tomar en cuenta las soluciones mj. 


Ejemplo 2.20. Hallar la solución general de la ecuación diferencial: 


d 6 y d 5 y d 4 y 
-4 - 10—4 + 59—4 
dx° dx° dx 4 


- 196 


d 3 y 
dx 3 


+ 419 


d 2 y 

dx 2 


442^ + 169 y 
dx 


0 . 


Solución. La ecuación característica asociada es: 


m 6 — 10m 5 + 59m 4 — 196m 3 + 419m 2 — 442m + 169 = 0. 


Las soluciones de esta ecuación son m=l,ffl = 2 + 3iym = 2- 3i todas con 
multiplicidad 2 (verificarlo). Por un lado, las funciones definidas por e x , xe x son 
soluciones linealmente independientes asociadas a la solución real múltiple m — 1. 
Por otro lado, las funciones e 2x sen(3x), e 2x cos(3a;), xe 2x sen (3a;) y xe 2x cos(3x) son 
las 4 soluciones linealmente independientes asociadas a la raíz compleja múltiple 
m = 2 + 3i (la raíz m — 2 — 3i no se toma en cuenta - ver observación 2.12). 

De esta forma, la solución general de la ecuación diferencial es: 
y(x) = cie x + c 2 xe x + C3 e 2xsen( ' 33: + c 4 e 2x cos(3x) + c 5 a:e 2:rsen( ' 3;c + c e xe 2x cos(3x). 


Ejemplo 2.21. Hallar la solución general de la siguiente ecuación diferencial 

d 4 y , A d z y d 2 y dy 

4 x? +i d^ + 13 i^ +6 x: +9y - 0 ‘ 


Solución. La ecuación característica asociada es 4m 4 + 4m 3 + 13m 2 + 6 m + 9 = 0. 

Las soluciones de esta ecuación son m = — - + ¿ y m — -—-—z cada una 

4 4 4 4 

con multiplicidad 2 . Así, la solución general de la ecuación diferencial es: 


y{x) — c\e 4X sen 


\/23 


-x 


c 3 e 4 eos 


V23 


+ c 2 xe 41 sen 


V23 


-x 


x + C/íXe 4X eos 


V23 
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Ej ercicios 

En los siguientes ejercicios hallar la ecuación diferencial lineal homogénea con 
coeficientes constantes cuyo polinomio característico asociados tiene las siguientes 
raíces: 


1. Ai = 2 raíz simple, A 2 = —3 raíz con multiplicidad 3, A 3 = 2 + 3* raíz simple, 
A 4 = 2 — 3 i. 

2. Ai = 2/5 raíz simple, A 2 = —1 raíz simple, A 3 = 0 raíz con multiplicidad 
2, A 4 = —3/4 raíz simple, A 5 = 1 + 2 i raíz simple, A 6 = 1 — 2 i. 

3. A = 1 raíz con multiplicidad 5. 

4. Ai = 3 + 2 i raíz con multiplicidad 2, Ai = 3 — 2 i raíz con multiplicidad 2. 


Hallar la solución general de las siguientes ecuaciones: 


5. B = 0. 

dx z dx 

„ „ d 3 y _ d 2 y dy 

6 . 2 —— + 5—— + 2— = 0. 

dx ó dx z dx 

7 «A 11A , n d*y L , d y 0 _ n 

■ 6 dx 4 U dx 3+2 dx 2+5 dx y 

o A 

8 ‘ d^~ V ~ ' 

O o A 7 A , 0 d 3 y d 2 y dy _ 

9 ' ^ + 2 díí “ + 7 Tx ~ 2V ~ °- 

10 fl-2^-+2^--i^- + ^--2v = 0 

dx 5 dx A dx 3 dx 2 dx J ' 


2.3 Ecuaciones no homogéneas con coeficientes 

constantes 


En esta sección, estudiamos dos métodos para hallar una integral particular de la 
ecuación no homogénea con coeficientes constantes. El primer método es bastante 
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fácil de aplicarlo pero no es aplicable a cualquier ecuación. El segundo método se 
puede aplicar a cualquier ecuación pero, en general, es más difícil que el primer 
método sobre todo porque este involucra la resolución de sistemas de ecuaciones e 
integración. 

2.3.1 Método de los coeficientes indeterminados 

En esta parte, estudiamos un método explícito para hallar una solución particular de 
la ecuación diferencial lineal con coeficientes constantes no homogénea 

dv d n v 

aoy + ai-f- + ... + a n — = F(x), a n ^ 0. (2.15) 

ax ax n 

El método, llamado método de los coeficientes indeterminados, desgraciadamente 
se puede aplicar solamente a una clase bastante restringida de funciones F(x) pero 
tiene la ventaja (cuando se puede aplicar) de ser relativamente sencillo. Antes de 
explicar el método necesitamos dar la siguiente definición. 

Definición 2.9. Se llama función de tipo CI a cualquier función de la forma: 

1. x n para todo n > 0. 

2. e ax donde a es una constante diferente de cero. 

3. sen(ax + ¡3) donde a y ¡3 son constantes y a ^ 0. 

4- eos (ax + (3) donde a y ¡3 son constantes y a 0. 

5. Cualquier función que se pueda conseguir como una suma finita de productos 
finitos de funciones de las clases anteriores. 

Las funciones x 3 , e 3x , sen(2a;), cos(4a;) son ejemplos de funciones de tipo CI. 
Además, la función 5x 3 e 3x + sen( 2 x) cos(4a;) — 4e 3x cos(4;r) también es una función 
de tipo CI puesto que se obtiene como una suma finita de productos finitos de las 
funciones de tipo CI anteriores. 

Observación 2.14. Aunque una función de tipo CI pertenece a una clase muy 
restringida de funciones, esta cubre una gama bastante amplia de funciones que 
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además tienen la ventaja de aparecer en muchas aplicaciones físicas. Resulta que el 
método de los coeficientes indeterminados se puede aplicar únicamente cuando la 
función F(x) es una función de tipo CI. 

La siguiente definición nos dice qué se entiende por conjunto asociado a una 
función de tipo CI. 

Definición 2.10. Sea f una función de tipo CI. Se llama conjunto CI asociado a la 
función f al conjunto CI/ formado por la función f y todas las funciones de tipo CI 
que se consiguen de f por derivación. 

Consideremos la función / definida como f(x) = x 3 . Se tiene que / es una función 
de tipo CI, luego tiene sentido hallar el conjunto CI asociado a x 3 . Al derivar se 
consiguen las funciones de tipo CI siguientes: x 2 ,x, 1. Por tanto, el conjunto CI T 3 es: 

CL^ = {a; 3 , a; 2 , x , 1}. 

Ejemplo 2.22. Hallar el conjunto CI asociado a la función f definida de la siguiente 
manera: f(x) = e 3x sena: + 7x 2 + cosa;. 

Solución. Primero, identifiquemos las funciones de tipo CI elementales (por decirlo 
de alguna manera) que forman parte de la función /. Es inmediato que las funciones 
e ¿x se ri a;, x 2 y cosa: son las funciones elementales de tipo CI que hacen parte de la 
función /. Derivando cada una de estas funciones se tiene: 

• e 3x sena; genera la nueva función (de tipo CI) e 3x cosx. 

• x 2 genera las nuevas funciones x y 1. 

• cosa: genera la nueva función sena:. 

Luego el conjunto CI asociado a la función / es el siguiente conjunto: 

CI/ = {e 3x sena:, e 3x cosa:, a; 2 , a;, 1, cosa;, sena:}. 
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El método CI 


Sea 


dy 

«o y + a i ~r~ 

dx 


+ ... + a n 


d n y 

dx n 


0 


la ecuación diferencial lineal homogénea asociada la ecuación diferencial 2.15. Sean 
fi, f 2 ,, f n las n soluciones linealmente independientes de la ecuación homogénea 
asociada (ver Sección 2.2). 


El método (algoritmo) de los coeficientes indeterminados consta de los siguientes 
pasos: 


Paso 1. Identificamos las funciones elementales de tipo CI que forman parte de la 
función F (ver ecuación 2.15). Sean gi,g 2 , ■ ■ ■ ,g m estas funciones. 

Paso 2. Calculamos los conjuntos de tipo CI asociados a cada una de las funciones 
gi para 1 < i < m. Sean CI 9l , CI 92 ..., CI 3m estos conjuntos. 

Paso 3. Si existen i,j, 1 < i, j < m tal que CI 9i C CI 9 ., entonces eliminamos el 
conjunto más pequeño, es decir, eliminamos el conjunto CI 9i . 


Sean CI 9l , CI S2 ,..., CI gfc , con k < m, los conjuntos que quedan después del paso 
3. 

Paso 4. Si fi G CI 9t para algún ¿ 6 { 1 , 2 ,...,?*} y algún t G { 1 , 2 ,..., k}, entonces 
multiplicamos cada elemento del conjunto CI St por la menor potencia positiva 
de x de tal forma que fi f Cl gt . 

Paso 5. Repetimos el paso 4 hasta que fi ^ CI gí para todo i G {1,2,..., ?z} y para 
todo í G {1,2,..., k}. 

Observamos que después del Paso 4 (posiblemente) algunos conjuntos CI fft habrán 
sido modificados. Con todos los conjuntos modificados y los que no necesitaron ser 
modificados pasamos al siguiente paso: 
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Paso 6. Ponemos y p , la integral particular (ver definición 2.7), como una combinación 
lineal de todos los elementos de todos los conjuntos que quedaron después 
del paso 5. 


Paso 7. Hallamos los valores de los coeficientes de la combinación lineal del paso 6. 

Observación 2.15. El algoritmo para hallar una integral particular y p de la ecuación 
2.15 parece demasiado artificioso y muy complicado de utilizar, pero (como cualquier 
procedimiento nuevo), después de practicarlo un par de ocasiones se va haciendo más 
natural y fácil de utilizar. 


Ejemplo 2.23. Hallar una integral particular de la ecuación 

d2 y o d V , o 2 x 
--3— + 2y = xe . 

dx ¿ d 


Solución. Primero hallamos las soluciones linealmente independientes de la ecuación 

d 2 y „ dy 

— - 3 - +2y = 0. 

Estas funciones solución están dadas por: f\ (x) = e x y / 2 (:r) = e 2x (comprobarlo). 
Ahora pasamos al algoritmo para determinar una integral particular y p . 


El paso 1 pide identificar las funciones elementales de tipo CI que hacen parte 
de la función definida por x 2 e x . En este ejemplo, la única función de tipo CI es la 
función x 2 e x . Esto termina el paso 1. 


El paso 2 pide calcular los conjuntos de tipo CI asociados a las funciones 
determinadas en el paso 1. Tenemos un sólo conjunto de tipo CI: 


CL 


= {x 2 


, xe 


C }- 


El paso 3 se omite en este ejemplo pues solo contamos con un conjunto. 

El paso 4 pide modificar el conjunto Cl x 2 e x pues la solución f 1 es elemento de 
este conjunto. Tenemos que multiplicar cada elemento del conjunto por la menor 
potencia de x de tal forma que /j ^ Cl x 2 e x. El conjunto que se obtiene de esta forma 
es el conjunto: 

CI x 2 e * = {x 3 e x , x 2 e x , xe x }. 
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El paso 5 no necesita ser aplicado, pues ni /j ni /2 son elementos del conjunto 
(en este caso modificado) Cl x 2 eX . 

El paso 6 nos dice que debemos poner la integral particular como una combinación 
lineal del único conjunto que tenemos (en este ejemplo). Así, y p esta dada por: 

y p = Ax 3 e x + Bx 2 e x + Cxe x , 


donde A, B, C son coeficientes indeterminadas (de ahí el nombre del método) que se 
pide calcular en el paso 7. 

Si realizamos los cálculos se tiene que 

y' p = ( Ax 3 + (3A + B)x 2 + (2 B + C)x + C)e x . 

y p = ( Ax 3 + (6 A + B)x 2 + (6 A + 4 B + C)x + 2 B + 2 C)e x . 


Al reemplazar los valores de y p , y' p y y " en la ecuación diferencial 

d 2 y dy 

—— — 3— + 2 y — x 2 e x obtenemos la siguiente igualdad: 
ax 2 a 

(-3Ax 2 + ( 6 A - 2 B)x + 2B- C)e x = x 2 e x . 


Para que esta igualdad se cumpla los coeficientes A, B y C deben satisfacer el 
siguiente sistema: 

r 

-3A = 1 
< 6A - 2B = 0 
2B - C = 0 

Resolviendo el sistema (hacerlo) se tiene A — —1/3, B = —1 y C = — 2 . Así la 
integral particular es: 

y p = -x 2 - 2x^j e x . 

Por el teorema 2.2 y la Definición 2.7, la solución general de la ecuación diferencial 
es y(x) = y c + y p donde y c es la función complementaria y y p es una integral particular 
(ver definición 2.7), esto es: 
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Ejemplo 2.24. Hallar una integral particular de la ecuación 


d^y 
dx 2 



+ 3 y 


2x 2 + e x + 2xe x + 4e 3x . 


Solución. Primero hallamos las soluciones linealmente independientes de la ecuación 
homogénea asociada. Es decir, hallamos el conjunto fundamental de soluciones de la 


ecuación 


d 2 V A dy 

dx 2 A dx + V 


0 . 


Aplicando el método estudiado en la sección 2 . 2 , vemos que el conjunto 
fundamental de soluciones de la ecuación homogénea está formado por las funciones 
definidas por: fi(x) = e x y /^(x) = e 3x (hacer las cuentas). Ahora pasamos al 
algoritmo para determinar una integral particular y p . 

El paso 1 pide identificar las funciones elementales de tipo CI que hacen parte 
de la función definida por 2x 2 + e x + 2xe x + 4e 3x . En este ejemplo, las funciones 
elementales de tipo CI son x 2 , e x , xe x y e 3x . Esto termina el paso 1. 

El paso 2 pide calcular los conjuntos de tipo CI asociados a las funciones 
determinadas en el paso 1. Tenemos cuatro conjuntos de tipo CI: 


CI X 2 = {x 2 ,x, 1}, 

Cíe* = K}, 
CI xe „ = {xe x ,e x }, 
CIg 3 x = {e 3x }. 


El paso 3 pide eliminar los conjuntos más pequeños. En este caso, eliminamos el 
conjunto CI e ^. 

El paso 4 pide modificar los conjuntos que tienen algún elemento coincidente 
con alguna de las soluciones fundamentales de la ecuación homogénea asociada. En 
nuestro caso tenemos que modificar los conjuntos CI xe * y Cl e 3x pues f\ G C4 xe x y 
/2 G CIf.3:':. Tenemos que multiplicar cada elemento de los conjuntos CI Ke .T: y CI e 3x 
por la menor potencia de x de tal forma que f\ ^ Cl xe x y f 2 ^ CI e 3 x. Los conjuntos 
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que se obtienen de esta forma son los conjuntos: 

CI xe * = {x 2 e x , xe x }, 
Cl e 3 * = {xe 3x }. 


El paso 5 no necesita ser aplicado pues ni / 4 ni f 2 son elementos de los conjuntos 
(algunos modificados en este caso) CI X 2 , CI xe x y CI e 3 i. 

El paso 6 nos dice que debemos poner la integral particular como una combinación 
lineal de los conjuntos que obtuvimos en el paso 5. Así, y p esta dada por: 

y p = A\x 2 + A 2 x + A 3 + A 4 x 2 e x + A 5 xe x + A 6 xe 3x , 


donde Ai, A 2 , A 3 , A 4 , A 5 y A 6 son coeficientes indeterminados que se pide calcular¬ 
en el paso 7. 

Si realizamos los cálculos se tiene que 


y' p = 2A ± x + A 2 + A 4 x 2 e x + (2 A 4 + A 5 )xe x + A 5 e x + 3 A 6 xe 3x + A 6 e 3x . 
y" = 2Ai + A 4 x 2 e x + (4 A 4 + A 5 )xe x + (2 A 4 + 2 A 5 )e x + 9 A 6 xe 3x + 6 A 6 e 3x . 


Al reemplazar los valores de y p , y' p y y' p en la ecuación diferencial 
y dy 

—— — 4— + 3y = 2x 2 + e x + 2xe x + 4e 3x obtenemos el sigue nte sistema de 
dx ¿ dx 

ecuaciones: 

3A] = 2 

— 8 Ai + 3A2 = 0 
2Ai - 4 A 2 + 3A 3 = 0 
-4A 4 = 2 
2A 4 - 2A 5 = 1 
2A 6 = 4 

Resolviendo el sistema (hacerlo) se tiene Ai = 2/3, A 2 = 16/9, A 3 = 52/27, 
A 4 = — 1 / 2 , A 5 = —1 y A e = 2 . Así la integral particular es: 

2 2 16 52 A 2 \ ,, 0 sj, 

Vp= 3 X + Y X + 27 ~ \ 2 X +x ) e + 2xe ■ 
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Por el teorema 2.2 y la definición 2.7, la solución general de la ecuación diferencial es 
y(x) — y,, + y P donde y c es la función complementaria y y p es una integral particular. 
Así: 

y(x ) = gU + —x + — - í -x +x-c 1 je + (2x + C 2 )e. 

Ej ercicios 


Hallar la solución general de las siguientes ecuaciones: 


1. 3 ^ + 3 ^ -f -^zsenz + teV. 
dx ó dx z dx 

d A y d 3 y d 2 y dy 2x 3 

2 -i~~ 6 y = xe ~ 3x • 

dx dx ó dx z dx 

3 . pL + pL-7pL- d JL + 6y = e * + xe -* + 2e -* + 5 . 

dx 4 dx ó dx z dx 


4. -12^ + 35^-59^| + 327^|+441^ + 108?/ = sen(3x) + 2x cos(3a:) 
dx ü dx 4 dx ó dx z dx 


xe 


. 4 . 7 : 


d 6 y d 4 y d 2 y 
dx 6 dx A dx 2 


+ y — sen x + 2 eos x — 3xe x + 4x e 


2 —x 


R 0 d 2 V K dy 2 

o. 2—- — 5-—h y = x +x — 1. 
dx- dx 


2.3.2 Método de variación de parámetros 


En la subsección 2.3.1, estudiamos el método de coeficientes indeterminados, el cual 
involucra únicamente derivación y técnicas de resolución de sistemas de ecuaciones 
lineales. Se indicó también que el método CI se aplica únicamente al caso de funciones 
de tipo CI. Por ejemplo, el método CI no se aplica a la siguiente ecuación: 

d 2 y 

+ y — tan x. 
dx z 

En esta sección estudiamos un método más general para hallar una integral particular 
de una ecuación diferencial de n-ésimo orden que además, tiene la ventaja de aplicarse 
a ecuaciones lineales con coeficientes variables. 
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Primero explicamos el método de variación de parámentros al caso particular 
de una ecuación de segundo orden y posteriormente lo generalizamos al caso de 
ecuaciones de orden n > 3. 


Consideremos la ecuación diferencial lineal de segundo orden 
a a (x)y + Oi(i)^ + a 2 (x)^¡ = F(x). 


(2.16) 


Sean y 1 , 1/2 dos soluciones linelmente independientes de la correspondiente ecuación 
homogénea asociada 

a 0 (x)y + a Á x )^ + = O’ ( 2 ‘ 17 ) 

La función complementaria de la ecuación 2.16 está dada por 

y c (x) = ciyi(x) + c 2 y 2 (x) 
donde ci,C 2 son constantes arbitrarias. 

El método de variación de parámetros asegura que una integral particular de la 
ecuación 2.16 tiene la forma 

y p (x) = d(x)yi(x) + c 2 (x)y 2 (x). 

Observación 2.16. Vemos que el método de variación de parámetros supone que 
una integral particular se puede conseguir haciendo variar los parámetros (de ahí el 
nombre) c\ y c 2 de la función complementaria. 

Para hallar las funciones c\(x) y c 2 (x) procedemos de la siguiente manera: 
Derivamos la función y p : 

y' p (x ) = yi{x)c\(x) + ci(x)r/í(a;) + y 2 (x)c' 2 (x) + c 2 (x)y' 2 (x). (2.18) 


Ahora exigimos que 


yi( x )c 1 (x) + y 2 (x)c 2 (x) = 0. 


(2.19) 


( 2 . 20 ) 


Con esta suposición, la ecuación 2.18 toma la forma 

y' P {x ) = dix^ix) + c 2 (x)y 2 (x). 
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Derivamos la ecuación 2.20: 

y''(x) = c' 1 (x)y[(x) + a(c)y"(x) + c 2 (x)y 2 (x) + c 2 (x)y 2 (x). 


Reemplazamos estos valores en la ecuación 2.16: 

a 0 (x)[c 1 (x)y 1 (x) + c 2 (x)y 2 (x)\ + a 1 (x)[c 1 (x)y[(x) + c 2 (x)y 2 (x)\ + 
a 2 (x)[c' 1 (x)y' 1 (x) + a(c)y"(x) + c 2 (x)y' 2 (x) + c 2 (x)y 2 (x)\ = F(x) 

Reagrupando términos tenemos: 

ci(x)[ao(;x)yi{x) + ai(x)y[(x) + a 2 (x)y"(x)} + 
c 2 (x)[a 0 (x)y 2 (x) + a x {;x)y' 2 (x) + a 2 (x)y 2 (x)\ + 
a 2 (x)[y[(x)c 'i(x) +y 2 (x)c 2 (x)} = F(x) 


Ya que yi, y 2 son soluciones de la ecuación 2.17, la última ecuación se reduce a 
la ecuación: 

Vi^y^x) + y 2 (x)c 2 {x) = ( 2 . 21 ) 

a 2 (x) 

Las ecuaciones 2.19 y 2.21 forma el sistema de ecuaciones 


2 /i(a;)ci(a;) + y 2 (x)c' 2 (x) = 0 


y'^xy^x) + y' 2 (x) c ' 2 {x) = 


F(x) 

a 2 (x) 


( 2 . 22 ) 


donde las incógnitas son las funciones c\ y c' 2 . 

Observación 2.17. Puesto que a 2 (x) ^ 0, la ecuación 2.21 está bien planteada. 
Luego, tiene sentido preguntarnos sobre la solución del sistema 2.22. 


Puesto que el wronskiano W(yi,y 2 ) de yi,y 2 es diferente de cero (ver definición 
2.5 y teoremas 2.4), este sistema tiene solución única y esta viene dada por: 


0 

V2 (x) 


F(x) 

a 2 (x) 

y' 2 (y 

y 2 (x)F(x) 

yi(x) 

y 2 (x) 

a 2 (x)W(y 1 ,y 2 )' 

y'iix) 

y 2 (x) 



94 











Leónidas Cerda y Janneth Morocho 


Vi(x) 

0 


y'A x ) 

F{x) 

a 2 (x) 

yi(x)F(x) 

yi( x ) 

V2(x) 

a 2 ( x )W(y l ,y 2 ) 

y[( x ) 

y 2 ( x ) 



Así, las funciones c\ y C 2 están dadas por: 


c\{x) 


y 2 (x)F(x) 
a 2 (x)W(y 1 ,y 2 ) 


dx, 


C 2 (x) 


y 1 {x)F{x) 

a 2 (x)W(yi,y 2 ) 


dx. 


Ejemplo 2.25. Hallar una integral particular de la siguiente ecuación: 

d 2 y 

+ y = tan x. 
dx 


Solución. Primero hallamos las soluciones linealmente independientes de la ecuación 
homogénea asociada. Realizando los cálculos (hecerlo) se tiene que y\{x) = cosa;, 
y 2 (x) = seno; son las soluciones buscadas. Ahora pedimos que una integral particular 
de la ecuación no homogénea tenga la forma y p (x) = Ci(x) eos x + c 2 (x) seno;. A partir 
de este punto podríamos repetir los pasos que se hizo para hallar los valores de las 
funciones Ci y C 2 sin embargo, vamos a utilizar el último resultado de los cálculos 
anteriores pues este nos da una fórmula explícita para las funciones C\ y c 2 . 


El wronskiano de las funciones seno;, coso; es: 


IP^sen x, eos x) = 


sen x eos x 
eos x — sen x 


= - 1 . 


Por lo tanto: 


= í i h(x)F(x) dx = _ 
J a 2 (x)W(y 1 ,y 2 ) 


sen 2 X 


dx 


cosx 

(seco; — eos x)dx 


= — ln(sec x + tan x) + sen x. 

, \ f Vi{x)F(x) f 

c 2 (x) = / — . - rdx = — sen xdx 

J a 2 (x)W ( 2/1 ,2/2) J 

= coso;. 

De esta forma tenemos que una integral particular es: 


y p (x) = — coso; ln(seco; + tana;) + sen(2a;). 
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Observación 2.18. La solución general de la ecuación diferencial esta dada por la 
función definida como 

y(x) = ci cosa: + C 2 sena: — eos x ln(sec x + tana:) + sen(2x), 
donde Ci,C 2 son constantes arbitrarias. 

Ahora vemos el método de variación de parámetros para el caso de una ecuación 
diferencial de n-ésimo orden. 

Consideremos la ecuación diferencial lineal de n-ésimo orden 



(2.23) 


y su correspondiente ecuación homogénea asociada 



(2.24) 


Sean a/i, 2 / 2 , • ■ ■ ,y n soluciones linealmente independientes de la ecuación 2.24. La 
función complementaria de la ecuación no homogénea 2.23 está dada por 


y c (x) = ciyi(x) + c 2 y 2 (x) + ... + c n y n (x ) 


donde ci,c 2 , ■ ■ ■ ,c n son constantes arbitrarias. 

Según el método de variación de parámetros, una integral particular tiene la 
forma 


y p (x) = ci(x)yi(x) + c 2 (x)y 2 (x) + ... + c n (x)y n (x). 


Derivamos la función y p \ 


Vp = yi( x ) c 'i( x ) + ci(x)y[(x) + ... + y n (x)c' n (x) + c n (x)y' n (x). (2.25) 


Pedimos que 


2/i ( x ) c i( x ) + ■ • • + y n {x)c' n (x) = 0 


Luego, la ecuación 2.25 se simplifica a 


y' P ( x ) = ci(x)y[(x) + ... + c n (x)y' n (x) = 0. 


(2.26) 
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Derivamos y' p de la ecuación 2.26 para obtener la segunda derivada de la integral 
particular: 

y'p O) = y'i^y^x) + c^y'^x) + ... + y' n (x)c' n (x) + c n (x)y”(x). 

Ahora pedimos que 

y'Áxy^x) + ... + y' n (x)d n (x) = 0 


Repitiendo el mismo proceso hasta la derivada de orden n de la función y p , 
exigiendo en cada derivada que las partes que contienen la derivada de Ci(x) para 
cada 1 < i < n sean igual a cero y pidiendo que la última derivada (y las derivadas 
de orden menor a n) satisfaga la ecuación 2.23, obtenemos el siguiente sistema de 
ecuaciones: 


y 1 (xy^x) + y 2 {x)d 2 (x) + ... + y n (x)d n (x ) = 0 
Z/íO) c í(X) + y' 2 ( x )c' 2 (x ) + ... + y' n (x)d n (x) = 0 

< 


y[ n X) (aj)cíi(a:) + y { 2 1] (x)d 2 (x) + ... + l \x)d n (x) 


F{x) 

a n (x)' 


(2.27) 


Si ponemos W¿(/i,..., f n ) al determinante que se obtiene del wronskiano 

/ F(x)\ T 

de /i,... , f n reemplazando la i-ésima columna por 0, 0 ,..., ——— ) , donde el 

V a n {x)J 

superíndice T indica traspuesta, entonces la i-ésima solución d^x) está dada por: 

,/ \ Wi(j/i,...,j/ n ) 

CÁX) = W-tol.»,)■ 

Observación 2.19. La solución del sistema la estamos calculando por la regla de 
Cramer; sin embargo, en muchos casos en más rápido buscar otra forma de hallar 
las expresiones para c((x). 


Para terminar, vemos que las funciones c¿ están dadas por la siguiente expresión: 


Ci (x) 


Wi(y u .. ., 2 / n ) 

W(y u ...,y n ) dX 


El siguiente ejemplo muestra la forma de aplicar el método de variación de 
parámetros a una ecuación de segundo orden con coeficientes variables. 
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Ejemplo 2.26. Hallar una integral particular de la ecuación 

2 d2 y R d y o 4 , R 3 

x — ox -— h 10 y = 3x + bx , 

dx ¿ dx 

dado que y\ = x 2 , y 2 = x 5 son soluciones linealmente independientes de la 
correspondiente ecuación homogénea asociada. 


Solución. Se tiene que una integral particular y p está dada por 

y p (x) = x 2 ci{x) + x 3 c 2 (x). 


Para hallar los valores de Ci(x) para i = 1, 2 tenemos que resolver si sistema 


{ x 2 d 1 {x) + x 5 c' 2 (x ) = 0 
2xc\ (x) + 5x 4 c' 2 (x ) = 3x 2 + 6x. 
Utilizando la regla de Cramer, se tiene 

0 x 5 
3x 2 + 6x [ 

C 'Á X ) = - 

X X o 

2x 5x 4 


x 2 0 

2x 3x 2 + 6íc 

3x 4 + 6a; 3 3a; 4 + 6a; 3 


2 5 


5a: 6 - 2a; 6 3a: 6 


X X 




2x 5x 4 




3x 7 + 6x 6 3x 7 + 6x d 

5x 6 - 2x e ~ 3^6 


Integrando, se tiene que Ci(x) 

V P {x) ~- 


-- -x 2 + 2x, c 2 (x) = - 

^ x 4 + 2x 3 — x 4 — x 3 
— ^x 4 + x 3 . 


1 

X 


—. Por lo tanto 
x ¿ 


Observación 2.20. Notamos que el método de variación de parámetros presenta un 
gran problema. Efectivamente, para poder aplicar el método, necesitamos conocer el 
conjunto fundamental de soluciones de la ecuación homogénea asociada. En general, 
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para una ecuación con coeficientes variables, esto es bastante difícil de averiguar. Sin 


embargo, para el caso de ecuaciones con coeficientes constantes averiguar cuál es este 


conjunto fundamental de soluciones no representa ningún problema (en la mayoría 
de casos). 


(fiy 

Ejemplo 2.27. Hallar una integral particular de la ecuación 


y = tana;. 


Solución. La ecuación homogénea asociada a esta ecuación es —- — y = 0, su 

dx 2 

ecuación característica asociada m 2 — 1 = 0 tiene raíces mi = 1, m 2 — — 1. Luego, 
el conjunto fundamental de soluciones de la ecuación homogénea está formado por 
las funciones definidas como yfix) = e x , y- 2 (x) = e~ x . Por el método de variación de 
parámetros, una integral particular tiene la forma y p (x) = e x C\{x ) + e~ x c 2 (x). 


Las derivadas de las funciones c\ , c 2 vienen dadas por 


0 e~ x 
tana; — e~ x 

sen x eos x 

eos x — sen x 


e x tan x 
— sen 2 x — eos 2 x 


—e x tana;, 



e x 

0 




e x 

tana; 


e x tan a; 

sena; 

cosa; 


— sen 2 x — eos 2 x 

cosa; 

— sena; 



—e x tan x. 


Observación 2.21. Este ejemplo muestra que el último paso en la aplicación del 
método de variación de parámetros puede resultar bastante complicado. En efecto, 
las integrales que se generan con este método pueden ser muy difíciles de calcular. 
Invitamos al lector intentar hallar los valores de cifx), c 2 (x) del ejemplo anterior. 


Ej ercicios 


Hallar la solución general de las siguientes ecuaciones: 


i o d3y o d2y 

' dx 3 dx 2 


dy 

— - y = sena;. 

dx 
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4. 


2 . 


3. 




2.4 Ecuaciones lineales homogéneas con 


coeficientes variables 


En la sección 2.2 vimos que la solución de una ecuación diferencial lineal de n-ésimo 
orden homogénea con coeficientes constantes se pueden expresar como la combinación 
lineal finita de funciones elementales. Para las ecuaciones diferenciales homogéneas 
de orden superior con coeficientes variables sus soluciones, en general, no se pueden 
expresar de manera tan sencilla. 

En esta sección, vamos a ver que alguna solución de la ecuación 



(2.28) 


donde ü 2 { x ) 0, se puede expresar como una serie de potencias. Primero tenemos 

que dar algunas definiciones y enunciar algunos resultados preliminares. 

Observación 2.22. Se llama serie de potencias de x — x o a cualquier suma de la 


OO 



donde los c n son coeficientes arbitrarios. 


Ahora recordamos el concepto de serie de Taylor de una función / centrada en 
un punto xq. 


Definición 2.11. Sea f una función que admite derivada de cualquier orden en xq. 
Se llama serie de Taylor de f centrada en Xq a la siguiente serie de potencias de 


x — Xq: 



n =0 
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Por ejemplo, la serie de Taylor de la función e x , centrada en xq = 0 es 



de la función sen ay centrada en xq = 0 es 



de la función eos x centrada en = 0 es 



1 

720 


x 6 + ... 


La siguiente definición clasifica una clase importante de funciones. 

Definición 2.12. Una función f se dice analítica en xo si su serie de Taylor centrada 
en Xq existe y es convergente a f(x) para todo x en algún intervalo que contiene a x 0 . 

Notamos que la función exponencial e x , así como las funciones sen x, cosa: y los 
polinomios a o + a±x + ... + a n x n son funciones analíticas en todos los reales. Además, 

p( x ) 

las funciones racionales son analíticas en todos los puntos x tales que Q(x) ^ 0. 

Q\ x ) 

Por ejemplo, la función l/(a; 2 — 1) es analítica en todo x ^ ±1. 

La siguiente definición es importante en nuestro proceso de hallar las soluciones 
de una ecuación diferencial de segundo orden con coeficientes variables. 

Definición 2.13. Se dice que una ecuación diferencial lineal homogénea de segundo 
orden esta expresada en forma normal si tiene la forma (o se puede escribir de la 


forma) 



(2.29) 


Puesto que 02 ( 2 ;) 7 ^ 0 (ver ecuación 2.28), podemos hablar de normalizar una 
ecuación homogénea de segundo orden. En efecto, lo único que tenemos que hacer es 
dividir la ecuación 2.28 por a 2 (x). 

Ejemplo 2.28. Normalizar la ecuación 
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Solución. Dividiendo la ecuación por x 2 obtenemos la forma normal de la ecuación. 


2 d 2 y r d V 


Es decir, la forma normal de la ecuación x ¿ —- — 5 ——\- (x — l)y = 0 es la ecuación 

dx 2 dx 

d 2 y 5 dy x - 1 _ 
dx 2 x 2 dx x 2 ^ 


Observación 2.23. El proceso que realizamos para hallar la forma normal de la 
ecuación diferencial del ejemplo 2.28 tiene sentido siempre que ü 2 (x) ^ 0. Luego 
tenemos que tener mucho cuidado al normalizar una ecuación. 


Definición 2.14. Se dice que un punto Xq es un punto ordinario de la ecuación 2.28 
si las funciones q± y q 2 de la ecuación 2.29 son ambas funciones analíticas en xq. 


Si el punto Xq no es un punto ordinario, diremos que es un punto singular. 


Ejemplo 2.29. Hallar los puntos ordinarios de la siguiente ecuación: 


[x 


d 2 y dy 2 

i)— - 5x -—b x y 
dx 2 dx 


0 . 


Solución. Normalizamos la ecuación: 


d 2 y 

dx 2 


5x dy x 2 

~3 TT ' 3 T y = 

x ó — 1 dx x ó — 1 


Ya que x 3 — 1 = 0 si y solo si x — 1 (las dos raíces complejas del polinomio x 3 — 1 
no las tomamos en cuenta puesto que los únicos puntos ordinarios que nos interesan 
son reales) tenemos que los puntos ordinarios de la ecuación son todos los (puntos) 
reales excepto el punto x — 1. Así, el punto x = 1 es el único punto singular de la 
ecuación. 


2.4.1 Solución alrededor de puntos ordinarios 


El siguiente teorema da una condición suficiente para la existencia de soluciones en 
series de potencias de una ecuación diferencial. 

Teorema 2.8. Si x o es un punto ordinario de la ecuación 2.28, entonces esta ecuación 
tiene dos soluciones linealmente independientes de la forma 

OO 

^c n (x - x 0 ) n . 

n=0 
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Observamos que el teorema 2.8 no solo asegura la existencia de soluciones 
linealmente independientes, sino que también nos asegura que estas series de potencias 
convergen en algún intervalo |x — x 0 | < R (donde R > 0) alrededor de x 0 . 

El siguiente ejemplo muestra la forma de hallar dos soluciones linealmente 
independientes, expresadas en series de potencias alrededor de un punto ordinario de 
una ecuación diferencial. 


Ejemplo 2.30. Hallar dos soluciones linealmente independientes de la ecuación 

(2.30) 


d2 y , .... d V , / 2 


dr- +x di + ix +2)v = 0 ' 


Solución. En este caso, se tiene que todo punto xq es un punto ordinario de la 
ecuación 2.30. Para facilitar los cálculos elegimos xq = 0. 

Así, suponemos que la solución que estamos buscando tiene la forma 


y = J2‘ 


= > C n x 
n =0 


(2.31) 


Hallamos la primera y segunda derivada de la ecuación 2.31: 

dy 


dx 

d 2 y 

dx 2 


E 


= > nc n x 


71—1 


n= 1 
oo 


íl = Y / n (n-l)cn x 


71—2 


77—2 


Reemplazando estos valores en la ecuación diferencial 2.30 obtenemos 

OO OO OO OO 

^ n{n — 1 )c n x n ~ 2 + ^ nc n x n + ^ c n x n+ 2 + 2 ^ c n x n = 0. (2.32) 


71=2 


77=1 


77=0 


77=0 


Para que las soluciones tengan la forma adecuada necesitamos que las potencias 
de x sean todas iguales a n. Luego la primera y tercera suma tienen que reescribirse. 

OO OO 

y^n(n - 1 )c n x n ~ 2 = ^(n + 2)(n + 1 )c n+2 x n , 

77=0 
OO 

£■ 


77=2 


^ ^ c n x y c n — 2 X . 


77=0 


77=2 
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Reemplazando estas últimas expresiones en la ecuación 2.32 obtenemos 

OO OO OO OO 

y^(n + 2)(n + 1 )c n+2 x n + ^ nc n x n + ^ c n - 2 x n + 2 ^ c n x n = 0. (2.33) 

n =0 n= 1 n =2 72=0 


Para realizar operaciones algebraicas necesitamos que todas las sumas de la 
ecuación 2.33 comiencen desde el mismo valor. Para lograr esto soltamos algunos 
términos de las sumas que comienzan con índice de sumación más bajo. Así, tenemos 


que: 

OO OO 

y^(n + 2)(n + 1 )c n+2 x n = 2 c 2 + 6c 3 x + ^(n + 2)(n + 1 )c n+2 x n , 

n =0 n =2 

OO OO 

nc n x n = Cjx + nc n x n y 

n=l 72=2 

OO OO 

c íl a: íl = c 0 + cqx + c n x n . 

72=0 71=2 


Reemplazando estas igualdades en la ecuación 2.33 obtenemos 

OO OO OO 

2 c 2 + 6c 3 x + ^(n + 2 )(n + l)c n+2 x n + c±x + ^ nc n x n + ^ c n - 2 x n 

n =2 72—2 72—2 

OO 

+2c 0 + 2c 3 x + 2 c n x n = 0. 

72—2 


Agrupando términos, la última ecuación se puede escribir como: 


(2c 0 + 2c 2 ) + (3ci + 6c 3 )x + ^ ^ [(n + 2)(n + l)c n _|_ 2 + (n + 2)c n + c n _ 2 ]x n — 0. (2.34) 

72—2 

Para que la ecuación 2.34 se satisfaga se deben cumplir las siguientes condiciones: 


Condición 1. 2 cq + 2c 2 = 0, 


Condición 2. 3cq + 6c 3 = 0 y, 


Condición 3. (n + 2)(n + l)c n+2 + (n + 2)c n + c n _ 2 = 0 para todo n > 2. 


Desde la condición 1, tenemos 


c 2 — — Cq. (2.35) 


A partir de la condición 2, tenemos 

c 3 = --ci. (2.36) 
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Finalmente, desde la condición 3, tenemos 

(n + 2 )c n + c n _ 2 , . 

Cn+2 = — 7 — -7777 — —7^, valida para n > 2. (2.37) 

(n + 2 )(n + 1 ) 

La ecuación 2.37 nos proporcionan una fórmula recursiva para los valores de la 
constante c íl+2 . Algunos valores son: 


n 

£n +2 

2 

4c 2 + c 0 1 

C 4 = - —— = -Co, ver ecuación 2.35 

3 

5^3 “1“ 3 

° 5 = -ññ— = — Cj., ver ecuación 2.36 

4 

c 6 =- ~ 4 2 =-c 0 , ver el valor de c 4 y ecuación 2.35 

b 30 60 ’ 

5 

C7 = — —LÍ—4 —- c ver valor de c 5 y ecuación 2.36 

7 42 1680 J 


Observación 2.24. Es inmediato ver que todos los coeficientes c n con n par están 
expresados en términos del coeficiente co, mientras que los coeficientes c n con n impar 
están expresados en términos del coeficiente c\. 


Reemplazando estos valores en la ecuación 2.31 tenemos: 


y = C 0 + CiX — CnX 2 - CiX 3 H- CnX 4 H- C il 5 - CnX 6 - CiX 7 + 

y ° ° 2 4 ° 40 60 1680 


Agrupando términos, tenemos 

y = c 0 ^1 - x 2 + ^ x 4 + .. .^ + ci - ^x 3 + ^x 5 + .. .^ (2.38) 

Las series, entre paréntesis, de la ecuación 2.38 son la expanción en series de 

potencias de dos soluciones linealmente independientes de la ecuación 2.30 donde las 

constantes Cq y C\ son constantes arbitrarias; luego, la ecuación 2.38 es la solución 

dy 

general de la ecuación diferencial —— + x— + (x 2 + 2)y = 0 . 

dx ¿ dx 
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Ejemplo 2.31. Hallar la solución general de la ecuación 


2 <Py 

dx 2 


dy 

+ x-f- + y 
dx 


0. 


(2.39) 


o 1 . . T r ! , ! • . o on d 2 y , x dy 1 

Solución. La forma normal de la ecuación 2.39 es —- H-1 —y 

dx 2 2 dx 2 y 

punto Xq es un punto ordinario de la ecuación 2.39, elegimos Xq 


0. Luego, todo 
0 para que se 


faciliten los cálculos. 


Ahora procedemos igual que el ejemplo 2.30 para hallar dos soluciones linealmente 
independientes. 

Suponemos que la solución que estamos buscando tiene la forma 

OO 

y = ^2 c nX n . (2.40) 

72=0 


Hallamos la primera y segunda derivada de la ecuación 2.40: 

dy 


dx 

d 2 y 

dx 2 


= y nc ” x 


n— 1 


n= 1 
oo 


d 2 y v-^ , 1 ^ 

= ^ n ( n -l) CnX 


n —2 


n =2 


Reemplazando estos valores en la ecuación 2.39 obtenemos 

oo i OO 1 71 

y, n{n - 1 )c n x n ~ 2 + - y nc n x n + - y c n x n = 0. 

n =2 72—1 72—0 


(2.41) 


Para que las soluciones tengan la forma adecuada necesitamos que las potencias 
de x sean todas iguales a n. Así, la primera suma tienen que reescribirse. 


y n(n - 1) 


c n x n 2 = 


y(n + 2)(n + 1) 


C n+ 2 X 


72=2 


72=0 


Reemplazando esta última expresión en la ecuación 2.41, obtenemos 

X ^ OO 1 n 

y (n + 2 )(n + 1 )c n+2 x n + ~y nc n x n + ~y c n x n = 0. (2.42) 


72=0 


72=1 


72=0 


Necesitamos que todas las sumas de la ecuación 2.42 comiencen desde el mismo 
valor. Para lograr esto saltamos algunos términos de las sumas que comienzan con 
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índice de sumación más bajo. Así, tenemos que 

OO OO 

y^(n + 2 )(n + 1 )c n+2 x n = 2 c 2 + ^(ra + 2)(n + l)c„ + 2 x" 


72=0 

n 


n =1 


= Co + C n X n . 


72=0 ? 2=1 

Reemplazando estas igualdades en la ecuación 2.42, obtenemos 


oo ^ oo ii n 

2c 2 + y^(n + 2)[n + 1 )c n+2 x n + - ^ nc n x" + 2 C ° + 2 Cn2;íl = 

72—1 72 — 1 72—1 


Agrupando términos, la última ecuación se puede escribir como: 


1 , ;V 

2c 2 + -c 0 + 

” 72—1 


. .. 11 

(n + 2)(n + l)c n+2 + -nc n + -c n 


x n = 0. 


(2.43) 


Para que la ecuación 2.43 se satisfaga, se deben cumplir las siguientes condiciones: 


Condición 1. 2c 2 + -c 0 = 0 y, 


Condición 2. (n + 2 ){n + l)c n+2 + -nc n + -c n para todo n > 1. 


Desde la condición 1, tenemos 



(2.44) 


A partir de la condición 2, tenemos 


^ 72+2 


(\ U +\) C n 

(n + 2)(n + 1) ’ 


válida para n > 1. 


(2.45) 


La ecuación 2.45 nos proporcionan una fórmula recursiva para los valores de la 
constante c íl+2 . Algunos es estos valores son: 
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n 

^n+2 

1 

i 

C 3 = yCi 

6 

2 

C4 = ~¡2 C ° 

3 

Cí = ¿ Cl 

4 

1 

Cq — Cn 

b 384 ° 


Observación 2.25. Al igual que el ejemplo 2.30, todos los coeficientes c n con n par 
están expresados en términos del coeficiente Cq, mientras que los coeficientes c n con 
n impar están expresados en términos del coeficiente c\. 


Reemplazando estos valores en la ecuación 2.40, tenemos: 


1 2 1 3 

y = co + CíX - -c 0 x + -cix 
4 o 


32, C ° X + 



1 

384 


c 0 x 6 + ... 


Agrupando términos, tenemos 

y = c 0 1 1 — -x 2 ——x 4 -— x 6 — ) + ci f x + -x 3 + —x 5 + ... ) (2.46) 

y V 4 32 384 ) \ 6 60 ) 

Las series, entre paréntesis, de la ecuación 2.46 son la expansión en series de 

potencias de dos soluciones linealmente independientes de la ecuación 2.39 donde las 

constantes Co y C\ son constantes arbitrarias; luego, la ecuación 2.46 es la solución 

d 2 y dy 

general de la ecuación diferencial 2—- + x— + y = 0. 

dx ¿ dx 


Ej ercicios 


El siguiente bloque de ejercicios se ha tomado de Ross [RS]. 

Resolver las siguientes ecuaciones diferenciales alrededor de algún punto ordinario: 
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t , 2 i\ d2 y , o d v , n 

1. (x - 1)—r + 3x— + xy = 0. 

ax z dx 

d 2 ydy 2 

2. —j+£— + (£ — A)y = 0. 
dx z dx 

3. (x + 3)^-| + (x + 2)^ + y = 0. 

ax 2 aa; 


4. 


d 2 2/ dy 

~ x i—y = °> 

dx z dx 

1/(0) = 1, 

2/(0) = 0. 


r d 2 y dy 

~Á+ X ~r ~ 2y = 0 ’ 

dx z dx 


5. 


2 /( 0 ) = 0 , 
2 /( 0 ) = 1 . 


6 . 


(o 2 o \d 2y o dy , n 
< 21 - 3) d^~ 2x ^ + v = 0 ’ 

2 /( 0 ) = - 1 , 

2/'(0) = 5. 


2.4.2 Solución alrededor de puntos singulares 


El método empleado en la subsección 2.4.1 para hallar soluciones linealmente 
independientes de la forma 

OO 

2/ = y^c n (x - x 0 ) n 

n =0 

para la ecuación diferencial 


a 0 (x)y + ai(x) 


dy_ 

dx 


+ a 2 (x) 


d 2 y 

dx 2 


0 


no funciona cuando el punto Xo es un punto singular de la ecuación diferencial. En 
efecto, el teorema 2.8 solo es aplicable para el caso de puntos ordinarios. 


Cuando Xq es un punto singular, suponemos que las soluciones tienen la forma 

OO 

y = \x-x Q \ r ^c n (x - x 0 ) n , 

n =0 
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donde r es una constante (real o compleja). Antes de explicar el método, llamado 
método de Frobenius , necesitamos dar la siguiente definición. 


Definición 2.15. Sea xo un punto singular de la ecuación 


a 0 (x)y + ai (a;) 


dy 

dx 


+ a 2 (x) 


d 2 y 

dx 2 


0 . 


Se dice que xq es un punto singular regular si las funciones definidas como 


(x - x 0 ) 


ai(x) 

a 2 (x) 


y (x -x 0 ) 2 


ao{x) 

a 2 (x) 


son analíticas en xq. 


Observación 2.26. Utilizando la terminología de la definición 2.13 podríamos decir: 
Xo es un punto singular regular de la ecuación diferencial 2.28 (ver pag. 100) cuando 
las funciones definidas como (x — x 0 )q\{x) y (x — x 0 ) 2 q 2 (x) son analíticas en x 0 . 
Además, diremos que Xq es un punto singular irregular de la ecuación diferencial 
2.28 cuando Xq no es un punto singular regular. 


Ejemplo 2.32. Hallar los puntos singulares regulares (si existen) de la siguiente 
ecuación: 


(x 


^d 2 y dy 

l)q-^ + (7x + 3)— + 4(x 
dx ¿ dx 


3 )y = 0. 


Solución. La forma normal de la ecuación es 

d 2 y 7x + 3 dy 4(g - 3) = Q 

dx 2 x — 1 dx x — 1 


El punto Xo = 1 es un punto singular de la ecuación diferencial. Para averiguar si este 
punto singular es regular calculamos los productos (x — Xo)<?i(x) y (x — Xo) 2 q 2 (x). 

(x — xo)qi(x) — (x — 1) ——— 

x — 1 

= 7x + 3. 

(x - x 0 ) 2 g 2 (x) = (x - l) 2 ^—^ 

x — 1 

= 4x 2 — 16x + 12. 


Como las funciones definidas por (x — xo)<?i(x) y (x — xo ) 2 q 2 (x) son analíticas en 
x 0 = 1 se tiene que x 0 = 1 es un punto singular regular de la ecuación diferencial. 
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El siguiente teorema nos da una condición necesaria para que la ecuación 
diferencial 2.28 tenga soluciones en series de potencias alrededor de un punto singular. 

Teorema 2.9. Si Xq es un punto singular regular de la ecuación 2.28, entonces la 
ecuación 2.28 tiene al menos una solución de la forma 

OO 

I x ~ x 0 r '^2 l c n {;x - x 0 ) n , (2.47) 

71=0 

donde r es una constante. Además esta solución es valida para 0 <| x — xq \< R. 

Observación 2.27. La constante r puede ser real o compleja y el proceso que se 
sigue para hallar el valor de los coeficientes c n es, más o menos, el mismo que se 
realizó para el caso de puntos ordinarios. Este proceso se conoce comúnmente como 
método de Frobenius. Además, el método de Frobenius supone que Cq ^ 0 para hallar 
el valor de la constante r. 

Ejemplo 2.33. Hallar una solución, alrededor de un punto singular regular, de la 
ecuación diferencial 

2 d 2 y dy 

2x —r — x- —I- (x — 5 )y = (J. 
dx 2 dx 

Solución. Se puede demostrar que xq = 0 es un punto singular regular de la ecuación 
diferencial (hacerlo). Por el teorema 2.9, existe al menos una solución de la forma 

OO 

V =| x r 

n=0 

Observación 2.28. El método que vamos a explicar halla una solución válida para 
todo x tal que 0 < x < R, si queremos hallar soluciones válidas para —R < x < 0 
simplemente ponemos —x > 0. 


Ya que \ x \ — x tenemos 


V = ^2 c n x 

71=0 


7i+r 


donde c 0 ^ 0. 

Calculamos la primera y segunda derivada de y: 


dy_ 

dx 

d 2 y 

dx 2 


J^(n + r) 


c n x 


71 + 7 *— 1 


71=0 

OO 


y~~^(n + r — 1 )(n + r)c n x 


7i+r—2 


71=0 


(2.48) 
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Sustituyendo en la ecuación diferencial tenemos: 

OO OO OO OO 

2 + r — l)(n + r)c n x n+r — y~](n + r)c n a; n+? ’ + ^ c n a; n+í ’ +1 — 5 ^ c n a; n+í ’ = 0. 

n— 0 77=0 77=0 77=0 

Para que la solución tenga la forma deseada, todas las potencias deben ser iguales 
a n + r, luego la tercera suma la tenemos que reescribir. Realizando los cambios 
necesarios tenemos 

OO OO OO OO 

2 ^(n + r — l)(r¿ + r)c n x n+r — y^(n + r)c n x n+r + ^ c n _ix” +r — 5 ^ c n x n+r = 0. 

77.—0 77=0 77=1 77=0 


Ahora necesitamos que todas las sumas comiencen desde un mismo valor. Es 
decir, todas las sumas deben comenzar desde n — 1. 

OO 

2(r — l)rco:r r + 2 + r — l)(n + r)c n x n+r — rcoX r 

77=1 

OO 

+ Cn-lX U+r - 5 C 0 X r 

77=1 

Agrupando términos, tenemos 

OO 

(2r 2 — 3 r — 5)c 0 a; r + ^[(2(n + r — l)(n + r) — (n + r) — 5)c n + c n -i]x n+r = 0. 

77=1 

Para que esta ecuación se cumpla, recordando que co ^ 0, se debe tener: 

{ 2r 2 — 3r — 5 = 0, 

(2 (n + r — l)(n + r) — (n + r) — 5 )c n + c n _i = 0, para n > 1. 

Observación 2.29. La primera ecuación se le suele llamar ecuación índice asociada 
al punto singular regular xq. 


OO 

~ y^Á n + f)c n x nJrr 

77=1 

OO 

- 5 °nX n+r = 0. 

n= 1 


Las soluciones de la ecuación índice son r\ — - y r 2 — — 1. Reemplazando el valor 

ri en la segunda ecuación obtenemos la ecuación n(2n + 7)c n + c n _i = 0 (realizar los 

cálculos), válida para todo n > 1. Luego, el valor de las constantes c n están dadas 

por c n = - 1 —-, para todo n > 1. Algunos de sus valores son: 

n(2n + 7) 
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n 

Cn 

1 

Co 

ci = - 9 ’ 

2 

Ci Co 

C2 _ _ 22 “ 198 ’ 

3 

C2 Co 

C3 “ 39 _ 7722 


Reemplazando estos valores en la ecuación 2.48, obtenemos una solución de la 
ecuación diferencial. En efecto, tenemos: 


yi = c 0 x 5//2 + Cix 7 / 2 + c 2 x 9 ^ 2 + c 3 x 11//2 + ... 

= CnX 5 / 2 - C 4x V2 + ^X 9 / 2 - ^-X 11 ' 2 + . . . 


= C 0 \X 


, 5/2 


X 


7/2 


+ 


198 

x 9/2 


7722 

x ll/2 


198 7722 


+ ... 


Ahora, reemplazamos el valor de r 2 = — 1 en la ecuación 


2 (n + r — l)(n + r) — (n + r) — 5)c n + c n _i = 0, 


válida para n > 1. Despejando el valor de c n obtenemos una fórmula recursiva para 
estos coeficientes. En efecto, tenemos que 

_ Cn— 1 

Cn = ~2^7’ 


para todo n > 1. 

Algunos valores de c n son: 
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n 


c 


n 


i 



2 



Co 

15’ 


3 


c 3 = c 2 = 


Co 

15 


Reemplazando estos valores en la ecuación 2.48, obtenemos una nueva solución de la 
ecuación diferencial. En efecto, tenemos: 

Vi = — T ci T c 2 x + c 3 x 1 2 + ... 
x 

co ci c 2 c 0 2 

=--- x H- x + ... 

x 5 15 15 

(l 1 x x 2 
= Co U + 5 + l5 + Í5 + ’” 

Observación 2.30. En el ejemplo 2.33, el método de Frobenius nos proporciona dos 
soluciones linealmente independientes, pero este no es siempre el caso. Así, surgen 
de manera natural las siguientes preguntas: 


1. Bajo qué condiciones podemos asegurar que la ecuación 

. . . dy d 2 y 

a 0 (x)y + ai(x)— + a 2 ( x = 0 

tiene dos soluciones linealmente independientes de la forma 


(2.26) 


I X - x 0 r ^C n (x - x 0 ) r 

n =0 

alrededor del punto singular regular xq. 


(2.37) 


2. Si la ecuación diferencial 2.26 no tiene dos soluciones linealmente 
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El siguiente teorema da respuesta a las dos preguntas anteriores. 

Teorema 2.10. Sea Xq un punto singular regular de la ecuación 2.26 y sean r 1; r 2 
(donde, podemos suponer, Re(ri) > Re(r 2 ) ) soluciones del la ecuación índice asociada 
al punto Xo- 

1. Si r\ — r 2 ^ N, entonces la ecuación 2.26 tiene dos soluciones linealmente 
independientes y\, y 2 cuya forma es 

OO 

7 /i =| x - x 0 \ ri ^2 Cn ( x ~ ^o)", 

71=0 

donde c 0 7 ^ 0 , y 

OO 

y 2 =| x - x 0 p ^2 c u( x ~ x o)"’ 

n=0 

donde Cq 7 ^ 0 . 

2. Si Ti — r 2 G N, entonces la ecuación 2.26 tiene dos soluciones linealmente 
independientes y 1 , y 2 cuya forma es 

OO 

7/i =| X - X 0 | ri ^ °n( X ~ X oT , 

71=0 

donde Co 7 ^ 0 , y 

OO 

y 2 =\x-xq I'' 2 ^ c* (x - x 0 ) n + Cyí(x) ln(x - x 0 ), 

71=0 

donde Cq 7 ^ 0 y C es una constante que puede ser cero. 

3. Si ri — r 2 = 0, entonces la ecuación 2.26 tiene dos soluciones linealmente 
independientes y\, y 2 cuya forma es 

OO 

yi=\x-x 0 | ri ^2 c "( x ~ x o) n , 

n=0 

donde Co 7 ^ 0 , y 

OO 

y 2 =| x - x ü | n+1 - x 0 ) n + 7/1 (x) ln(x - x 0 ). 

71=0 

Las soluciones en los Ítems 1, 2 y 3 son todas válidas en el intervalo 0 <| x — xo |< R 
alrededor de Xq. 
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Para terminar esta parte observamos que: siempre se pueden realizar los cálculos 
suponiendo que Xq = 0. En efecto, cuando el punto singular regular es diferente de 
cero, utilizamos la sustitución t = x — x o para trasladar al origen todos los cálculos. 

Ej ercicios 

El siguiente bloque de ejercicios han sido tomados de Ross [RS]. 

Emplear el método de Frobenius para hallar soluciones en un entorno de x = 0 
en cada uno de los siguientes ejercicios: 








2.5 Aplicaciones de las ecuaciones diferenciales 

lineales de orden superior 

Quizá los problemas de naturaleza mecánica sean los más claros ejemplos de aplicación 
de una ecuación diferencial de orden superior. En esta sección estudiamos muy 
someramente el comportamiento de un cuerpo sujeto al extremo de un resorte cuanto 
se le proporciona a este una fuerza. Más aplicaciones se pueden encontrar en los 
circuitos eléctricos, la flexión de vigas, etc. 
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Vibraciones mecánicas 

El objetivo de esta parte es determinar la posición de un cuerpo sujeto al extremo de 
un resorte cuando este abandona su posición de equilibrio como efecto de aplicarle 
una fuerza al sistema resorte-cuerpo. Impondrémos una hipótesis (bastante fuerte) 
sobre el posible movimiento del cuerpo para que las ecuaciones sean más sencillas de 
construir y trabajar. 

Consideremos un resorte de longitud L que se encuentra suspendido verticalmente 
desde el techo (ver figura 2.1). 



Figura 2.1: Resorte suspendido desde un techo, longitud natural L 


Coloquemos un cuerpo de masa m sujeto a un extremo del resorte y dejemos que 
alcance nuevamente la posición de equilibrio (ver figura 2.2). 





L 

m 


1 





Figura 2.2: Masa en equilibrio, longitud del resorte L + l 
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Finalmente saquemos de la posición de equilibrio al sistema resorte-cuerpo 
estirándolo o comprimiéndolo una longitud x (ver figura 2.3). 



Figura 2.3: Masa a una distancia x por debajo de la posición de equilibrio; longitud 
del resorte L + l + x 

Observación 2.31. El modelo matemático que vamos a construir para describir 
la posición de un cuerpo de masa m considera que el movimiento de este ocurre 
únicamente en un dimensión (vertical). Esta hipótesis es bastante fuerte pero muy 
importante, pues nos permite construir las ecuaciones de forma mucha más sencilla. 
Además, consideramos que la dirección positiva es hacia abajo. 

Para construir el modelo matemático utilizamos la segunda ley de Newton: u La 
sumatoria de todas las fuerzas que actúan sobre un cuerpo es igual a la masa del 
cuerpo por la aceleración de este." 

A continuación enunciamos las fuerzas que actúan sobre el cuerpo, determinando 
en cada caso su valor: 

1. Fi, la fuerza de la gravedad. Si g es la aceleración de la gravedad (ya que esta 
actúa en la dirección positiva, hacia abajo), entonces F l = mg. 

2. F-¿, la fuerza de resistencia del resorte. Según la ley de Hooke, la fuerza que actúa 
sobre un cuerpo es directamente proporcional a la elongación (compresión) del 


118 








Leónidas Cerda y Janneth Morocho 


resorte. Como esta fuerza siempre se opone al movimiento del cuerpo tenemos 
F -2 = —k{x + l) -ver figura 2.3. Notamos que si x = 0 (el cuerpo se encuentra 
en posición de equilibrio), entonces la fuerza F 2 debe ser igual a la fuerza de la 
gravedad y su dirección es hacia arriba, luego —mg = —k( 0 + /). Así, mg = kl. 
Reemplazando este valor en la expresión para F- 2 tenemos F 2 = —kx — mg. 

3. F 3 , la fuerza de resistencia del medio. No se conoce con exactitud cuál es 

la magnitud de esta fuerza, pero, para velocidades pequeñas, esta fuerza es 

dx 

proporcional a la velocidad, es decir, F 3 = —a —. donde la constante a > 0 se 

dt 

llama constante de amortiguamiento. 

4. F 4 . cualquier fuerza externa que actúa sobre el cuerpo. Como no conocemos 
con exactitud cual es la naturaleza de esta fuerza externa, nos permitimos 
denotarla con F(t). 


A partir de la segunda ley de Newton tenemos: 


E* 


i = rn 


i= 1 


d 2 x 

dP’ 


reemplazando los valores de cada descritas arriba, tenemos 

j2, 4 


d 2 x 

m ~dF = 




i =1 


dx 


d 2 x dx . 

+ a ~rr + kx — F (t). 
dt z dt 


= mg — kx — mg — a— + F(t) 

, dx . 

= —kx — a ——b F(t). 

CIL 

Así, la ecuación diferencial (modelo matemático) que describe la posición de un 
cuerpo de masa m en cualquier instante t está dada por 

(2.38) 

Si además imponemos las condiciones x(t 0 ) = x 0 y x'{t\) = aq obtenemos el siguiente 
problema (con condiciones iniciales o con condiciones en la frontera dependiendo de 
los valores de to y ti): 

d 2 x dx . . 

m—— + a— + kx = t (t) 
dt ¿ dt 

x(t 0 ) = x 0 
x'{ti) = Xi. 
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Observación 2.32. Cuando F(t ) = 0 se dice que el cuerpo tiene un movimiento 
libre; en otro caso, se dice que el movimiento es forzado. Además, cuando la constante 
a = 0 se dice que el movimiento es no amortiguado y cuando a ^ 0 se dice que el 
movimiento es amortiguado. 

Naturalmente, se pueden tener combinaciones de los tipos de movimientos que se 
menciona en la observación 2.32. Por ejemplo, se puede hablar de un movimiento 
libre amortiguado o de un movimiento forzado no amortiguado, etc. en cada caso la 
ecuación 2.38 toma una forma específica. 

Para finalizar esta parte, mencionamos que, para hallar la posición de un 
cuerpo que tiene un movimiento forzado, primero tenemos que hallar la función 
complementaria del cuerpo que tiene un movimiento libre. 

Circuitos electrónicos 

En esta parte construimos un modelo matemático para estudiar el comportamiento 
de la corriente i que fluye a travéz de un circuito en serie como el que se muestra en 
la figura 2.4: 


R 


WW\A 



c 


Figura 2.4: Circuito en serie 


En la Figura 2.4 se utiliza la siguiente simbología: 


Para la fuerza electromotriz (pila o generador) de voltaje E 
—VWV— R Para un resistor de resistencia R 

l- Para un inductor de inductancia L 




Para un capacitor de capacitancia C 
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Recordamos que la caída de voltaje E a (donde a puede ser el resistor, el capacitor o 
el inductor) a través de cada elemento del circuito está dada por: 


1 . 

2 . 

3. 


Eji = Ri, donde R es la resistencia del resistor. 


Ec = —q, donde C es la capacitancia del capacitor y q es la carga instantánea 


del capacitor. 


dx 

Er = L —, donde L es la inductancia del inductor. 
dt 


Por la primera ley de Kirchhoff: “La suma algebraica de las caídas de voltage a través 
de un circuito cerrado en una dirección específica es igual a cero.” y puesto de las 
caídas de voltaje a través del resistor, capacitor e inductor tienen signo opuesto al 
generado por la fuerza electromotriz se tiene que: 

„ 1 ^ di „ 

Rl+ c q+ 5 


Además, recordando que la corriente i se define como la derivada de la carga con 
respecto al tiempo, la última ecuación se puede escribir como: 


T d 2 q dq 1 

L de + R dt + c q ~ E ■ 


(2.39) 


La ecuación 2.39 es una ecuación diferencial de segundo orden donde la variable 
independiente (incógnita) es la carga q. 


Para averiguar cómo es el comportamiento de la corriente i (una vez resuelta la 
ecuación 2.39), tenemos que derivar la carga q. Es decir, 


i(t) = 


dq 

dt 


Flexión de vigas 

Consideremos el siguiente problema: 

Problema 2.1. Determinar la flexión de una viga rectangular sometida a una carga. 


Para hallar la solución al problema 2.1 consideramos las siguientes hipótesis: 
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1. Inicialmente la viga es recta y su eje central coincide con el eje X (ver figura 
2.5). 

2. El eje central de la viga se flexiona debido a la acción de una carga (suma de 
fuerzas aplicadas a la viga) -ver figura 2.6. 

3. Todas las fuerzas que se aplican a la viga están sobre un mismo plano, el cual 
además contiene al eje central. 



Figura 2.5: Viga horizontal 



Figura 2.6: Aplicación de una carga a una viga 


La curva punteada en la figura 2.6 (la deformación del eje central), llamada curva 
elástica, nos indica cuánto se ha deformado la viga. Por lo tanto, se desea obtener la 
ecuación de la curva elástica. 


Sea M (x) el momento flexionante en una sección transversal vertical de la viga 
en x. En mecánica se demuestra que el momento flexionante en x , debido a todas las 
fuerzas exteriores que actúan sobre la viga (ver figura 2.7), está dado por 


M(x) 


El 

~R' 


122 


















Leónidas Cerda y Janneth Morocho 


donde E es el módulo de elasticidad de Young que depende del material y del diseño 
de la viga, / es el momento de inercia de la sección transversal de la viga en a; y i? 
es el radio de curvatura de la curva elástica. 

y 


Fi F 2 1 

"3 

Fi J 

M(x) 









X 


Figura 2.7: Momento flexionante en x 


[1 + (y') 2 f 2 

El radio de curvatura R estás dado por R = -. Si asumimos que la 

y" 

viga se dobla solo levemente, y' del radio de curvatura es tan pequeño que su cuadrado 
es despreciable con respecto a 1; es decir, para fines prácticos, se puede considerar 
que R = —. Reemplazando este valor en la expresión para M(x ) se tiene que: 

y" 


M(x) = EI C p f. 

dx ¿ 


(2.40) 


La ecuación 2.40 es una ecuación diferencial completa de segundo orden con 
coeficientes constantes que modela la flexión de una viga. 


2.6 Ejercicios del capítulo 2 


1. Hallar la solución general de las siguientes ecuaciones homogéneas: 


a) 

b ) 

c) 

d) 


upL + + ll8 P - 68 ? - 2 % = o. 

ax° dx 4 dx ó dx z dx 

d 3 y d 2 y dy 

~r~% + + j h 3y = 0. 

dx ó dx ¿ dx 

d Ay 0 d 2 y dy = 

dx 4 dx 2 dx 

2—^ + 5—|-7 ^- + 3y = 0. Sugerencia: Usar el método de Cardano para 

dx 6 dx 2 dx 

hallar las raíces del polinomio característico. 
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2. Hallar la solución general de la ecuación —+ 3—^ + —^ + 3 y = xe 2x . 

ax ó dx 2 dx 

C fy dy 

3. Hallar la integral particular de la ecuación 3—- + 5—-2 y — ln x. 

dx 1 dx 

4. Hallar la solución general, en series de potencial alrededor del punto Xo = 1, de 

, •. d 2 y dy 

la ecuación x—— + -—h ¿y = (J. 
dx ¿ dx 

5. Hallar dos soluciones (en series de potencias) linealmente independientes, 

alrededor de x 0 = 0, de la ecuación x 2 —^ + x \ x -]—) — + xy = 0. 

dx ¿ \ 2 J dx 


Los siguientes ejercicios se han tomado de O’Neil [OP] 


6 . Una pesa de 16 libras queda suspendida de un resorte, alargándolo 6 pulgadas. 
La pesa es jalada otras 3 pulgadas abajo de esta posición de equilibrio y liberada. 
En este instante, se aplica al sistema una fuerza igual a - cos( 6 í). Calcule y 
trace una gráfica de la función de desplazamiento, suponiendo que no hay 
amortiguamiento. 

7. Calcule el movimiento de la pesa del problema 6 si hay una fuerza de 
amortiguamiento de 8 v libras, donde v es la velocidad de la pesa. 


En cada uno de los siguientes ejercicios obtenga la corriente i(t) en el circuito RLC 
que se muestra a continuación, con los valores de i?, L y C y suponiendo que una 
corriente inicial y carga del capacitor iguales a cero. 



W 




c 


8 . R = 400 ohms, L = 0.12 lienry, C = 0.004 farad, E(t) = 120sen(20í) volts. 

9. R = 200 ohms, L = 0.1 henry, C = 0.006 farad, E(t) = íe -í volts. 

10. R = 450 ohms, L = 0.95 henry, C = 0.007 farad, E(t) = e -í volts. 

11. R = 150 ohms, L = 0.2 henry, C = 0.05 farad, E(t) = 1 — e -í volts. 
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Capítulo 3 


Sistemas de ecuaciones diferenciales lineales 

En el capítulo 1 y el capítulo 2, estudiamos algunos fenómenos que se pueden 
expresar por medio de una sola cantidad (ver sección 1.1 y sección 2.5). Por ejemplo, 
la población P(t) se expresa por medio de una sola cantidad. Además, vimos que 
para estudiar la razón de cambio de esta población necesitamos únicamente algunos 
supuestos sobre la misma población. 

Para fenómenos más complejos, necesitamos más de una cantidad para estudiar 
su razón de cambio. Por ejemplo, para estudiar la razón de cambio de una población 
de zorros necesitamos saber cual es el tamaño de la población de presas (por ejemplo, 
cual es la población de conejos) y la disponibilidad de alimento. 

Por un lado, el estudio de los fenómenos cuya razón de cambio depende de una 
sola cantidad nos condujo a la idea de ecuación diferencial ordinaria de n-ésimo orden 
(donde n > 1). Por otro lado, el estudio de fenómenos cuya razón de cambio depende 
de más de una cantidad nos llevará a la idea de sistema de ecuaciones diferenciales. 
El siguiente ejemplo nos ilustra un fenómeno de este último tipo. 

Consideremos un ecosistema donde coexisten dos especies con la peculiaridad de 
que una especie se come a la otra. Es decir, una especie es la depredadora y la otra 
especie es la presa. 

Las hipótesis que vamos a considerar para construir nuestro modelo son las 
siguientes: 

1. Si no hay depredadores, entonces la población de presas crece a una tasa 
proporcional a su población y, además esta población no está afectada por la 
sobrepoblación. 

2. La tasa a la que las presas son devoradas es proporcional a la tasa a la que los 
depredadores y las presas interactúan. 

3. Sin la presencia de presas, la población de depredadores disminuye a una razón 
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proporcional a ella misma. 

4. La tasa de nacimientos de depredadores es proporcional al número de presas 
devoradas que, por la hipótesis 2 , es proporcional a la razón a la que las presas 
y los depredadores interactúan. 

Para la formulación del modelo matemático utilizamos las siguientes variables y 
constantes: 

D es la población de depredadores al tiempo t. 

P es la población de presas al instante t. 

a es el coeficiente de la razón de crecimiento de las presas. 

¡3 es el coeficiente de la razón de muerte de los depredadores. 

7 es el coeficiente de proporcionalidad de las interacciones entre el 
depredador y la presa (en las que la presa es devorada). 

5 es la constante de proporcionalidad del beneficio de la población 
de depredadores por una presa devorada. 

Suponemos que los coeficientes a, ¡3, 7 y ó son todos positivos. Desde la hipótesis 
1 concluimos que la variación de la población de presas contiene al término aP. 
Desde la hipótesis 3, la población de depredadores contiene al término —/3D donde, 
el signo negativo indica que la población decrece. 

Ahora debemos hallar un término que exprese la hipótesis 2 , es decir, un término 
que crezca si D o P aumentan pero, que se anule cuando I) = 0 o I’ = 0. La forma 
más fácil de expresar estas ideas es elegir DP como la expresión para indicar la 
interacción entre los depredadores y las presas. Luego la variación de la población de 
presas está afectada por el término — 7 DP. De la misma manera la variación de la 
población de depredadores es afectada por el término 5DP. 

Con estas consideraciones, se tiene que: el modelo matemático (sin resolver) para 
describir la dinámica depredador-presa está dado por 

\ii = aP ~ lDP ' 

\f = . fi D + SD p . 
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Este conjunto formado por dos ecuaciones diferenciales es un ejemplo de lo que se 
conoce como sistema de ecuaciones diferenciales ordinarias. Más específicamente, se 
tiene la siguiente definición: 

Definición 3.1. Se llama sistema de ecuaciones diferenciales a cualquier 
conjunto de dos o más ecuaciones en el cual aparezcan las derivadas de una o 
más variables independientes, las incógnitas de la ecuación, respecto a una o más 
variables dependientes. 

Son ejemplos de sistemas de ecuaciones diferenciales los siguientes: 




El sistema de ecuaciones diferenciales se denomina sistema ordinario cuando 
todas las derivadas que aparecen en el sistema son derivadas ordinarias (el sistema 
1 es un ejemplo de sistema ordinario); en caso contrario, el sistema se denomina 
sistema parcial (el sistema 2 es un ejemplo de sistema parcial). 

Observación 3.1. En este texto tratamos únicamente sistemas ordinarios y cada 
vez que mencionemos sistema de ecuaciones diferenciales entenderemos que se trata 
de un sistema ordinario. 

El objetivo de este capítulo es estudiar sistemas de ecuaciones diferenciales para 
los que se conoce métodos analíticos de solución. 


3.1 Generalidades sobre los sistemas lineales 


Una primera clasificación que se hace a los sistemas de ecuaciones diferenciales es en 
sistemas de ecuaciones diferenciales lineales y sistemas de ecuaciones diferenciales no 


127 





INTRODUCCIÓN A LAS ECUACIONES DIFERENCIALES ORDINARIAS 


lineales. La siguiente definición establece con toda precisión lo que se entiende por 
sistema de ecuaciones diferenciales lineales. 


Definición 3.2. Se llama sistema de ecuaciones diferenciales lineales al siguiente 
conjunto de ecuaciones 


Mu j i M12 

E i / N d X\ \ —> 9 , \ (i Xo 


M\n 


i =0 
M 21 


¿=0 

M22 


dt l 


i =0 i =0 


dt l 


■ + y^ a n( ¿ ) 

i =o 
M2n 

■+J2<Át) dt¡ 


d l x n 

dt l 

d l x n 


i =o 


h(t) 

h (t) 


^ml 77 ^m 2 77 Mmn 77 

E «í-w^ + E + • • • + E = *•<*> 


k 2=0 zO z =0 

Observación 3.2. Para simplificar el lenguaje diremos sistema lineal en vez de 
sistema de ecuaciones diferenciales lineales. 


En el caso particular de = 1 para todo 1 < i < m, 1 < j < n, el sistema se 
llama lineal de primer orden. De esta forma, en un sistema lineal de primer orden, 
cada ecuación contiene a lo máximo la primera derivada (la derivada de orden cero, 
cuando i = 0, corresponde a la misma función). Luego, un sistema lineal de primer 
orden se puede escribir como 








. , dx 1 . . dx 2 

an (t)—+a n ( t )— + .. 

dr EE 

■ + ain(t) + y y bu(t)xj = /i(í) 

i=n 


. . dx 1 , . dx 2 

“ 2,(í) ^r + “ 22(í) ^ó + " 

dr EE 

• + 0>2n{t) ^ b2i(t)Xj — f2 (t) 

i=n 

(3.1) 

. . dx\ , . dx 2 

+ . 

dr EE 

• • + CLmnif) + y^ b mi (t)Xi = f m (t) 



i=n 


Cuando m = n. el sistema 3.1 tiene igual número de ecuaciones y de incógnitas y 
podemos llamarlo cuadrado 1 . Este texto trata únicamente sistemas cuadrados y cada 

: La denominación no es estándar; sin embargo, nosotros utilizamos este nombre para que guarde 
relación con los sistemas algebraicos. 
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vez que mencionemos sistema de ecuaciones diferenciales entenderemos que se trata 
de un sistema cuadrado. 

Ahora clasificamos el sistemas 3.1. 


1. El sistema se llama homogéneo cuando f t (t) = 0 para todo 1 < i < m. 

2. El sistema se llama con coeficientes constantes cuando a^(x) y bij(t) son 
funciones constantes para todo 1 < i < m, 1 < j < n. 

3. El sistema se dice que es un sistema con coeficientes variables cuando alguna 
función a,ij(x ) o bij(t) no es una función constante para algún 1 < i < m, 
1 < j < n. 


Por ejemplo, el sistema 


2 dx dy 

x — - 3— — x — 1 

dt dt 

dx dy 

——b sena;— = 0, 
dy dt 

es un sistema con coeficientes variables no homogéneo. Mientras, el sistema 

9 dx dy 

x 2 --?>4- = 0 


dt 


dt 


dx dy 

——b seno;— = 0, 
dy dt 

es un sistema con coeficientes variables homogéneo. 


Definición 3.3. Se llama solución del sistema 3.1 en el intervalo I a cualquier 
vector x(t) = (oq(í), a; 2 (í), ..., x n (t)) tal que las funciones x*(í) y sus derivadas están 
definidas en el intervalo I, para todo 1 < i < n, y las funciones Xi(t ) satisfacen el 
sistema en el intervalo I. 


Ejemplo 3.1. Comprobar que el vector x(t) = (eos 2t, sen 2t) es solución del sistema 

' dx 


dt 2y 

dy _ 0 

l dt 2 ’ 


en el intervalo I = 
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Solución. Las funciones definidas por x(t) = eos 2 1, y(t) = sen2t están definidas en 
dx dy 

todo t real. Ademas — = — 2sen2í = —2y. — = 2 eos 2 1 = 2x. Luego, las funciones 
dt dt & ’ 

x, y satisfacen el sistema. 

Definición 3.4. Se dice que un sistema de ecuaciones diferenciales es normal cuando 
tiene la forma (o se puede reducir a esta forma) 
dx\ 

~r,-~ — o-u(t)xi + ai 2 Ít)x 2 + ... + a\ n x n (t) + b\{t) 

^ = a 2 i(t)x 1 + a 2 2 {t)x 2 + • • • + a 2 nX n {t ) + b 2 (t) 


dx n 
, dt 


= a nl (t)xi + a n2 (t)x 2 + ... + (í) + bi(t). 


Un sistema normal es importante por su conección con las ecuaciones de orden 
superior. En efecto, consideremos la ecuación de n-ésimo orden 


. . , .dx . . d 2 x d n x , . . 

a 0 (t)x + + a 2 {t)~^p + • • • + = b(t). 


(3.2) 


Realicemos las siguientes sustituciones: 


dx d 2 x 

Xi = x,x 2 = = 


dt 


á¿2 ....,X n _ 1 - dt n-2^n 


d n 2 x d n l x 

x n = 


dt 


n —1 


De esta manera conseguimos la siguiente cadena de derivadas: 

dx dx i d 2 x dx 2 d n ~ 1 x dx n _ i d n x dx 

dt dt ’ dt 2 dt ’ ’ dt 


n —1 


dt dt n 


dt 


Así, podemos considerar el siguiente sistema normal 

r dx i 
dt 

dx 2 


= x 2 


dt 

dXn-l 

dt 


= x 3 


= X r 


dx n 
< dt 


= b n - a 0 (t)xi - a l {t)x 2 - ... - a n _i(t)x n , 

cuya solución x{t) = (xi(t),X 2 {t),... ,x n (t)) es claramente solución de la ecuación 
diferencial de n-ésimo orden 3.2. 
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3.2 El operador diferencial 

En esta sección estudiamos un método de resolución de sistemas de ecuaciones 
diferenciales con coeficientes constantes. Este método depende fuertemente del uso 
del operador diferencial que introducimos en la siguiente definición: 

Definición 3.5. Se llama operador diferencial D a la función definida como la 
operación de derivación. 

A partir de la definicón 3.5 se tiene que el campo de aplicación del operador 
diferencial D es el conjunto de todas las funciones x que admiten primera derivada 
(sobre algún intervalo I). 

Por ejemplo, D aplicado a la función senx es la función cosa:, es decir, 

dísenx) 

sena:) =---= cosa:. 

dx 

Observación 3.3. Por recursión se puede definir el operador D n . En efecto, se tiene 
que 

D n = D(D íl - 1 ). 

d n f 

A partir de la definición del operador D n se tiene que D n (f ) = —— siempre que 
/ sea una función que admite derivada hasta el orden n. Además, admitimos que 

D°(f) = f. 

Al operador diferencial D n se llama operador diferencial de orden n u operador 
diferencial de n-ésimo orden. 

Sea Q el conjunto de todos los operadores diferenciales. Se define en Q las 
operaciones suma y producto por escalar. En efecto, sean Di, D 2 dos operadores 
diferenciales y A un número real. La suma D\ + D 2 y el producto por escalar A D\ se 
definen de la siguiente manera: 

1. (D 1 + D 2 )(f) = D 1 (f) + D 2 (f). 

2. (XDi)(f) — XDfif). 
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A partir de las operaciones suma y producto por escalar se define el siguiente 
operador diferencial. 

Definición 3.6. Se llama operador diferencial lineal de orden n con coeficientes 
constantes al siguiente operador: 

L: = oq + a\D + a 2 -D 2 + ... + a n D n , (3-3) 

donde G K para todo i G {0,1,..., n}. 

Es inmediato, a partir de la Definición 3.6, que el operador diferencial lineal de 
orden n es un operador lineal, es decir, L(af + ¡3g) = aL(f) + ¡3L(g ). 

Definición 3.7. Sea I un intervalo abierto, n un número natural. Denotamos por 
C n (I) al conjunto de todas las funciones f tales que f y todas sus derivadas de orden 
hasta n son continuas en I. 


En la Definición 3.7 se entiende que C°(I) es el conjunto de todas las funciones 
continuas en I y C°°(I ) es el conjunto de todas las funciones infinitamente derivables 
con continuidad. 

Observación 3.4. Por definición, L es una función con dominio C n (I) y recorrido 

C°(I). 


Sean L\ y L 2 los operadores diferenciales dados por 

L\ = oq T a\D, 

L 2 = b\ D. 


Se tiene que 


{L 1 oL 2 ){f) = L 1 (L 2 (f)) 

= Li(6i£>(/)) 


T [h df 

— L\ ( bi — 
dx 


, df / df 

— o-obi ~—b a\D I b\ —— 

dx \ dx 

, df , d 2 f 

— —b afiji—— 

dx dx ¿ 

— iflobiD + aib 1 D 2 )(f). 
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Se puede comprobar fácilmente que L\ o Lo = L 2 ° L\. 

Las operaciones realizadas arriba nos muestran que para hallar L\ o L 2 podemos 
multiplicar L\ con L 2 como si fueran polinomios. En realidad se puede demostrar el 
siguiente teorema. 

Teorema 3.1. Sean Li = a 0 + a\D + ... + a n D n , L 2 = b 0 + b\D + ... + b m D m dos 
operadores diferenciales lineales. La composición del operador Li con el operador L 2 
está dada por 

m n 

a.¡b.jD i+3 . 

j =0 i =0 

Observación 3.5. El teorema 3.1 nos dice que la composición entre operadores 
tiene el mismo comportamiento que la multiplicación entre polinomios, luego podemos 
escribir LiL 2 en lugar de Li o L 2 . 

La analogía que existe entre la composición de operadores diferenciales lineales 
con coeficientes constantes y la multiplicación de polinomios nos permite factorar un 
operador diferencial. En efecto, el operador 3.3 se puede factorar como 

L = a n (D - n)(D - r 2 ) ■ ■ ■ (D - r n ), 

donde r¿ para 1 < i < n son las raíces del polinomio ao + Oír + + ... + a n r n . 

En la práctica, cuando existen raíces complejas, digamos r 3 y rj, se las combina 
(multiplica) para formar un polinomio de segundo grado irreducible. 

Un método operacional para resolver sistemas lineales 


En esta parte, explicamos como aplicar la noción de operador diferencial a la 
resolución de sistemas lineales de ecuaciones diferenciales con coeficientes constantes. 
Nos limitamos a explicar el método, que ciertamente funciona, abandonando 
completamente el sustento teórico del mismo. 


Consideremos el sistema 


dx dy , , 

M + M- 2X ~ iy = e ’ 


dx dy 
dt dt 


V = e 


„ 4 í 


(3.4) 
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Utilizando el operador diferencial lineal, el sistema 3.4 lo podemos escribir como 


{(D-2)(x) + (D-4)(y)=e t , 
[D(x) + (D- 1 )(y) = e 4í . 


(3.5) 


Trabajamos ahora en el sistema 3.5. Multiplicando por D la primera ecuación y 
por D — 2 la segunda ecuación, obtenemos el sistema 

í D{D - 2){x) + D(D - 4) (y) = D{é), 

(^(U - 2)(x) + (D - 1)(D - 2 )(y) = (D — 2)(e 4í ). 

Restamos la primera ecuación con la segunda para obtener (—D — 2)(y) = e ¿ —2e 4í , 
la cual representa la ecuación diferencial 

§ + 2 a = 2e 4 ‘ - e*. (3,6) 


La ecuación 3.6 es una ecuación diferencial lineal de primer orden cuya solución 

es y = -e 4í — -e t + kie~ 2t (ver sección 1.4, pag. 52) donde k\ es una constante 
O o 

arbitraria. 


Ahora multiplicando por D — 1 la primera ecuación y por D — 4 la segunda 
ecuación obtenemos el sistema 

í (D - l)(D - 2)(x) + (D - l)(D - 4 )(y) = (D - l)(e í ), 

[l)(D - 4)(x) + (D — 1 )(D - 4)(y) = (D - 4)(e 4í ). 

Restamos la primera ecuación de la segunda para obtener (D + 2) (re) = 0, la cual 
representa la ecuación diferencial 

+ 2x = 0. (3,7) 

La ecuación 3.7 es una ecuación diferencial en variables separables cuya solución 
es x = k 2 e~ l (ver sección 1.2, pag. 34) donde k-¿ es una constante arbitraria. 

A este punto podemos verificar que las funciones x, y halladas con este 
procedimiento aún no satisfacen el sistema 3.4 pero, esta inperfección se debe a que 
estas funciones contiene demasiadas constantes arbitrarias. 
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Para averiguar cuantas constantes arbitrarias deben contener las soluciones del 
sistema calculamos el determinante 


A 


D-2 D-4 
D D-l 


Realizando el cálculo del determinante se tiene que A = D + 2. Resulta que 
A es un operador diferencial de primer orden. Luego debe existir únicamente una 
constante arbitraria en las soluciones del sistema de ecuaciones diferenciales. 

Para hallar la única constante arbitraria reemplazamos los valores hallados para 
x, y en cada una de las ecuaciones del sistema 3.4. 

Desde la primera ecuación, obtenemos (realizar los cálculos) 


Qk\ + Ako = 0 . 


(3.8) 


A partir de la segunda ecuación obtenemos la ecuación (realizar los cálculos) 


3fci + 2k 2 = 0 . 


(3.9) 


Las ecuaciones 3.8 y 3.9 constituyen un sistema de ecuaciones en las incógnitas 
ki, k 2 . Además, la matriz de coeficientes tiene rango 1, luego el sistema tiene infinitas 
soluciones, en particular una de las incógnitas se puede expresar en términos de la 
otra. 

3 

Desde la ecuación 3.9 (o desde la ecuación 3.8) se tiene que k 2 = — -k\. Así, las 
soluciones del sistema 3.4 están dadas por 

x = —-/cie _£ , y = -e 4í - e t + £qe~ 2í . 

2 ’ y 3 3 


Ejemplo 3.2. Hallar la solución general del siguiente sistema de ecuaciones 
diferenciales lineales: 


{ dx 

di 

dx 

dt 


du 

+ dt +2y = Sení ’ 

, d y n 

+ -- x - y = 0 . 
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Solución. Escribimos el sistema utilizando el operador diferencial. 

D(x) + (D + 2 )(y) = sen t, 

(D — l)(x) + (D — l)(y) = 0. 

Multiplicamos la primera ecuación por D—1 y la segunda ecuación por D+ 2, luego 
restamos la primera ecuación de la segunda y obtenemos (—2D + 2)(x) = eos x — sen x. 
Claramente, esta expresión representa la ecuación de primer orden lineal 

dx sen x — eos x 

di ~ X = 2 


(3.10) 


Ahora multiplicamos por D — 1 la primera ecuación y por D la segunda ecuación. 
Restando la primera ecuación de la segunda se obtiene (2 D — 2) (y) = eos x — sen x, 
esta expresión representa la ecuación lineal de primer orden 

dy eos x — sen x 


(3.11) 


Las ecuaciones 3.10 y 3.11 tienen solución x = 


sení k\e l sen t k^e 1 

~2~ + ~2~’ y = ~Y~ + ~Y~ 

respectivamente. Nos falta saber cuantas constantes arbitrarias tienen las soluciones 
(por el momento aparecen las constantes arbitrarias k\ y k-¿)- Con el propósito 
de averiguarlo calculamos el determinante cuya primera fila está formada por los 
operadores diferenciales de la primera ecuación y cuya segunda fila está formada por 
los operadores de la segunda ecuación. Es decir, calculamos el determinante 


A = 


D D+2 
D-l D-l 


= -2D + 2. 


Puesto que este es un operador de orden 1, las soluciones deben tener una sola 
constante arbitraria (esto significa que las soluciones deben compartir la misma 
constante arbitraria). Para obtener los valores de k \, reemplazamos las ecuaciones 
de x, y en el sistema. 


Desde la primera ecuación del sistema se obtiene la ecuación k\ + = 0. 

La segunda ecuación no da ninguna información sobre los valores de k\ y fc 2 
(reemplazando se tiene 0 = 0, que es una verdad absoluta pero, no da información 
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sobre k\ y fc 2 )- Luego se tiene que fc 2 = —k\. Finalmente, reemplazamos el valor de 
k <2 en las ecuaciones que obtuvimos para x, y. 


x = 


sení k\e l 


y = 


sen t k\e t 


Ej ercicios 


Los siguientes ejercicios han sido tomados de Ross [RS]. 

Utilizar el método operacional para hallar la solución general en cada uno de los 
sistemas siguiente: 


1 . 


2 . 


dx dy , 

Tt + Tt~ 2x - iv = e ' 

dx dy 

Tt + Tt~ y = e 

dx dy 

Tt + H~ x = - 2t ' 

dx dy 9 

d¡ + d¿~ 3x - v = t ' 

„dx dy . 

2 ~dt + ~dt~ X ~ y ^ e ’ 

dx dy . 

- + - + 2x + y = e. 

dx dy 

- + - + 2 y = sent , 
dx dy 

dt + dt~ X ~ y = 0 ' 

„ dx dy 

2 dt + dt ~ X ~ y= ’ 
dx dy 

Tt + Tt Jr2x - y = L 


3.3 Método matricial para sistemas normales 
homogéneos con coeficientes constantes 


En esta sección, explicamos cómo utilizar matrices para resolver un sistema normal 
homogéneo con coeficientes constantes (ver definición 3.4, pag. 130). Es decir, sistemas 
normales para los cuales las funciones son nulas para todo 1 < i < n. 
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Suponemos que el lector tiene conocimiento de la teoría básica de matrices. Sin 
embargo, comenzamos esta sección recordando algunos conceptos y resultados básicos 
sobre matrices. 

Conceptos básicos de la teoría de matrcices 


Se llama matriz de orden n x m a cualquier arreglo de n filas con m columnas 2 . Para 
una matriz A de orden n x m, utilizamos la siguiente representación: 


«11 

«12 

• • n 

«21 

022 

• • ^2 n 

«mi 

«m2 • 1 

■ • ^mn 


Los números reales (o funciones reales) a¿j se llaman entradas de la matriz donde, 
el subíndice i indica la fila y el subíndice j indica la columna, luego 1 < i < n, 
1 < j < m. 

Es costumbre denotar a la matriz A simplemente con A = ( Oij) mxn o con A = (a¿j) 
cuando está claro cual es el orden de la matriz. 

Clases especiales de matrices son las llamadas matrices cuadradas, la matriz 
traspuesta de una matriz dada y las matrices columna cuyas definiciones son como 
sigue: 


Definición 3.8. 1. Se dice que una matriz A es cuadrada cuando m = n. Es 

decir, una matriz de la forma 


A 


a ii « 12 
«21 022 


Oln 

02n 


Onl O n 2 


2. Sea A = (a¿j) mxn una matriz. Se llama matriz traspuesta de A a la matriz 


2 Aunque esta noción de matriz es correcta, es demasiado general para nuestros propósitos. En 
este texto, las filas y columnas de una matriz están formadas únicamente por números reales o 
funciones reales. 
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A T = ( aji) nX m■ Es decir, a la matriz que se obtiene de A cambiando sus filas 
por sus columnas. 


3. Se dice que una matriz B es una matriz columna cuando n = 1. Es decir, una 
matriz de la forma 


B = 


a n 
«21 

®ml 


Además de las matrices definidas anteriormente se puede mencionar a la matriz 
identidad y a la matriz nula. La matriz cuadrada / = (a¿j) se dice matriz identidad 

si üij = 1 cuando i = j y = 0 en otro caso. Se llama matriz nula a la matriz cuyas 

entradas son todas iguales a cero. 

Observación 3.6. Una matriz columna de orden m x 1 se llama vector de orden m 
o simplemente vector cuando no haya ambigüedad sobre el orden de este vector. 

Definición 3.9. Sea A = (a¿j) una matriz de orden m x n. Se dice que la matiz 

A es constante si üíj G R para todo 1 < i < m, 1 < j < n. Si la matriz A no es 

constante se dice que es una matriz función. 


Son ejemplos de matrices las siguientes: 



1 

-4 

3 

2 


t -1 

-í 2 + 5 

sen t + 7 

lní 

A = 

-5 

0 

7 

-3 

, ^0) = 

t 2 

t — 2 

tañí 

secí 


1 

0 

-3 

10 


i -t t 2 +1 -1 

1 

t 

t 


A es una matriz constante y tp(t) es una matriz función. 


Observación 3.7. Sean tp(t) una matriz función, c G R. Se definen las matrices 
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d<p(t) 

dt 


d<p(t) 

dt 


de <¿?(t) ; 

y J v(u)du 

¥>ii (*) 

(p' 12 (t) ... 

^21 (0 

^22 {t) 

V'ml (0 

^m 2 (í) ••• 


tflrSt) 

<P2 ni*) 
V'mnif) 


ip(u)du = 


rt rt 

ipn(u)du / <¿>i 2 (u)du 


<¿>21 (u)du / (p 2 2 Íu)du 


ipin(u)du 


tp 2 n(u)du 


J </>mi(íí)<¿w J íp rn 2 {u)du ... J ip mn (u)du 
Sea Ai el conjunto de todas las matrices, se define las siguientes operaciones en 


M: 


1. Suma de matrices. Sean A = () y B = ( b^ ) dos matrices del mismo orden. 
La suma de A con B se define como 

A + B = ( a,ij + bij). 

2 . Producto de una matriz con un número real. Sea A = [a l j) una matriz y A un 
número real. El producto de A con A se define como 


XA — (Xdij). 


3. Producto de matrices. Sean A = (a^) y B = () dos matrices de orden m x n 
y n x p respectivamente. El producto de A con B se define como 

n 

AB (c¿j) m xp 5 donde ^ ^ díkbkj- 

k =i 

Definición 3.10. Sea A una matriz cuadrada de orden n. Se dice que A es invertible 
si existe una matriz B tal que AB = BA = /. 
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Naturalmente la matriz B de la definición 3.10 tiene que ser del mismo orden 
que la matriz A. Además, la matriz B se la llama matriz inversa de A y se la denota 
con A -1 . 


Por ejemplo, la matriz A 


1 -4 
5 0 


tiene matriz inversa 


A " 1 


0 1/5 

-1/4 1/20 


y la matriz (p(t) 


sen t eos t 
eos t — sen t 


tiene matriz inversa 


ip l (t) 


sen t eos t 

eos t — sen t 


Observación 3.8. Cuando A es una matriz función, el concepto de matriz inversa 
de A está supeditado a un intervalo I. Es decir, se habla de la matriz inversa de A 
en el intervalo I. Si no existe confución, se dirá simplemente matriz inversa aún en 
el caso de una matriz función. 


Ahora exponemos un método para hallar la matriz inversa de una matriz A. 
Antes necesitamos dar las siguientes definiciones: 

Definición 3.11. Sea A = ( aij) nxn , Mij es la matriz de orden (n — 1) x (n — 1) que 
se obtiene de la matriz A eliminando la i-ésima fila y la j-ésima columna de A. El 
determinante de la matriz A, denotado como det A, se define recursivamente como 
sigue: 


• Si n — 1, entonces det A = a n , 


n 


• Si n > 2, entonces det A = (—1)q^- det M t j, 

3= 1 

arbitrario pero permanece fijo en todo el cálculo. 


donde el subíndice i es 


Ejemplo 3.3. Calcular el determinante de la matriz A = 


a n «i 2 
«21 a 22 
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Solución. Para realizar el cálculo del determinante de A fijemos la primera fila 3 . Es 
decir mantenemos fijo el valor de i — 1. Se tiene que Mu = [ 022 ], M 12 = [ 021 ]. Luego 

2 

det A = l) 1+J aij det M\j 
3 = 1 

= (—l) 1 + 1 a,n det Mu + ( — l) 1 + 2 Oi 2 Mi 2 
= 011022 — a '12 a 21- 

El ejemplo anterior muestra la forma clásica de calcular el determinante de una 
matriz cuadrada de orden 2 . 

Observación 3.9. En realidad, el determinante es una función cuyo dominio es 
el conjunto de todas las matrices cuadradas y cuyo recorrido es el conjunto de los 
números reales. Es decir, si 9Jl c denota el conjunto de todas las matrices cuadradas, 
entonces el determinante es la función definida como: 

det: ÜJl c —y M 
A 1 ^ det A, 

donde det A es como en la definición 3.11. 

Definición 3.12. Sea A una matriz cuadrada de orden n, a¿j el elemento de la 
i-ésima fila y la j-ésima columna de A. 


1. El menor de aij, denotado rriij, es el determinante de la matriz Mij (ver 
definición 3.11). 


2. El cofactor c^- de a y está definido como Cij = 

3. La matriz de los cofactores de los elementos de A, denotada cof A, está definida 
como la matriz que se obtiene de A reemplazando los elementos a t] por 0 ^ para 


1 < i,j < n. 


Ejemplo 3.4. Sea A = 


2 1-1 
4 3-2 

-6 2 5 


hallar cof A. 


3 Se deja como ejercicio comprobar que: fijando la segunda fila, se obtiene el mismo resultado. 
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Solución. Primero realizamos el cálculo del cofactor de an: 


cu = (—l) 2 mn 


= det 


3 

2 


-2 

5 


= 19. 


De la misma manera vemos que (realizar los cálculos) C \2 — ~ 8 , C 13 = 26, C 21 — 7, 
c 22 = 4, c 23 = -10, c 3 i = 1, c 32 = 0, c 33 = 2. Luego 


cof A = 


19 -8 26 

-7 4 -10 

10 2 


Definición 3.13. Sea A una matriz cuadrada, cof A la matriz de los cof actores de 
A. Se llama matriz adjunta de A a la siguiente matriz: 


Adj A = (cof Á) T 


El siguiente teorema nos proporciona una condición necesaria y suficiente para la 
existencia de la inversa de una matriz. Además, nos dice cómo calcular la inversa de 
una matriz (cuando existe). 


Teorema 3.2. Sea A una matriz cuadrada. Si det A 7 ^ 0, entonces existe la inversa 
de A. Además, la matriz inversa de A está dada por 


A - 1 


1 

det A 


Adj A. 


Ejemplo 3.5. Calcular (si existe) la matriz inversa de A 


3 -2 
2 5 


Solución. Primero veamos si tiene sentido calcular A -1 . Como det A = 19, existe 
la matriz inversa de A. Para hallar A -1 primero calculemos los cofactores de A: 
Cu = 5, C 12 = —2, C 21 = 2, C 22 = 3. Con estos valores se tiene que la matriz de los 


cofactores de A es la matriz cof A = 


5 -2 
2 3 


Luego, la matriz adjunta de A es 
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Adj A = 


5 2 
-2 3 


. Finalmente calculamos la matriz inversa 


A " 1 


1 

det A 
1 

19 


Adj A 

5 2 
-2 3 


5/19 2/19 

-2/19 3/19 


Valores y vectores propios 

Sea A una matriz constante de orden n. Cosideremos la ecuación 

Ax = Xx, (3-12) 

donde A es un número real y x es un vector (la incógnita). 

Es trivial ver que el vector nulo es solución de la ecuación 3.12. Estamos 
interesados en dar solución al siguiente problema: 

Problema 3.1. Dada una matriz constante A y un número real A. Hallar una 
solución no trivial de la ecuación Ax = Xx. 


Por ejemplo, el vector 


x = 


-3 

2 


es una solución (no trivial) de la ecuación 3.12 para A = 


6 -3 
2 1 


y A = 4. 


Antes de resolver el problema 3.1 necesitamos dar la siguiente definición: 


Definición 3.14. Un valor propio (o eigenvalor) de una matriz A es un número X 
para el cual la ecuación Ax = Xx tiene solución x no trivial. 


Sea A un valor propio de la matriz A. La solución no nula x de la ecuación 
Ax = Xx se llama vector propio (o eigenvector) asociado a A. Por ejemplo, el vector 
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x = 


-3 


es un vector propio asociado a A 


2 

propio de la matriz A 



4 ya que este valor de A es un valor 


2 1 

Observación 3.10. La definición 3.14 toma en cuenta matrices de orden m x n. 
Sin embargo, en este texto solamente consideramos matrices cuadradas. 


Para hallar los valores propios de una matriz A se procede de la siguiente manera: 
sea Ao un valor propio de la matriz A, luego existe un vector no nulo Xq tal que la 
igualdad Ax o = XoXo = XqIxq (donde / es la matriz identidad) se satisface. Luego 
(. A — ÁoI)xo = 0. Por tanto, el sistema lineal (A — A I)x = 0 tiene al menos dos 
soluciones, a saber: la solución trivial y la solución no nula Xq, luego det(Al — XI) = 0. 
Puesto que det(Al — XI) representa un polinomio de grado n (el orden de la matriz 
A), denotado por pa, se tiene que: los valores propios de la matriz A son todos y 
solamente las soluciones de la ecuación algebraica Pa = 0. 

Observación 3.11. Puesto que pA es un polinomio de grado n se tiene que: una 
matriz de orden n tiene exactamente n valores propios. 


Para hallar los vectores propios asociados al valor propio A se tienen que calcular 
las soluciones no nulas del sistema lineal (A — XI)x = 0. Notamos que existen infinitos 
vectores propios asociados a un único valor propio A. 


Ejemplo 3.6. Hallar los valores y vectores propios de la matriz A = 


1 2 
2 1 


Solución. Primero vemos cual es el polinomio pa'- 


Pa = det(AL — XI) 



= (1- A) 2 -4 
= A 2 - 2A - 3. 
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Ahora calculamos las raíces del polinomio pa, es decir, calculamos las soluciones 
de la ecuación A 2 — 2A — 3 = 0. Realizando los cálculos tenemos Ai = 3, A 2 = —1, 
como los valores propios de la matriz A son las soluciones de la ecuación Pa = 0 
tenemos que los valores propios de la matriz A son 3 y —1. 

Calculamos los vectores propios asociados al valor propio Ai = 3. Para hacer 
clara la exposición, denotamos al vector x por x = (xi,X 2) t - Los vectores propios 
son las soluciones no triviales del sistema lineal (A — 3 1)x = 0. Es decir, los vectores 
propios asociados al valor propio Ai = 3 son las soluciones no nulas del sistema 

J —2xi + 2 x 2 = 0 
I 2 a; 1 — 2x2 = 0 . 

Puesto que las ecuaciones del sistema anterior son similares 4 se tiene que el valor 
de una incógnita depende de la otra. En efecto, tenemos que X\ = X 2 - Luego los 
vectores propios asociados al valor propio Ai = 3 tienen la forma x = (a, a) donde 
cl £ K . \ { 0 ]-. 

De la misma forma (realizar los cálculos) se tiene que los vectores propios asociados 
al valor propio A 2 = — 1 tienen la forma x = (a, —a) donde a € M \ {0}. 

El método matricial 


Ahora presentamos una técnica para hallar las soluciones de un sistema normal de 
ecuaciones diferenciales (ver definición 3.4, pag. 130) que además es homogéneo y 
tiene los coeficientes constantes. Es decir, sistemas de la forma 


" dx 1 


dt 


a n xi + CL12X2 + • • • + a\ n x n {i) 


dx 2 
< dt 


0,2lX\ + 0,22^2 + • • • + d2n%n(t) 


(3.13) 


dx n 
, dt 


1 + ®n2^2 + ... + a nn x n {t). 


4 Dos ecuaciones se dicen similares si y solo si tienen las mismas soluciones. Luego, una ecuación 
se obtiene de la otra multilicándola por un número real diferente de cero. 
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Observamos que: si escribimos 


«11 Üi2 


«1 n 


A 


«21 « 22 


«2 n 



«ni «n2 


entonces el sistema 3.13 se puede escribir como 


ix m A -¿ 


(3.14) 


La ecuación 3.14 se llama ecuación diferencial vectorial asociada al sistema 3.13. 

Observación 3.12. Se puede demostrar que: si ipi, p 2 , ■ ■ ■, p m son vectores solución 
linealmente independientes (ver definición 2.3 del capítulo 2, pag. 66) de la ecuación 
3.14, entonces o¿\p\ + ... + a m p m también es solución de la ecuación 3.1 4- 

La siguiente definición es una generalización de la definición 2.4 del capítulo 2 
(ver pag. 68). 

Definición 3.15. Sean ■ ■ ■ i<Pn vectores solución de la ecuación diferencial 

vectorial 3.14 en el intervalo I (ver definición 3.3). El conjunto {p i, p 2 , ■ ■ ■, <En} se 
llama conjunto fundamental de soluciones de la ecuación 3.14 en el intervalo I si 
y solo si <pi, P 2 -, • i'-Pn son vectores linealmente independientes en el intervalo I y 

cualquier solución de la ecuación diferencial vectorial se puede expresar de la forma 


CiPi + C 2 P 2 + • • • + c n ip n , 


donde ci, C 2 ,..., c n son constantes arbitrarias. Además, la suma anterior se llama 
solución general de la ecuación diferencial 3.14 en el intervalo I. 

Se puede demostrar que el conjunto fundamental de soluciones de la ecuación 
3.14 tiene cardinalidad n. Luego, para hallar la solución general, es suficiente con 
encontrar un conjunto de n soluciones linealmente independientes. 

Ahora ya estamos en condiciones de enunciar el resultado principal de esta 
sección. 
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Teorema 3.3. Consideremos la ecuación diferencial vectorial 


donde A = 


dx(t 

dt 


= Alt, 


«11 

a\2 

021 

022 

Onl 

a n 2 

un valor propio 


^1 n 
n 


y x(t) = [x 1 (t),x 2 {t),... ,x n {t)] T . 


(3.15) 


un vector propio asociado al valor propio X, es solución de la ecuación diferencial 
vectorial. 


En la literatura matemática se puede hallar el siguiente resultado referente a la 
teoría de los valores propios de una matriz y sus correspondientes vectores propios 
asociados: 

Teorema 3.4 (Valores propios diferentes). Sean Ai, A 2 valores propios de la matriz 
A con vectores propios asociados <pi,<p 2 respectivamente. Si Ai 7 ^ A 2 , entonces y 
P 2 son linealmente independientes. 


Naturalmente, cuando tp 1 y <^ 2 son vectores linealmente independientes también 
los vectores ifie M y </ 9 2 e Aí son linealmente independientes para todo A € M. Luego se 
tiene el siguiente corolario: 

Corolario 3.1. Si los valores propios Ai,..., A n de la matriz A son todos diferentes, 
entonces el vector 

n 

<p(t) = 'Y^CiOtie Xit 
i= 1 

es la solución general de la ecuación diferencial vectorial 3.15, donde cp son vectores 
propios asociados a los valores propios A¿ y Ci son constantes arbitrarias. 


Ejemplo 3.7. Hallar la solución general del siguiente sistema: 

' dx 1 

—- = 7x\ - x 2 + 6x3 
dt 

dx2 

—— = — 10xi + 4x9 — 12x 3 
dt 

dx 3 

—— = —2xi + x 2 - x 3 . 
dt 


148 









Leónidas Cerda y Janneth Morocho 


Solución. Ponemos A 


7-16 
-10 4 -12 

-2 1 -1 


Se tiene que los valores propios de 


la matriz A son Ai = 2 , A 2 = 3, A 3 = 5 (realizar los cálculos). Los vectores propios 
asociados a Ai. A 2 y A 3 son los vectores oq = ( 1 , — 1 , — 1 ) T , a 2 = ( 1 ,— 2 ,— 1) T y 
«3 = (3, — 6 , — 2) t (recordamos que estos vectores no son únicos; así, el lector puede 
tener otros vectores propios. Lo importante es que estén calculados de forma correcta). 


Por los teoremas 3.3, 3.4 y el coroloario 3.1, la solución general del sistema es: 


fit) = 


— / ' CiOiiC 
i =1 


A it 


= Cl (l, -1, — l) T e 2í + c 2 (l, -2, —lf e 3í + c 3 (3, -6, -2)' T e 5t 
= Cl (e 2í , —e 2í , —e 2í ) T + c 2 (e 3í , -2e 3í , -e 3í ) T + c 3 (3e 5í , -6e 5í , -2e 5í ) T 


Observamos que otra manera de escribir la solución general del sistema es: 

X\(t) = Cie 2í + c 2 e 3í + 3c 3 e 5í , 
x 2 (t) = —C\e 2t - 2 c 2 e 3t - 6c 3 e 5í , 
x 3 (t) = —ce 2t — c 2 e 3< — 2 c 3 e 5í . 

Observación 3.13. El corolario 3.1 hace referencia a valores propios diferentes de 
la matriz A, pero estos pueden ser valores complejos. En tal caso, se tendrá que, 
si \j = a + bi es un valor propio de la matriz A para algún 1 < j < n, entonces 
Xj = a — bi también es un valor propio de la matriz A (recordamos que las raíces 
complejas de un polinomio siempre vienen en pares: una conjugada de la otra). 


Para el caso de valores propios complejos se puede demostrar que el vector 
solución de la ecuación diferencial vectorial, asociado al valor propio A = a + bi, tiene 
la forma <p(t) = (</?i(t), <p 2 (t),..., <p n (t)), donde las componentes de <p(t) están dadas 
por 

Pkft) = e at (ci(ak eos bt — /3k sen bt) + c 2 (ak sen bt + /3k eos bt)), (3.16) 

para 1 < k < n. 

Observación 3.14. En la ecuación 3.16 aparecen las constante arbitrarias ci,c 2 y 
las constantes akiPk Que se tiene que determinar (ver ejemplo 3.8). 
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Ejemplo 3.8. Hallar la solución general del siguiente sistema de ecuaciones 
diferenciales: 


i : - - -■ 


Solución. Para este sistema de ecuaciones diferenciales se tiene que A 


3 2 
-5 1 


Luego, los valores propios de la matriz A son Ai = 2 + 3i, X 2 = 2 — 3i (realizar 
los cálculos). Para hallar las constantes ai, «2 Y Pufo precedemos de la siguiente 
manera: suponemos que las soluciones del sistema son de la forma 5 : 


Vh(¿) = 2le Alí , 


i¡} 2 (t) = 03e Alí . 


Derivando y reemplazando en el sistema obtenemos: 

2lAie Alí = 32le Ali + 203e Al \ 

Q3Aie Alí = -52le Alí + 03e Alí . 

Reemplazando el valor de Ai, las últimas ecuaciones nos conducen al sistema (en 
las incógnitas 21, 03): 

{ (1 - 3¿)2l + 203 = 0 
-521 - (1 + 3¿)03 = 0 

1 3 

Se puede verificar fácilmente que 21= ly03 = — - — -i son soluciones del 
sistema anterior. Sustituyendo los valores de 21 y 03 en las soluciones ipi(t), 1 ^ 2 (í) 
tenemos: 

,/,,(() = e< 2+3i >‘, 

Mt) = (-i - ^P 2+3,) ‘. 

5 Estas soluciones son complejas pero, en el transcurso de los cálculos, se verá como se generan 
las soluciones reales del sistema de ecuaciones diferenciales. 
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Usando la fórmula de Euler, estas soluciones toman la forma: 

0i (t) = e 2t (eos 3 1 + i sen 3 t ), 

1 3 / i 3 

0 2 (t) = e 2t — eos 3 1 -\— sen 3 1 + i ( — sen 3 t -eos 3 t 

l ' J 2 2 V 2 2 


Se puede demostrar que las partes real e imaginaria de las soluciones complejas 
0i, 02 son también soluciones linea lm ente independientes del sistema. Luego, la 
combinación lineal de estas soluciones también es una solución del sistema 6 (ver 
observación 3.12). Así, las soluciones reales del sistema están dadas por: 


ipi (t) = e 2t (ci eos 3 t + c 2 sen 3 1) , 

Mt) = e 2t 


, 1 3 \ / 1 3 

ci ( — - eos 3í + - sen 3t ] + C 2 ( — - sen 3 1 — - eos 3 t 


Podemos concluir que a i = 1, «2 = — -, 0i = 0, 02 = —- 

Observación 3.15. La conclución del corolario 3.1 nos dice como es la solución 
general de la ecuación 3.15, todo esto bajo la condición de que todos los valores 
propios (reales y/o complejos) de la matriz A sean diferentes. Cuando los valores 
propios tienen multiplicidad no queda garantizado que las n soluciones del sistema 
fundamental de soluciones tengan la forma ae xt donde X es un valor propio múltiple 
y a es un vector propio asociado a X. 


Se puede demostrar que: si A es una matriz de orden n y A es un valor propio de 
multiplicidad m, donde 1 < m < n, entonces existen p vectores propios asociados al 
valor propio A donde 1 < p < m. Nosotros consideramos los dos casos: (1) p — m y 
(2) p < m. 

En el caso (1), p = m, existem m vectores propios linealmente independientes 
ai,...,a m asociados al valor propio A; luego, los vectores a\e xt ,..., a m e xt son 
soluciones linealmente independientes de la ecuación 3.15. 

Observación 3.16. Se puede demostrar que un tipo especial de ecuación diferencial 
vectorial que siempre lleva a este caso (p = m) es aquel en el cual la matriz A es 
simétrica, es decir, cuando la matriz A es tal que A = A T . 

6 La combinación lineal se debe hacer entre las partes real y la imaginaria de cada solución 
compleja. 
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Ahora explicamos muy someramente el caso (2), p < m. Para una explicación 
más detallada recomendamos la lectura de un texto más especializado de ecuaciones 
diferenciales. Por ejemplo, Braun [BR, cap. 3]. 

Sea A un valor propio de la matriz A cuya multiplicidad es m. Si p — 1, entonces 
existen m vectores linealmente independientes (que son solución de la ecuación 
diferencial vectorial 3.15) de la forma tpi(t) = a\e xt tal que (A — A/)aq = 0 y para 
todo j tal que 2 < j < m, 

( £ 771 — 1 £ 772—2 \ 

ai (j _ !)! + a2 (j -2)! + ' ‘ ‘ + aj ) eM ’ 

tal que (A — A I)otj-i = aj. 

Observación 3.17. Para el caso 2 < p < m se procede de manera similar, pero no 
se trata en este texto introductorio. 


Para finalizar esta sección, consideramos el siguiente problema: 


Problema 3.2. Hallar la solución del sistema de ecuaciones diferenciales 


dx{t 

dt 


= Aifr tal que x(t 0 ) = xq. 


Notamos que el problema 3.2 tiene la misma estructura que el problema con 

condiciones iniciales (o condiciones de frontera) formado por la ecuación diferencial 
dx 

— = ax y la condición inicial (condición en la frontera) x(to) = Xo cuya solución es 
x(t) = Xoe a( ' t ~ t °' ) (ver la sección 1.2 del capítulo 1, específicamente la definición 1.6 de 
la pag. 34 y los párrafos que siguen a esta definicón). Luego, por analogía, la solución 
del problema 3.2 tendría que ser el vector x(t) definido por la siguiente ecuación: 



4e (í " ío)A . 


(3.17) 


La ecuación 3.17 contiene un término de la forma e 21 , donde 21 es una matriz 
cuadrada. El término e 2t se llama exponencial de una matriz y es una función que 
generaliza la función exponencial (real o compleja) usual. Es decir, es una función 
que satisface las siguientes propiedades (como mínimo): 

1. Para todo s,t 6 C se cumple e s ^e m = e^' s+ ^ 2t . 
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2 . 

3. 


e° = /, donde o es la matriz nula, / es la matriz identidad. 

4- (e a ) = te m . 

dt v ’ 


En la literatura matemÁjtica se puede hallar la siguiente definición para la exponencial 
de una matriz. 


Definición 3.16 (Exponencial de una matriz). Sea 21 una matriz cuadrada con 
entradas complejas, definimos la exponencial de 21 como la matriz 

°° oi k i i 

^ = EF = í + a + S a> + 5¡ a,+ ■ 

k=0 

Se puede demostrar que la definición 3.16 garantiza el cumplimiento de las tres 
propiedades exigidas a la exponencial de una matriz. Para el lector interesado en este 
tema le recomendamos Klain [KLD] y Apóstol [AT|. 

Ej ercicios 


Resolver los siguientes sistemas de ecuaciones diferenciales: 


1 . 


2 . 


3. 


' dx\ 
dt 
dx 2 
dt 
dx 3 
dt 

' dx i 
dt 
dx 2 
dt 
dx 3 
dt 

dx i 
dt 
dx 2 
dt 
dx 3 
dt 


= Xi + x 2 - x 3 , 

= 2xi + 3x 2 - 4x 3 , 

= 4xi + x 2 - 4x 3 . 

= Xl — x 2 — x 3 , 

— X 1 + 3^2 + X 3 , 

= —3xi — 6 x 2 + 6 x 3 . 

— X\ + 2x2 + 2 x 3 , 

= 2xi + 3 x 2 , 

= 2xi + 3 x 2 . 
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4. < 


dx i 
dt 
dx 2 
dt 
dx 3 
dt 


x\ + 9x 2 + 9 x 3 , 
19x 2 + 18x 3 , 
9x 2 + 10x 3 . 
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La teoría de las ecuaciones diferenciales es una aplicación directa 
del cálculo diferencial e integral. El cálculo diferencial es una herra¬ 
mienta poderosa para la construcción de modelos matemáticos de 
problemas en los cuales están involucrados la variación de una o 
varias variables, las dependientes, con respecto a otra u otras 
variables, las independientes. 

Determinar el crecimiento poblacional, calcular el decaimiento 
radioactivo, hallar la edad de un fósil son problemas en los cuales se 
utiliza el cálculo diferencial para el planteamiento de un modelo 
matemático. 

El modelo matemático que se construye para hallar la solución de 
un problema particular deja planteado un nuevo problema: ¿qué 
herramientas se deben utilizar para hallar la solución del modelo 
matemático? Esta herramienta son las ecuaciones diferenciales 
ordinarias. 


Leónidas Cerda Romero. Profesional con más de 20 años de 
experiencia en docencia, doctor en Matemáticas por la Universidad 
de Concepción, Chile. En la actualidad su área de investigación es la 
teoría de números y su relación con la lógica matemática. 

Janneth Morocho Yaucán. Profesional con más de 20 años de 
experiencia en docencia, con una maestría en Estadística Aplicada 
de la Universidad de Granada, España. En la actualidad su investi¬ 
gación está encaminada en cuantiftcar el impacto de la ¡mplemen- 
tación de nuevas metodologías de enseñanza-aprendizaje. 



































